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Preparation and Uses of Silica Replicas in Electron Microscopy 


CHARLES H. GEROULD 
The Dow Chemical Company, Midland, Michigan 


(Received November 21, 1946) 


A detailed description of the preparation of silica replicas and substrates and their many 
and varied uses is presented with photographs and electron micrographs. A method of pre- 
paring silica replicas of specimens which cannot be subjected to the temperatures and pressures 
of the ordinary technique is described. It consists of applying a polystyrene lacquer to the 
specimen surface in place of the conventional molding. The question of “up and down”’ in 
stereo pictures of silica replicas is definitely settled by adding a metal oxide smoke (MgO) to the 
polystyrene replica before the silica evaporation. A successful powder dispersing technique con- 
sists of dispersing the powder in an ethyl cellulose lacquer and depositing am a silica sub- 
strate. Several other techniques are also included. 





INTRODUCTION 


ITHIN the last two years the field of 

electron microscopy has expanded many 
fold resulting in a great need for the development 
of new and improved electron microscopic tech- 
niques as well as the optimum use of those 
available. The great diversity of the recent 
entrants into the microscopic field has made an 
even greater need for literature on techniques 
which are of general application. 

The polystyrene-silica replica technique of 
Heidenreich and Peck' is one which has succeeded 
in solidly establishing itself and is now being 
widely used on many and varied research prob- 
lems involving the electron microscope. Among 
these we might mention the application of silica 
replicas to such metallurgical specimens as steels, 
copper alloys, and magnesium alloys.2-> We 


1R. D. Heidenreich and V. G. Peck, J. App. Phys. 14, 
23 (1943). 

?C. S. Barrett, J. App. Phys. 15, 691 (1944). 

7R. D. Heidenreich, SAE J. (trans) 53, 588 (1945). 

*R. D. Heidenreich, L. Sturkey, and H. L. Woods, J. 
App. Phys. 17, 127 (1946). 


might mention also their application to non- 
metallic materials such as teeth,* bones, minerals, 
glasses, ceramics, resins, woods, and textile 
fibers.? In addition, we must include the use of 
structureless silica substrates’ as supports for 
many different types of electron microscope 
specimens such as bacteria, viruses, powders, 
emulsions, metallic oxide smokes, thin evapo- 
rated films, microtome sections, etc. 

These many and varied uses of silica replicas 
have resulted in a demand for more explicit 
information on their preparation and uses. It is 
the purpose of this paper to fill this need and at 
the same time make available to the field some 
of the more recently developed, and as yet 
unpublished, microscope techniques in which 
silica replicas are an essential part. 

These new techniques include a new method 
~ §R. D. Heidenreich, C. H. Gerould, and R. E. McNulty, 
Tech. Pub. 1979, Metals Tech. 13 (April, 1946). 

6 C. H. Gerould, J. Dent. Res. 23, 239 (1944). 


7R. B. Barnes, C. J. Burton, and R. G, Scott, J. App. 
Phys. 16, 730 (1945). 


®R. D. Heidenreich, J. Opt. Soc. Am. 35, 139 (1945). 
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Fic. la. Block of clear polystyrene containing two 
specimens of a magnesium alloy sheet whose surfaces have 
been covered with a corrosion resistant coating. Silica 
replicas are to be prepared of this coating. The molding 
was formed from molding granules at a _— of 
160°C and a pressure of about 1500 'b./sq. 





Fic. 1b. Polystyrene impression on left and coated metal 
specimen on right. Mechanical separation has been 
accomplished by sawing away the excess polystyrene from 
the edges of the specimen and then pulling the two surfaces 
apart. Specimens may be dissolved with acids or alkalies 
if necessary. 


of preparing replicas for stereoscopy which makes 
it possible to settle definitely the question of “‘up 
and down” in stereo pictures, a procedure for 
preparing silica substrates, an ethyl cellulose 
lacquer particle dispersing technique for use with 
such substrates, a method of preparing silica 
replicas utilizing a room temperature produced 
polystyrene impression and several other miscel- 
laneous techniques. 
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Fic. 1c. Polystyrene impression of specimen mounted 7 


to 10 cm above conical tungsten filament containing 
approximately 1 mg of silica in the form of quartz. 





*. Fic. 1d. Actual evaporation of silica with filament 
incandescent under a vacuum of about 10-* mm of Hg. 
Silica is completelyfevaporated in 5 to 15 seconds, part of 
which condenses on surface of polystyrene to give a 
continuous film of approximately 200A thickness. 


All micrographs were obtained with either 
the RCA type B or U electron microscope. 


PREPARATION OF SILICA REPLICAS 


The cleanliness of the specimen surface is one 
of the most important points to be considered in 
preparing silica replicas. In the case of an etched 
metal the prepared surface must be free of 
reaction products, grease, stains, oxide, or hy- 
droxide films or any material which may either 
adhere to the replica or prevent it from repro- 
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ducing the true metal surface. When working 
with non-metallic specimens, the surface must 
likewise be free of any foreign material which 
would contaminate the replica. It is our belief 
that many of the difficulties met in preparing 
silica replicas by the polystyrene-silica technique 
arises from improperly prepared surfaces, mainly 
dirty surfaces. 

The polystyrene-silica technique can best be 
illustrated by a series of photographs of the 
various steps involved in preparing the poly- 
styrene impression, evaporating the silica upon 
the plastic, and removing the silica replica film. 
These various steps are described in detail in 
order that one may experimentally follow the 
technique and obtain good replicas with a 
minimum of difficulty. 

The first step, after the specimen has been 
given a suitable surface preparation, is to mold 
polystyrene against the surface. This is accom- 
plished quite simply by imbedding the specimen 
in polystyrene as illustrated in Fig. 1a. Poly- 
styrene molding granules of medium viscosity 
or molecular weight* are best for this. 

The molding may be accomplished in an ordi- 
nary electrically heated metallographic specimen 
mounting press or a steam-heated compression 
molding die. The specimen is placed in the 
molding die with the prepared surface upward 
and sufficient molding granules added to give a 
final clear molding at least } inch thicker than 
the specimen. The actual molding is accom- 
plished by heating the molding granules to a 
temperature of 160°C, care being taken that no 
pressure is applied until they are thoroughly 
softened. Otherwise they may be pressed against 
the specimen surface causing deformation. After 
the granules have had a chance to soften thor- 
oughly (5 to 10 minutes usually is required), 
pressure is applied up to about 1500 lb./sq. in. 
and the cooling of the die commenced at once. 
The pressure is maintained until the molding has 
cooled below 80°C since release of the pressure 
at higher temperatures may result in the forma- 
tion of bubbles in the molding. 

The specimen is now separated from the 
plastic impression by first sawing away the 
excess material from the edges of the specimen 


* Such a material is available commercially as Styron 
No. 379-clear. 


VOLUME 18, APRIL, 1947 


and then separating the two, mechanically if 
possible. Figure 1b shows a magnesium alloy 
specimen which has just been mechanically 
separated or stripped from the polystyrene. 

If mechanical separation is impossible without 
fracturing or otherwise distorting the polysty- 
rene, chemical reagents may be used to dissolve 
the specimen. Mineral acids or strong caustic 
solutions may be used on most specimens without 
damage to the plastic impression. 

Figures 1c and 1d show the next step in which 
silica is evaporated in vacuum upon the surface 
of the plastic impression. The molding is mounted 
7 to 10 cm above a spiral tungsten filament 
(cone shaped of about 5-mm length with a 35° 
angle) containing about 1 mg of clear quartz in 
the form of small slivers weighing roughly 0.2 mg 
each. The filament of the type shown in the 
photographs is heated to incandescence for a 
period of 5 to 15 seconds, using 25 to 30 amp. 
at 10 to 15 volts. The silica is completely evapo- 
rated in a few seconds, providing the individual 
splinters of silica are not too large and the 
vacuum is 10-* mm of Hg or better. 

The tungsten filament used in this evaporation 
is prepared by winding 20-mil tungsten at a dull 
red heat on a tapered mandrel having about a 
35° included angle. The quartz used may be 
prepared quite easily by breaking up a natural 
clear quartz crystal by means of mortar and 
pestle and then screening to obtain particles of 
the proper size. Other laboratories have success- 
fully used cloudy quartz, fused silica, or Santocel. 

The block of polystyrene supporting the silica 
film is now prepared for the removal of the film. 
The sides and bottom are first scraped thoroughly 
to remove all undesirable silica film and ‘any dirt 
which may contaminate the solvents. The silica 
film is scored into slightly over }’/’ squares and 
the polystyrene block lowered to the bottom 
of a dish of high purity ethyl bromide and 
benzene with the replica side up as shown in 
Fig. 2a. The block will stick to the bottom of the 
dish as the plastic softens. About 10 percent of 
benzene is added to the ethyl bromide to decrease 
its volatility and thus prevent moisture conden- 
sation upon the silica films as they are removed 
from the solvent. Such moisture ‘condensation is 
invariably accompanied by precipitation of the 
polystyrene. Ethyl bromide is used for this step 
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Fic. 2a. Silica film has been scored into squares and the 
whole polystyrene block is about to be immersed in ethyl 
bromide plus 10 percent benzene. 





Fic. 2b. Silica films are seen as they start floating off 
from the polystyrene block stuck to the bottom of the 
dish of ethyl bromide and benzene. Silica films generally 
come free from the polystyrene 1 to 5 min. after immersion 
in the solvents. 


in the replica technique for the following reasons: 
1. It is an excellent solvent for polystyrene; 
2. its density is sufficiently high to cause the 
films to float about in the solvent; 3. its great 
volatility tends to produce convection currents 
‘which aid in washing the films free from poly- 
styrene; 4. its refractive index is such that the 
silica films are generally readily visible in the 
solvent. 

The silica films will generally free themselves 
from the polystyrene 1 to 5 minutes after im- 
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Fic. 2c. Silica films are being transferred on a piece of 
200-mesh stainless steel screen to a rinse consisting of 
ethyl bromide plus a few percent of iso-amyl alcohol. 





Fic. 2d. Catching the silica films on a standard electron 
microscope specimen screen after they have washed in the 
solvent rinse for at least 5 min. After catching, the screens 
are blotted on a filter paper and viewed directly in the 
electron microscope. 


mersion. Side illumination with a dark red back- 
ground is of great help in seeing them. Figure 2b 
shows the films as they come free from the 
polystyrene. They are transferred almost at once 
to a rinse containing ethyl bromide plus a few 
percent of iso-amyl alcohol added to reduce the 
solvent’s volatility. Benzene is not used because 
an effective concentration is probably upwards 
of 10 percent. Transfer of the silica films to the 
rinse is easily accomplished by using a piece of 
200-mesh stainless steel screen as shown in Fig. 
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Fic. 3. Silica replica of ‘magnesium alloy showing twinning _ F'16- 6. Silica replica of cleaved enamel from human tooth, 
and fine globular precipitation, 4000. 3500 X. 





Fic. 4. Silica replica of ee eee of quartz crystal, Fic. 7. ta wie ¢ internal 





Fic. 5. Silica replica of glass etched with HF, 10,500. Fic. 8. Silica replica of maple wood cleavage, 4500 X. 
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Fic. 9. Comparison between silica replicas of a magnesium alloy prepared by 
means of: (a) a polystyrene molding using both heat and pressure, 6000. 
(b) a polystyrene lacquer film using no heat or pressure, 6000. 


2c. The films readily loose themselves from the 
screen when it is lowered into the rinse provided 
they are not allowed to dry during transfer. 

The silica films are allowed to wash at least 
five minutes and then caught on the standard 
microscope specimen screens as shown in Fig. 2d. 
The screen with the films are dried by blotting 
on a filter paper. Such films are believed to have 
a thickness of about 200A. 

Figures 3-8 show examples of silica replicas as 
applied to a magnesium alloy, quartz crystal, 
etched glass, tooth, bone, and wood. 

One often finds in replica work a need for 
obtaining silica films from specific areas of a 
specimen. This is not difficult since the poly- 
styrene molding with the silica film may be 
examined quite readily with the light microscope 
at magnifications up to 1000 X using illumination 
through the specimen. Regions which are of 
particular interest may be outlined readily by 
scribing the surface with a needle point and 
then scraping off all the silica film except that 
of particular interest. This film may then be 
removed in the usual manner. Such a technique 
is often found very useful in studying interfaces 
between dissimilar materials. 

As has been mentioned before, many of the 
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- structure in the original specimen. It appears that 


Courtesy A.I.M.E, 


difficulties experienced when preparing silica 
replicas arise from improperly prepared surfaces. 
Dirty surfaces often give silica replica films 
which tend to roll or fold up and which may 
refuse to straighten out by the customary method 
of ‘‘teasing’’ the films with a specimen screen as 
they float in the rinse. Another difficulty often 
experienced in preparing silica replicas is fragility 
of the films caused by lack of appreciable fine 


fine detail in a specimen adds considerably to the 
mechanical strength of silica replicas. 

On the other hand, silica replicas of single 
phase or electrolytically polished metals, glass, 
or other microscopically smooth surfaces are very 
often extremely fragile. Care must be taken to 
handle these films delicately and even then part 
of the film may break up on the specimen screen. 

Occasionally one wishes to prepare replicas of 
surfaces which for some reason cannot be sub- 
jected to the temperatures and pressures neces- 
sary to obtain a polystyrene impression. In such 
cases an alternate method is used in which a 
polystyrene lacquer is applied to the specimen 
to give a heavy coating, stripped off and then 
used to prepare silica films in the ordinary 
manner. This technique has given quite reliable 
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results in practice using a 1 percent solution of 
polystyrene in benzene for the initial application 
and then backing this film up with several 
applications of a more viscous lacquer of 2 to 3 
percent polystyrene in benzene. By means of the 
more concentrated solution one can build the 
original film up to the point where it may be 
mechanically stripped from the specimen or the 
specimen dissolved if it is felt to be advisable. 
Figures 9a and b show a comparison between the 
conventional method and this polystyrene lac- 
quer technique. 

The resolution possible with this technique 
appears to be very good though very likely not 
equal to that of the pressure molding technique 
whose resolution is apparently better than 100A. 
The resolution of this ‘‘cold’’ replica technique 
can probably be improved, as suggested by 
Matheson,* by vacuum boiling of the specimen 
in the lacquer as a means of outgassing the 
surface. The initial lacquer coating can then be 
built up as previously described. 


STEREOSCOPY 


The use of stereoscopy as an electron micro- 
scope technique has become of quite general 


application and is held to be invaluable in inter- 
preting electron micrographs by many micro- 
scopists. The stereoscopic technique is especially 
well suited for use with silica replicas. However, 
one of the difficulties experienced in such applica- 
tion has been the impossibility of definitely 
establishing the “‘up and down”’ of the stereo 
pictures unless something was known about the 
topography of the original specimen surface. 

A very simple procedure is now being used to 
remove this uncertainty. This consists of adding 
a metal oxide smoke to a known side of the 
silica replica film and using this oxide smoke to 


determine the correct viewing arrangement. This 


is best accomplished by depositing the oxide 
smoke on the polystyrene impression of the 
specimen surface. Magnesium oxide smoke works 
very successfully, an even distribution being 
obtained by wafting the molding through the 
smoke of burning magnesium. Zinc oxide smoke 
can also be used but is much more difficult to 
distribute evenly. 

Silica is evaporated upon the polystyrene 
impression and the films removed in the usual 
manner carrying the oxide smoke with them, 





Fic. 10. Stereo pair from a fine grain magnesium alloy showing the addition of magnesium oxide 
smoke (indicated by arrows) to the replica film. The question of “up and down” is answered by 
viewing the stereo pair arranged so the oxide cubes appear below the surface. Compound precipi- 
tation can be observed projecting upward in relief from the grain boundaries. Magnification 6000 X. 


*L. A. Matheson and R. D. Heidenreich, ‘Magnification calibration of electron microscopes,’ Presented at Chicago 


meeting of E.M.S.A. (November, 1944). 


*R. D. Heidenreich and L. A. Matheson, J. App. Phys. 15, 423 (1944). 
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Fic. 11. Silica smoke dispersed in ethyl cellulose on a 
silica substrafe, 7500. 


attached to the side reproducing the specimen’s 
surface structure. Electron micrograph stereo 
pairs of these films are then obtained and viewed 
through the stereoscope with the pictures ar- 
ranged so the oxide smoke appears below the 
surface, the proper arrangement. By this method 
the question of “up and down” is definitely 
settled. Figure 10 shows a stereo pair of a 
magnesium surface with MgO smoke on the 
replica. 


SILICA SUBSTRATES AND A POWDER 
DISPERSING TECHNIQUE 


Structureless silica substrates have been found 
very useful as supports for a variety of electron 
microscope specimens. They can be prepared 
having a thickness of from <100 to 200A and 
are completely structureless and without holes 
or fibers. The films are structurally strong as 
compared to organic films, and may be subjected 
to very intense electron beams with a minimum 
of electron scattering, “charging up,” and break- 
age. Microscope specimens may be heat treated 
at elevated temperatures for short periods of 
time while supported on silica substrates. They 

_also can be successfully used as the support for 
specimens for electron diffraction transmission 
study. 

These substrates are prepared by dipping a 
clean and dry microscope slide in a 2—5 percent 
solution of polystyrene in benzene and allowing 
to dry. The silica film is formed upon the lacquer 
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Fic. 12. CaCO; dispersed in ethyl cellulose on a silica 
substrate, 4500. 


film by the same technique as used for replicas 
except that 10 or 20 percent less silica is generally 
used. The films are removed in the same manner 
as other silica replica films and caught on stand- 
ard electron microscope specimen screens. How- 
ever, the silica films separate from the thin 
polystyrene film quite slowly, } hour to several 
hours usually being required. 

Difficulties are occasionally met in preparing 
these films and they seem to be largely caused by 
improper thickness of the polystyrene film. For 
some inexplicable reason, if this film is too thick, 
the silica films are fragile, and if too thin, the 
silica films come loose from the lacquered slide 
with great difficulty or may not come loose at all. 





Fic. 13. Polystyrene emulsion dispersed on a silica 
substrate, 7500X. 
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Fic. 14. Magnesium-aluminum alloy evaporated in vacuum 
and condensed on a silica substrate, 6500 X. 


It appears that a 2-5 percent lacquer solution 
gives a film of optimum quality. 

Precautions must also be taken to see that the 
lacquer is relatively moisture-free and the glass 
slide dry. Failure to observe these precautions 
may result in a cloudy lacquer film and sub- 
strates which show bubbles. 

Many different particle dispersing techniques 
have been published'®" which make it possible 
to determine accurately -particle size and shape. 
Most of these techniques can be successfully used 





Fic. 15. Silica replica of a polished steel surface shadowed 
with chromium, 3500. 


”M. L. Fuller, P. G. Brubaker, and R. W. Berger, J. 
App. Phys. 15, 201 (1944). 

11 A, M. Cravath, A. E. Smith, J. R. Vinograd, and J. 
N. Wilson, J. App. Phys. 17, 309 (1946). 
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Fic. 16. Lampblack dispersed in ethyl cellulose on a silica 
substrate and shadowed with chromium, 3500. 


with silica substrates with improved results. A 
powder dispersing technique has been in use here 
for several years which has given excellent results 
on a wide variety of powdered materials ranging 
in size from <0.01 micron to several microns. 
This technique consists of dispersing the powder 
in an ethyl cellulose lacquer as a thick slurry, 
diluting and adding to the silica substrate. This 
is accomplished in practice by placing a few 
milligrams of the powder in a mortar and adding 
a drop or two or an ethyl cellulose lacquer 
consisting of 2 percent ethyl cellulose in 70 
percent amyl acetate and 30 percent ethyl alcohol. 
The mixture is lightly ground with mortar and 
pestle to a thick slurry or paste during which 
time the dispersion of the particles actually takes 
place. Amyl acetate is then added slowly at first 
and the concentrated dispersion diluted to the 
proper concentration for adding to the silica 
substrate on a specimen screen. Tests have indi- 
cated that the effect of the grinding action upon 
the particle size and shape is negligible in speci- 
mens having a particle size of less than 5 microns. 

Ability to judge the proper quantities and 
dilutions is readily obtained with a little experi- 
ence. Generally the slurry dilution is about 1:50. 
Quantities should be varied so that the minimum 
thickness of dried ethyl cellulose is used, i.e., 
just sufficient to keep the particles firmly at- 
tached to the silica film. Figures 11 and 12 show 
typical particle dispersions obtained by this 
technique. 
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Courtesy A.I.M.E. 


Fic. 17. Silica replica of a magnesium-nickel alloy 
showing not only a replica of the globular Mg2Ni particles 
but also some of the actual particles themselves on the 
silica replica film. The faithfulness of the technique in 
reproducing the size and shape of the particles is evident. 
Figure 3 also shows this same’ phenomenon on a smaller 
scale. Magnification 4000 


Silica substrates may also be used for sup- 
porting a variety of other types of specimens as 
illustrated by Figs. 13 and 14 showing an emul- 
sion and an evaporated magnesium alloy. Emul- 
sions are readily dispersed by diluting to the 
proper concentration and then adding a drop of 
the liquid to the substrate on a specimen screen. 


MISCELLANEOUS USES AND TECHNIQUES 


Silica replicas and substrates have also been 
found to be useful in the field of heavy metal 
“‘shadow-casting.”"* This technique is probably 
not of great value on silica replicas for the 
purpose of improving contrast since silica replicas 
normally show considerable contrast. However, 
the technique has definite possibilities as a means 
of making elevation determinations without the 
- use of stereoscopy. 

Shadowing of silica replica films is probably 
best accomplished by evaporating the heavy 
metal upon the silica film while it is on the 
surface of the polystyrene molding. The films are 
then removed in the usual manner. Figure 15 
shows an example of the “‘shadow-casting”’ of a 
polished steel specimen. 





® R. C. Williams and R. W. G. Wyckoff, J. App. Phys. 
15, 712 (1944), 
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Shadowing of various electron microscope 
specimens supported on silica substrates is of 
definite value both from the point of view of 
improving contrast and of making elevation 
determinations. This technique is probably best 
accomplished by applying the polystyrene lac- 
quer film to a glass slide, evaporating the silica, 
adding the specimen, “shadow-casting,”’ and 
then removing the silica film in the usual manner. 
It seems advisable to use such a procedure rather 
than attempt the shadowing of the substrate 
and specimen while supported by the specimen 
screen. In this case the silica film may be uneven, 
making it a very poor target for accurate 
“‘shadow-casting.”” Figure 16 shows a specimen 
of lampblack ‘‘shadow-cast’”” in this above 
recommended manner. 

A very interesting and often valuable use of 
silica replicas lies in the isolating of impurities 
and certain second phase elements or compounds 
from some metals. This occurs occasionally in 
metals which contain one or more elements or 
compounds which are present as a second phase 
and which are relatively insoluble in the etchant 
used during the surface preparation. Such con- 
stituents of the metal are often either left lying 
on the metal surface after etching or are pro- 
truding from the surface and only loosely at- 
tached to the metal. When polystyrene is molded 
against the metal surface and separated, these 
particular constituents are imbedded in the poly- 
styrene and removed with it. When the silica 
evaporation is completed and the polystyrene 
dissolved, as the silica films are removed, the 
particles are often transferred to the silica replica 
and attached to it. As a result we find that in 
the electron microscope we see not only a replica 
of certain second phase constituents in the metal 
but also some of these particular constituents 
actually attached to the replica film. Figures 3 
and 17 show examples of this phenomenon. Such 
materials can often be identified by transmission 
electron diffraction. 

Occasionally one finds an electron microscope 
specimen in particle form which tends to sublime 
when subjected to the high vacuum and heat of 
electron bombardment in the microscope. This 
sublimation may be appreciably decreased by 
supporting the specimen on a silica substrate 
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and then evaporating more silica over the surface 
of the specimen. Sulfur is a material which may 
be treated in this manner. 


CONCLUSIONS 


1. Silica replicas of a variety of bulk materials 
such as metals, minerals, glasses, ceramics, resins, 
woods, etc., can be successfully prepared using a 
heat and pressure molded polystyrene impression 
plus evaporated silica. 

2. Silica replicas can be successfully prepared 
from bulk materials which cannot tolerate either 
the temperatures or pressures of (1), by applying 
a heavy polystyrene lacquer, removing and 
evaporating silica upon the lacquer impression. 

3. The question of ‘‘up and down” in stereo 
pictures of silica replicas can definitely be estab- 
lished by depositing a metal oxide smoke on the 
polystyrene impression before evaporation of the 
silica. 

4. Structureless silica substrates are well suited 
as the supporting membranes for specimens such 


as powders, emulsions, metallic oxide smokes, 
thin evaporated films, bacteria, viruses, micro- 
tome sections, etc. 


5. Excellent dispersions of powders from <0.01 


to several microns in diameter may be obtained 
on silica substrates using an ethyl cellulose 
lacquer dispersing technique. 


6. Heavy metal “‘shadow-casting” can be 


effectively applied to both silica replicas and 
silica substrates. 


7. Isolation of metallic phases directly on 
silica replicas of metal surfaces is often useful in 
minutely examining and identifying these con- 
stituents. 
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The X-Ray Storage Properties of the Infra-Red Storage Phosphor 
and Application to Radiography 


O. E. BerG AnD H. F. KalIser 
Naval Research Laboratory, Washington, D. C. 


(Received November 4, 1946) 


An application of the characteristics of infra-red storage phosphors to radiography is de- 
scribed, whereby the expense and trouble of x-ray film processing may be partially or com- 
pletely eliminated. The properties of these materials to store x-ray energy for periods of time 
and release that energy upon exposure to infra-red radiation have proved practical for either 
temporary and permanent recording of the latent x-ray image. One distinct feature of this 
method is the complete removal of harmful x-radiation to the technician. The possibility of a 
new field of radiography in ‘‘flash fluoroscopy”’ is also suggested by the authors. Graphs and 
charts exhibiting the characteristic behavior of infra-red phosphors to x-ray energy are included 
as well as actual photographs and radiographs comparing radiography with phosphorography. 





HE image storage phosphor was a research 

development required in war work as a 
practical and sensitive means for the detection 
of infra-red radiation. In the course of research! 
carried on to produce such phosphors it was found 
possible to produce phosphors which exhibited 
storage properties not only for visible light, but 


1Contract OEMSR 982, OSRD Report 16.5—119 
(N5643). 
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also for cathode rays, and x-rays. The optical 
characteristics of such materials have been 
described in papers of Roland Ward? and of N. 
F. Miller and C. E. Barnett.* 


2R. Ward, “Preparation and properties of infra-red 
sensitive strontium selenide and sulfide-selenide phosphor,” 
J. Opt. Soc. Am. 36, 351(A) (1946). 

3N. F. Miller and C. E. Barnett, “Infra-red sensitive 
phosphors from zinc and cadmium sulfides,” J. Opt. Soc. 
Am. 36, 352(A) (1946). 
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Since no attempt had been made to study 
these materials for x-ray energy storage char- 
acteristics, the authors have carried out a series 
of experiments on materials kindly furnished by 
the Brooklyn Polytechnic for this purpose. The 
ability of these materials to save x-ray energy 
suggests their application to radiography as a 
possible substitute for the photographic emulsion 
for the retention of a latent image. 

The ability of these materials to store radiant 
energy is not surprising in view of the fact that 
all phosphors exhibit such properties when cooled 
to liquid air temperature, exposed to an energy 
source to fill the electron traps and then warmed 
up, whereupon the freeing of electrons from the 
traps is indicated by bursts of light.‘ 

Phosphors may be produced which store 
energy from x-rays, beta-rays, gamma-rays, and 
visible light. and which hold this energy for 
indefinite periods of time. The energy is released, 
in the form of light varying from red to ultra- 
violet, by infra-red radiation with speed of release 
governed by intensity of the infra-red radiation. 
Figure 1A and B shows how these effects may 
be employed in radiography. The x-rays pass 
* through the subject into the phosphor plate or 
screen, and by differential intensity record in 
detail thereon producing a latent image. This 


* “Characteristics of luminescent materials for cathode- 
ray tubes,”’ J. Inst. Elec. Eng. 92, Part III, No. 20, 300 
(December, 1945). 
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latent image is quite stable for many phosphors. 

Then, in an adjoining room, away from harm- 
ful x-radiation, and at the convenience of the 
operator, the image is made visible to the eye by 
irradiating the plate with infra-red radiation. If 
a particular area of the image is of special 
interest to the technician, he may expose it 
separately with an infra-red spotlight (photo- 
graphic masking technique is also applicable 
here). The length of viewing time of a phosphor 
plate may be determined from its energy decay 
rate graph shown in Fig. 2A which indicates 
the intensity of the released light as a function 
of time and intensity of the infra-red source. 

When plotting light units against log time as 
in Fig. 2B, we get several linear portions indi- 
cating that the decay of luminescence is not of 
simple exponential character, but occurs in a set 
of several such processes. Light values are given 
as arbitrary light units measured by a 931A 
photo-cell multiplier with voltage regulators. 

If a permanent record of the x-ray image is 
desired, it is an inexpensive and a simple task to 
expose the phosphor image directly to a sensitized 
paper (sensitive in the color region of the partic- 
ular phosphor in use). Such a process eliminates 
the need for double emulsion x-ray film. 

A 4”X5” experimental phosphor viewing plate 
was prepared according to the following pro- 
cedure: 
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While it might have been desirable to apply 
the phosphor to a thin cardboard base as for a 
fluorescent screen, because of the high processing 
temperatures involved, it was found necessary 
and advantageous to use a 33”" steel plate with a 
thick silver plating. This base plate was not 
affected by high temperatures and its silver 
coating prevented chipping or peeling of the 
phosphor during experimentation. 

Furthermore, silver plated steel is preferable 
to a paper base in that it serves as an absorbing 
medium for transmitted x-radiation after the 
x-rays have transversed the subject and the 
phosphor, thereby reducing scatter. 

An emulsified mixture of phosphor, ethylene 
dichloride, and jasonite (plastic binder) was 
applied with a small artist’s spray gun, and dried 
at 850°C to remove excess binder and render the 
phosphor active.5 After a few experiments, it 
was found that a .020-inch coating of the 
phosphor was sufficiently thick to produce peak 
luminescence. 

Test exposures were made on the resultant 
plate at various kilovoltages and exposures. 


5 For further information on phosphor coating applica- 
tion see OSRD Report No. 5643. 
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Figures 3 and 4 are comparative results of an 
x-ray film image and a phosphor image, respec- 


tively, of a 3-inch ammeter. It is an indication 
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meter, X-ray exposure made 
on fine grain industrial x-ray 
film. 
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Fic. 3. Radidgram of 3 inch 
ing. 


the definition and latitude obtainable in 


“‘ phosphorogram.” 
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An interesting and advantageous feature con- 
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Fic. 4. Phosphorogram of 
3-inch meter. Dark areas in- 
dicate uneven phosphor coat- 


cerning latitude in this method of exposure is 
revealed in Fig. 5, another photograph of the 
same exposure after some diminution by previous 





Fic. 5. Secondary film expo- 
sure to phosphor image of me- 
ter. Notice increased latitude 
and preservation of detail. 
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Fic. 6, The storage sensitivity curve in the x-ray region of the spectrum for one particular infra-red phosphor. 
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infra-red exposure. The definition is as sharp as 
before and areas that were overexposed (too 
brilliant) in the previous photographic exposure 
now appear in readable detail. The latitude of a 
phosphorogram is thus extensive and _ flexible. 
An overexposed plate can be irradiated with 
infra-red till it reaches a readable intensity. 
Likewise, an underexposed plate can be increased 
in intensity to a readable density by increasing 
the infra-red radiation. 

Figure 6 shows the storage sensitivity curve in 
the x-ray region of the spectrum for a particular 
phosphor. Each type of storage phosphor has its 
own sensitivity curve and exhibits its own char- 
acteristic properties as well. Although a blue- 
green phosphor was used extensively in these 
experiments, it is now possible to obtain a 
phosphor for any color region in the visible 
spectrum. Blue-green is a comfortable color to 
the eye, but it is not the most efficient color for 
definition. 

In general, the results of these experiments 
have revealed numerous possibilities of infra-red 
phosphors in radiography. Their properties, how- 
ever, are not always advantageous. 

Referring back to Figs. 3 and 4, we recognize 
the same undesirable characteristics of graininess 
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and diffusion as in the fluorescent screen. This 
might be corrected to a large extent by more 
refined processing methods. Figure 2 shows an- 
other disadvantage—a reluctant release of stored 
energy, because of the exponential decay rate, 
which may handicap the speed of repeated expo- 
sures upon a given phosphor plate. For all prac- 
tical purposes, however, the energy remaining in 
the phosphor after reasonable exposure to infra-red 
radiation is negligible and the contrast of the re- 
corded image appears to diminish with time. 

‘Phosphorography will not immediately replace 
the x-ray film nor the fluorescent screen. It will 
supplement these methods and be particularly 
valuable in special cases. For example: flash 
fluoroscopy will be possible with the high current 
condenser-discharge x-ray tubes, whereby an 
image made with a very short x-ray exposure of 
high intensity will be made immediately visible, 
and retained long enough for inspection or 
recording on paper or film. 

Grateful acknowledgment is extended to Drs. 
R. E. Ward and R. T. Ellickson of the Poly- 
technic Institute of Brooklyn for their willing 
cooperation in supplying information and ma- 
terial for this work. 
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A Flux Plotting Method for Obtaining Fields Satisfying Maxwell’s 
Equations, with Applications to the Magnetron 


Prescott D. Crout 
Radiation Laboratory,* Massachusetts Institute of Technology, Cambridge, Massachusetts 
(Received September 13, 1946) 


In this report the flux plotting methods which have previously been applied to fields satis- 
fying Laplace's and Poisson's equations are extended to fields satisfying Maxwell's equations. 
Application is made to the hole and slot, and vane types of magnetron tubes. 


1. INTRODUCTION 


66 [LUX plotting” is a procedure for drawing 

two mutually orthogonal families of 
curves in such a way that the rectangles so 
formed have a specified geometric shape at each 
point. This shape is so determined that the curves 
represent a field which satisfies a specified par- 
tial differential equation. This equation, which 
characterizes the field, is not used directly, but 
it or its equivalent is used only indirectly to de- 
termine the geometrical condition that is used 
in actually constructing the plot. Flux plotting 
procedures for obtaining fields satisfying La- 


9.8 cm Magnetron 


~_m 


fcathode) 


\ 4 


~~ 





Fic. 1. Cross section of hole and slot type magnetron. 
Magnetic field is axial (perpendicular to the plane of the 
paper). Its directions in the anode wells at a certain instant 
are shown by the circles @ and ©. Electric field is across 
the gaps between the anode segments. Its directions at 
a certain instant are shown by the arrows. 


* This paper is based on work done for the Office of 
Scientific Research and Development under contract 
OEMsr-262 with the Massachusetts Institute of Tech- 
nology. 
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place’s equation, Poisson’s equation, and certain 
other equations have been devised and applied 
successfully to wide varieties of problems.'! The 
purpose of this paper is to extend these methods 
to fields characterized by Maxwell’s equations. 
In so doing we do not use these equations di- 
rectly, but use Ampere’s law and Faraday’s law 
instead. These laws together with Stoke’s theo- 
rem evidently give Maxwell’s equations. Some 
familiarity with flux plotting methods is assumed. 

For the sake of clarity the new method will 
be described in connection with the specific ap- 
plication which gave rise to it—that of an in- 
finite, hole and slot type magnetron, a cross 
section of which is shown in Fig. 1, operating in 
its main mode (+ — + — or II mode). Later in 
Section 6 application is made to the vane type 
magnetron. All fields are sinusoidally alternating. 
The electric current which links each anode well 
(cavity) and passes across the gap between two 
adjacent anode segments produces the axial 
magnetic field in the well (Ampere’s law); also 
the rate of change of this axial magnetic field 
produces the voltage which maintains the elec- 
tric current (Faraday’s law). Because of sym- 
metry we need consider only the field which lies 
between two axial planes which bisect, respec- 
tively, an anode segment and an anode well, Fig. 2. 


1H. Poritsky, “Graphical field-plotting methods in 
engineering,” Trans. A.I.E.E. 57, 727 (1938). B. Hague, 
“Méthodes Analytiques, Graphiques et Expérimentales 
Utilisées Pour Etude des Champs Magnétiques et 
Electriques dans Les Machines et Appareils Electriques,”’ 
Comptes rendus 4, 47 (1932). (This reference contains an 
extensive bibliography.) A. D. Moore, Fundamentals of 
Electrical Design, p. 31. P. D. Crout, “‘The Determination 
of Fields Satisfying Laplace’s, Poisson’s, and Associated 
Equations by Flux Plotting,”’ Radiation Lab. Rep. 1047. 

2 The field so obtained can be made the basis of calcula- 
tions for determining the frequency spectrums of finite 
strapped or unstrapped magnetrons. See Crout and Both- 
well, “A Method for Calculating Magnetron Resonant 
Frequencies and Modes,” Radiation Lab. Rep. 1039. 
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2. NATURE AND CONSTRUCTION OF 
THE FLUX PLOT 


The peak electric field is shown in Fig. 2 by 
means of a number -of electric lines of force, 
which divide the region into corresponding tubes 
of flux; and a number of lines of no work, which 
are orthogonal to the lines of force, and which 
divide the region into corresponding strips. The 
peak magnetic field is perpendicular to the plane 
of the paper. Let the various tubes and strips be 
numbered consecutively as indicated, and let 


L;;=length of elementary rectangle common to ith 
tube and jth strip measured along a line of 
force through the center; 

W;;= width of elementary rectangle common to 7th 
tube and jth strip measured along a line of no 
work through the center; 

A,;=area of elementary rectangle common to ith 
tube and jth strip; 

oi@=electric flux in the ith tube, ¢ being the arbi- 
trarily chosen flux per unit tube; 

E=electric intensity =electric flux density since 
dielectric constant = 1; 

H=magnetic intensity=magnetic flux density 
since permeability = 1; 

w=angular frequency; 

c=velocity of light. 


(1) 


Gaussian units are used throughout 


We shall first consider the equilibrium of an 
elementary rectangle, and shall leave till later 
the consideration of the elementary triangles 
which occur at the “kernel’’ or point of inter- 
section of the lines of no work. The total peak 
displacement current flowing outside (on the 








Line ot force of electric 
field (also a constant H) 


Line of no work of 


electric field my 





Kernel 


Fic. 2. Electric field in hole and slot type magnetron. 


cathode side) of the center of rectangle 77 is* 
i-1 ow 
soi +> oO; j-—> (2) 
k=1 4 


hence, noting Ampere’s law, we see that the value 
of H at this point is independent of j, and is 


given by 
’ i—1 ow 
g0i+ y Ox j—- (3) 
k=l c 


The total magnetic flux in the jth strip between 
the center of rectangle ij and the kernel is 
therefore 


ow :; i—1 ow n—l l-1 ow n—1 
‘S6;4+>. CK A ijt =e » toit>, CK Aijt+ = ajon +d Ok Aaj, (4) 
2c k=1 C l=i+1 k=1 c k=l 


where 
Aji; = LijWi; for rectangle ij, 
L,j= base of triangle in jth strip, 
W,,;=altitude of triangle in jth strip, 
A,;=area of triangle in jth strip= 
(Ln; )average Wj; (5) 


a;=fraction of W,; from the last line of force to the 
center of gravity of the last triangle (equals 4 for 
triangles with straight sides). 


By using the center of gravity as the point at 
which to evaluate B in an elementary rectangle 
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or triangle, multiplication of this value of B by 
the elementary area gives the flux exactly if the 
surface which represents B is a plane, regardless 
of whether or not that plane is horizontal (B 
constant). The first term in (4) gives the peak 
magnetic flux in half of rectangle 7j7, whereas the 
last gives that in the triangle. 





3 From here on the word “peak” will be omitted, it being 
understood that all quantities are peak values of sinu- 
soidally varying quantities. Complex notation can be used 
if desired, in which case the magnetic intensity will be 
purely imaginary. 
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Noting Faraday’s law we see that the electric flux in the ith tube can be obtained from the rate of 
change of the magnetic flux given by (4), thus 


ow i—1 n—1 n—1 l-1 wW i; 
oo=|“Y3(sortE ou) At (aso t¥ on )Aart * (to+E 01)An|| (6) 
Cc kml k=l l=i+1 k=1 cLij 


Here the brace is the magnetic flux given by (4), which with w gives the rate of change of this flux, 
or the induced voltage in abvolts.* Division by cL;; gives the electric intensity in statvolts/cm, and 
multiplication by W,; gives the total electric flux in the rectangle, which must equal ¢;¢. Cancelling 


¢@ and transposing, (6) becomes 





which is the desired equilibrium condition for 
the rectangle. 

We shall next obtain the equilibrium condi- 
tion for a triangular area at the kernel. Let the 
triangular area Fig. 3 be broken up into ele- 
mentary bands of constant width dx by a suit- 
able family of curves, x being measured along 
the orthogonal trajectory which passes through 
the center of gravity. We shall suppose that 
is essentially constant over the triangle, its value 
Ho being taken as that at the center of gravity; 
also we shall suppose that the length of the ele- 
mentary band at x is given closely by 


L= Kx", (8) 


where K and m are constants. The total mag- 
netic flux in the triangle from the vertex to the 
elementary band is 


Ay Kx 


f Hy) Kx"dx = 
0 m+1 


hence the electric intensity at the band is 


1 Hy Kx" fw 
wat » 
Kx m+1 c 


in which the last two factors give the voltage 
across the band. Integrating (9) we obtain o,¢, 





ee 
<4) Fic. 3. Tri- 
se = angular area at 
= kernel. 
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Li; i 2 i—l n—l n—l 
. - (=) {3(toi+E os ) Aart (ae. +E os )Aart > (104% os) if (7) 
Wi; c k=l k=l l=i+1 k=1 





the total electric flux i in the triangle, thus 








W nj wx HwW,,; 
ono = Fan, = §G) 
0 ae 2c(m+1) 
Noting the last term of (4) we replace Ho by 
its value 
ow n—l 
= (aon >» a), (11) 
c keel 
and obtain, after cancelling ¢ and transposing, 
2(m+1)on 
W,7= a ’ (12) 


@ 2 n—1 
Cc kel 


Dividing the left hand side by L,;W,;, and the 
right by its equal KW,,;"*', and noting that 


Anj=KW,j""/(m+1), 


we finally obtain 


Lain 1 w 2 n—1 
We; =3() (wtEe)4n 9 


which is the desired equilibrium condition for 
the triangle. We note that this expression con- 
tains neither K nor m; also, in comparing (13) 
with (7) we note the presence of the factor $ in 
(13). 

As usual with flux plotting methods, the actual 
process of making the plot is one of systematic 
trial and error. The choice of the o,’s is a matter 
of convenience. If o; is put equal to one, the ith 
tube is a unit tube; however, at the edge of the 
plot and in the weaker parts of the field it is 
desirable to use fractional tubes so that more 
lines show on the plot. 
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In constructing the plot a step-by-step process 
can be used whereby the values of L,;, Wi;, and 
o; taken from one approximate plot are placed 
in (7) and (13) to obtain improved values of 
Li3/W.; and o;, which are then made the basis 
of the plot which constitutes the next approxi- 
mation. The first plot would be largely assumed 
from physical considerations. One step of the 
step-by-step process will now be described. 

First we note that the o;,’s associated with a 
plot may be multiplied by any constant without 
altering the plot. This is evident from (7) and 
(13), and also from the fact that the effect is 
merely to magnify the field strength. Values of 
Li;, Wi;, and o; taken from the given approxi- 
mate plot are now placed in the right hand sides 
of (7) and (13) to give for each rectangle or tri- 
angle the value of L;;0;/W:; to be used in con- 
structing the next approximate plot.‘ In so doing 
the unknown factor (w/c)? is replaced by an 
arbitrary constant multiplier k; hence the values 
obtained are of the form k8,;; where the §;;’s are 
computed numbers. Dividing by the measured 
values of L;;/W,; we obtain a set of values of 
kBi;W.;/L.;, which we denote by 4;;. If k can be 
chosen so that ¢;;=o; for all rectangles and tri- 
angles, then this value of k can be taken as 
(w/c)*, all conditions are satisfied, and the given 
plot is final. Actually this cannot be done except 
with the final plot ; hence in general we choose the 
new set of o;’s so that o; approximates all the 
é;;'s in the ith tube. Since o; is of the form ky; 
this is equivalent to choosing 7; so as to approxi- 
mate all the values of 8;;W;;/Li; in the ith tube. 
This can be done by choosing an average value 
or, perhaps better, a weighted average with the 
areas as weighting factors. Dividing the com- 
puted values of L;;0;/W:;, namely RB;;, by these 
new values of o; we obtain the values §;;/7i;, 
which we use as a new set of L;;/Wi; values in 
constructing the new plot.‘ We note that these 
are independent of k, which may be given any 
convenient numerical value. 

The step-by-step process ends when a value of 
k can be chosen for which ¢;;= 0; for all rectangles 


‘Each computed value of Lijo;/Wi; is associated with 
the center of its rectangle or the center of gravity of its 
triangle, respectively. If the center of a new rectangle 
differs appreciably from that of the old, interpolation in 
this field of values may be necessary. The same is true of 
triangles. 
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and triangles, this value being (w/c)?, as stated 
above.’ The most accurate value of w/c is prob- 
ably given by equating the total electric energy 
to the total magnetic energy obtained from the 
final plot, as described in Section 4.* 


3. ALTERNATIVE PROCEDURE USING 
A PLANIMETER 


Let 
= distance from the kernel measured along a line of no 
work at the center of the jth strip, (14) 


#=0¢= total electric flux outside of the line of force at 
s;- Here ¢ is again the arbitrarily chosen ‘normal’ 
flux per tube. 

A;=area of jth strip from the kernel to s;; 


then proceeding as with (3), noting Ampere’s 
law, the magnetic flux density for any s; is 


$w/¢=cdbw/c. (15) 


The total magnetic flux in A; is therefore 


Aj 
= f siti 
c 0 


it being permissible to consider o a function of 
A; since both o and A; are functions of s;. Con- 
tinuing as with (6), noting Faraday’s law, we 
multiply by w/c to obtain the voltage 


o(°) J ii (16) 


across the strip at s;. Choosing for s; the value at 
the center of the ith tube, dividing by L,; to 
obtain the electric intensity, and multiplying by 
W; we obtain the flux oi¢ of the tube, thus 


Ww 2 Wi; Ai 
v= 0(=) 2 feds, 
Cc Li; 0 


where A,; is the strip area to the center of the 
ith tube. Cancelling ¢ and transposing we ob- 





5 From the derivation of (7) and (13) or, later, of (17) 
and (18) it is evident that the above method of successive 
substitutions is equivalent to using the first approximate 
electric field to determine the corresponding magnetic 
field, corresponding _ to which, in turn, is obtained an 
electric field, which is taken as the second approximation 
—and so on. A method of successive substitutions may 
not be desirable in some cases, in which the new set of 

Li;/Wé values may be chosen arbitrarily after comparing 
the new values computed as above with those of the 
preceding plots. 
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Fic. 4. Electric field in hole and slot type magnetron. 
Approximation 1. 


Lijoi w\? Aaj 
= (=) f odA jy 
Wi; c 0 


which is the desired equilibrium condition for a 
rectangle. 

We shall now derive a condition analogous to 
(17) which can be applied to the triangles at the 
kernel. Denoting by &; the voltage, given by 
(16), across the strip at the base of the jth tri- 
angle, and noting from (9) that for any approxi- 
mating function Kx” for triangle width the field 
intensity varies linearly with x, we have for the 
total electric flux in the triangle 

&; Wai x 6;Waj 


x . 


21s 


tain finally 





(17) 








or, substituting from (16), cancelling ¢, and 
transposing, 


Lajon 1fw\? (4* 
=— (*) f od A jy 
Waj 2\c 0 


where the upper limit of the integral is the area 





(18) 


of the triangle. This is the desired equilibrium . 


condition for the triangle. In contrast to (17) we 
note the factor } in (18), also the fact that L,; 
and W,, are the base and altitude of the triangle, 
respectively. 

In applying (17) and (18) a convenient set of 
o;'s for the various tubes is chosen, and A; is 
* obtained with a planimeter for values of s; cor- 
responding to lines of force on the* plot. These 
values of A; are then plotted against the corre- 
sponding known values of oc. This curve gives ¢ 
as a function of A;; hence a planimeter or ap- 
proximate integration can be used to give the 
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integrals in (17) and (18). Such a curve is drawn 
for each strip. 

The step-by-step process described at the end 
of Section 2 evidently applies here, (17) and (18) 
being used instead of (7) and (13). As before 
w/c is not determined until the final plot has 
been made, after which w/c is obtained by“equat- 
ing the total electric and magnetic energies as 
described in Section 4. 


4. ELECTRIC AND MAGNETIC FIELD 
ENERGY. DETERMINATION OF «w/c 


In any strip the electric flux passing through 
an elementary band bounded by two lines of 
force is ¢de. Dividing by ds; we obtain the elec- 
tric intensity. The total electric energy in the 
strip is therefore® 


Electric energy in jth strip 
¢? Aj da 2 
=— (=) dA; (19) 
8x 0 ds; 
Noting (15) we also have 


Magnetic energy in jth strip 


2 2 Aj 
-<(*) | o°d A ; 
8x \c 0 


Using the method of Section 3 a curve of ¢ asa 
function of s;can be plotted, and da/ds; obtained 
by approximate differentiation ; also a curve of ¢ 
as a function of A; is already available.’ The 
integrals in (19) and (20) can hence be obtained 


(20) 


6 The expression for electric energy obtained by means 
of the magnetic field using (16) merely duplicates (20), as 
is shown by the following. 


Electric energy in jth strip 


1(°T fo\ (45 oda 
=5. x [o(2) J. od A, | 
_  (o\? (Ai (45 de 
= -£(2) I, [U. odAs| Fr dAy, 


Integrating by parts noting that ¢=0 when A; takes its 


greatest value, this becomes 
Electric energy in jth strip 


2 2 Aj A Aj 
= - “) {[eJS, odA;|,'— J, ‘otda;\ 


2 2 Aj 
EC) feu, 
which is (20). 


7 Tables for approximate differentiation are given in 
Crout, “An application of polynomial approximation to 
the solution of integral equations arising in physical 
problems,” J. Math. and Phys. 19, 34 (1940). 
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using a planimeter or approximate integration. 
Adding the energies of the various strips the 
total electric energy and the total magnetic 
energy can be obtained, both being proportional 
to ¢. Equating these, ¢? cancels out, and we have 


i Aj f/da\2 4 
£ {ea 
0 ds; 
“. i pj 
> { o°dA; 
0 


the upper limits on the integrals being the strip 
areas, and the summations being over the strips. 
This method of determining w/c reminds one of 
Raleigh’s method in elastic vibration problems, 
and offers probably the most accurate way of 
using the flux plot to obtain the resonance 
frequency. 

We shall finally obtain expressions for the 
electric and magnetic energies which are ap- 
plicable if the individual rectangles and triangles 
are considered separately, as in Section 2. The 
electric energy in a rectangle is evidently 





(21) 


Electric —_ in rectangle 


(CA) (Cs) a, 


also, noting (3), the magnetic energy of a rec- 
tangle is seen to be 





Magnetic energy in rectangle 
¢? w 2 —1 2 
-© (2) (d+ on) LaWn (23) 
8c k=l 


From (9) it follows that in a triangle the elec- 
tric intensity is of the form K.x where K; is a 








APPROXIMATION 2 


Fic. 5. Electric field in hole and slot type magnetron. 
Approximation 2. 
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constant. K, can be determined by the relation 


Wai K, W,7 
On = Kixdx = “?<s 


0 
whence the electric intensity becomes 


2onbXx 
E.;=—— 


2 
nj 


The electric energy in the triangle is therefore, 
noting (8), 
Electric energy in triangle 

1 Wail 2onbXx 
” Se W,7 
= (onp)*K Wj” '/24(m+3), 





2 
Kx"dx 


; ae (24) 
Electric energy in triangle 


= (ond)? Lnj/24(m+3) Waj, 
where L,; and W,; are the base and altitude of 
the triangle, respectively, and K and m are the 
values applicable to the jth triangle. 
In the triangle the magnetic field is closely 


equal to the constant value (11) ; hence the mag- 
netic energy is closely 


Magnetic energy in triangle 


2 2 n—1 2 Wai 
8r ic bel 0 


Magnetic energy in triangle 


2 2 n—1 2 LnjWaj 
mses) (eZ) Gey 
TX 





k=l (m+1) 








APPROXIMATION 3 


Fic. 6. Electric field in hole and slot type magnetron. 
Approximation 3. 
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Fic. 7. Electric field in vane type magnetron. 


5. APPLICATION TO A SPECIFIC PROBLEM. 
HOLE AND SLOT TYPE MAGNETRON 


The procedure used to obtain the field in the 
cavity shown in Fig. 2 will now be outlined.* 
The approximation that was used for the electric 
field in order to get the step-by-step process 
started was a uniform electric field in the gap 
between two adjacent anode segments, the elec- 
tric field being taken as zero everywhere else, 
Fig. 4. Data from this first approximation are 
now inserted in the right-hand sides of (17) and 
-(18) as described in Section 3 to obtain corrected 
data for the plot which constitutes the second 
approximation. In so doing we note that o is 
zero in the anode-cathode region, and increases 
. linearly through the gap to a value which re- 
mains constant in the well. It follows that 
StiodA; and hence the computed values of 
Lijo;/Wi; are constant in the anode-cathode re- 

*The work outlined in this section including that 


involved in cinatrneeing, the flux plots shown in Figs. 5 
and 6 was done by Mr. F. S. Holt. 
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gion, and that in the well Lijo;/Wi; and Lajon/ 
W,j are proportional, respectively, to A,; and 
A,;/2. The flux plot which constitutes the second 
approximation was begun at the inner (cathode) 
side of the gap, the work being carried both in- 
ward and outward from this point. We note that 
in the anode-cathode region this plot can be a 
lattice of curvilinear squares, which gives a field 
satisfying Laplace’s equation since L,;o;/W4; is 
constant in this region; also we note that the 
choice of the o,’s is a matter of convenience. The 
second approximation is shown in Fig. 5.° 

Continuing the step-by-step process, data from 
this second approximation are now inserted in 
(17) and (18) to obtain data for a third approxi- 
mation. In so doing each new value of o; is taken 
as the average of the ¢4;;'s of the elementary 
areas in the 7th tube. After the third approxima- 
tion, Fig. 6, was obtained, it was found that the 
ratio of o; to W;;8;;/L:; was constant for all the 
elementary areas to within the accuracy ex- 
pected of the graphical process. This third ap- 
proximation, Fig. 6, was therefore taken as the 
final result, the constant ratio being (w/c)*. The 
frequency corresponding to the value so ob- 
tained for w/c is 2.818-10° c.p.s.; that obtained 
by equating electric and magnetic energies using 
(21) was 2.825-10° c.p.s. The various tube di- 
mensions were 


Cathode radius 0.1182” 
Anode radius 0.3175” 
Gap width 0.066” 
Radius to center of well 0.640” 
Well diameter 0.398”, 


The values of the o;’s for the final plot, Fig. 6, 
are shown in Table I. 


6. APPLICATION TO A SPECIFIC PROBLEM. 
VANE TYPE MAGNETRON 


If the shape of the electric field is known ap- 
proximately, the work required to obtain a plot 
can be greatly reduced.” For example in the 


® The labels shown in Figs. 5 and 6 are those that were 
used in making the plot. The strips are indicated by 
primes, the tubes are numbered, and subdivisions are 
indicated by letters or, if obvious, are not indicated. This 
notation differs from that used in the derivation of (7), 
(13), (17), and (18); but is convenient for actual con- 
struction. 

10 The work outlined in this section including that 
involved in constructing the flux plot Fig. 7 was done by 
Mr. F. E. Bothwell. 
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TABLE I. Values of oy. 


TABLE II. 











Tube ov Tube ov Tube Gi 

Triangle .0625 4 .999 11 .999 

Ic 125 5 .999 12 1.00 

1b .250 6 .999 13 1.00 

la .250 7 .985 14 1.00 

2b 499 8 .989 15 .836 

2a 499 9 .994 

3 .997 10 .997 








vane type magnetron, Fig. 7, the electric lines 
of force in most of the region between the vanes 
evidently approximate arcs of circles centered 
at the cathode; also in the anode cathode region, 
where the magnetic field is weak, the electric 
field approximates a solution of Laplace’s equa- 
tion. With these facts in mind Fig. 7 was drawn 
for the case where the cathode and all vanes but 
one are at the same voltage, no attempt being 
made to satisfy (7). Omitting the line of no work 
let us temporarily consider the field as occupying 
a single strip; then Eqs. (7) form a set of twelve 
homogeneous, linear, algebraic equations for the 
o;'s. These can be conveniently solved by suc- 
cessive substitutions by first placing all o,’s on 
the right-hand side of (7) equal to unity. The 
a; s so obtained are inserted in the right-hand 
sides of (7) to give the next approximation, and 
so on. In each step the k which replaces (w/c)? 
was given that value which would make the 
average of the computed a,’s equal to unity. 
This method of successive substitutions has the 
physical significance described in footnote 5; its 
convergence was rapid, and the results obtained 
are shown in Table II. When put equal to 
(w/c)? the last obtained value of k corresponds 
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Values of o¢ 





step 
4 1 2 3 4 5 
1 1.00000 .97989 .97990 .97998 .97999 
2 1.00000 -97444 .97484 .97489 .97491 
3. =1.00000 .96696 .96786 .96792 .96793 
4 1.00000 .96162 -96293 .96302 .96300 
5 1.00000 1.14617 1.14824 1.14833 1.14835 
6 1.00000 1.00887 1.01120 1.01130 1.01127 
7 1.00000 1.06811 1.07140 1.07148 1.07150 
8 1.00000 1.24744 1.25118 1.25127 1.25128 ° 
9 1.00000 1.32741 1.33002 1.33001 1.33002 
10 1.00000 1.14551 1.14404 1.14388 1.14386 
11 1.00000 .83281 .82418 .82390 -82388 
12 1.00000 .34074 .33420 33401 .33400 
Table of values of k 
step 1 2 3 4 5 
k 6.0374 5.7562 5.7517 5.7516 








to a wave-length of 2.6200 cm. The (probably) 
more exact value obtained by equating total 
electric and magnetic energies as described in 
Section 4 is 2.6456 cm. 

Let us finally try to improve the accuracy of 
the plot by subdividing and proceeding as in 
Section 5. We therefore draw the line of no work 
that was omitted above; and, placing the o;’s of 
Step 5 in the right-hand side of (7), determine a 
new set of values L;;/W;;. These new values, 
however, duplicate the old values in all tubes 
except the first four, in which there is a dis- 
crepancy of one percent. Since Li; and W,; 
cannot be measured to within one percent, no 
alteration of the plot is indicated, and it is there- 
fore taken as final. It is probable that the pro- 
cedure used in this section on the vane type mag- 
netron can be used to advantage on the hole and 
slot type magnetron considered in Section 5. 
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A Magnetron Oscillator for Dielectric Heating 


R. B. NELSON 
Research Laboratory, General Electric Company, Schenectady, New York 


(Received January 31, 1947) 


Magnetron oscillators can provide very high frequency power for dielectric heating. Factors 
determining the choice of frequency for particular jobs are discussed. The design of a magnetron 
having 5 kw continuous output at 1050 megacycles is described, and performance data on the 


tube are given. 


INTRODUCTION 


HE method of heating non-conducting ma-: 


terials by dielectric loss,in a high frequency 
electric field is becoming an important industrial 
process. Many of the materials heated by this 
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GEOMETRY OF CAVITY AND LOAD 








E2 
E= MAX. ' 
1 ! 
re) : i iL 
AMIS , WALL 
LOAD 


RADIAL VARIATION OF FIELD 


Fic. 1. Distribution of heating in cylindrical load. In a 
simple cylindrical cavity where the electric field is produced 
by radial currents in the walls, the r.m.s. field isa maximum 
at the center and decreases radially as a Bessel function. 
The rate of decrease depends upon the wave-length in the 
dielectric, giving the discontinuity in slope of the curve of 
EF plotted against radial distance. Rate of production of 
heat per unit volume is proportional to E*. 
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method, such as plastics, are of such low loss that 
it is hard to heat them as rapidly as desired. The 
rate of heating increases with the electric field 
strength and, in practically all materials, with the 
frequency. Since the usable field is definitely 
limited by arcing between the electrodes, the only 
remaining way to speed up the heating process is 
to increase the frequency. 

Commercial dielectric heaters have been re- 
stricted to frequencies below about 40 megacycles 
by the limitations of available vacuum tubes. 
However, the recent development of magnetron 
oscillators has made available considerable power 
at frequencies up to some 10,000 megacycles. 
Using magnetrons for heating thus allows the 
choice of almost any desired frequency. 


CHOICE OF FREQUENCY 


It must not be assumed that the highest fre- 
quencies are always best. While heating rates go 
up with frequency, a point is reached where non- 
uniformity of heating limits the usable frequency. 
Two effects contribute to non-uniformity. First, 
in the case of high loss materials such as most 
foods, there is the decreased heating toward the 
center of large objects because the power is 
absorbed by the surface layers before it can 
penetrate. Second, in the case of low loss ma- 
terials such as plastics, the wave-length of the 
heating power limits the volume of material that 
may be heated uniformly. To apply power 
efficiently to low loss material, it is necessary to 
surround the material with reflecting walls so 
that the electromagnetic wave passing through it 
is not radiated and lost after one passage. The 
reflecting walls set up a standing wave in and 
around the load material. This wave has in 
general a three-dimensional space variation of 
electric field, with periodic maxima and minima. 
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Since the dielectric heating at any point is pro- 
portional to the square of the electric field 
strength, widely varying temperatures will be 
produced in the material if its dimensions are 
comparable with a wave-length. As an example 
of how the standing-wave pattern limits the 
usable frequency, consider a cylinder of load 
material located axially in a cylindrical resonant 
cavity, as shown in Fig. 1. Exciting the cavity in 
its lowest mode with axial electric field, the field 
is a maximum at the center and falls off radially 
as 
E = EmaxJ0(6.28r/X), 


where E is the electric field strength at radius r 
from the axis, Emax is the field on the axis, and A 
is the wave-length in the material. If it is required 
to heat the edges of the load 90 percent as much 
as the center, 


E?/ Emax”? = J?(6.287r/r) =.90 at the edges, 


and solving the above equation for a load of 
radius a 
a=.0725iX. 


If, for example, the load is a piece of phenolic 
plastic with dielectric constant 4.0 and radius 
1.0 cm, 

\=)o/4.0}, 


where Xo is the free space wave-length corre- 
sponding to the frequency used. Solving for Ao 


\o= 27.6 cm. 


Thus, the maximum frequency allowed in this 
case would be 1087 megacycles. 


OSCILLATOR DESIGN CONSIDERATIONS 


On the basis of considerations such as the above 
example, it appeared that many dielectric heating 
processes could be done well at a frequency of 
about 1000 mc. Therefore, a magnetron oscillator 
was developed in the General Electric Research 
Laboratory to give 5 kw output at 1050 mc. In 
designing this oscillator, consideration had to be 
given to several special requirements of heating 
apparatus which differ from other uses of 
magnetrons. 

The most important specific requirement is 
that the heating oscillator must be capable of 
operating into any load impedance, including the 
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Fic. 2. Photograph of the 5-kw magnetron, showing the 
external appearance and also an axial section. 


completely unloaded condition. This is caused by 
the non-uniformity of successive charges of load 
material, to the certainty of bad adjustment of 
matching transformers by unskilled operators, 
and the possibility of accidentally running the 
oscillator with no load at all. To meet these con- 
ditions, the anode should be capable of dissipating 
the entire input power to the tube. The output 
seal and transmission line to the load must stand 
very high voltages and currents associated with 
high standing wave ratios. The cathode should 
not be damaged by the excessive power it re- 
ceives from back-bombardment by electrons, 
which is a function of the load impedance as well 
as the plate voltage and current. 

Although the multiple-cavity magnetron has 
many resonant modes, efficient oscillation occurs 
only in the ‘‘z’’ mode where each anode is 180° 
out of phase with those adjoining it. When this 
mode is overloaded, as may occur by poor load 
matching, the oscillator will tend to jump to 
another mode, which is usually completely 
unloaded, and may continue in the spurious 
mode, even if the mismatch is corrected, until it 
is turned off and restarted. The accepted cure for 
mode jumping is to strap the anodes together in 
two alternate sets with low impedance con- 
ductors. This raises the frequency of all modes 
except the “‘x’’ mode to a point where the 
oscillator does not readily excite them. The 
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Fic. 3. Construction of water cooled magnetron. Sec- 
tional, views showing the tubular anodes and water circu- 
lating channels. The O.D. of the copper shell is 34’’. 


heating magnetron, as distinguished from the 
radar transmitter, runs continuously, so once it 
has jumped modes it does not have a chance to 
restart properly. Also the wide range of load 
impedances favors mode jumps. These factors 
necessitate very tight strapping. It is desirable 
that the next higher frequency mode be at least 
twice the frequency of the operating mode. In 
_ the oscillator to be degcribed, this separation was 
attained by a new method of strapping. 


A 5-KW, 1050-MC MAGNETRON 


The magnetron is shown in the photograph, 
Fig. 2, and constructional details in Fig. 3. It has 
10 anodes, and so may be run at relatively low 
voltage and magnetic field. Rated operating con- 
ditions for the full 5 kw output are: plate voltage 
5.0 kv, plate current 1.5 amp., magnetic field 
1500 gauss. Water cooling was adopted for the 
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anodes, and this allowed making the anodes of 
U-shaped pieces of copper tubing through which 
the water flows. Channels in the copper tube shell 
form headers for water distribution. 

In addition to providing excellent cooling, the 
tubular anodes have a minimum of inter-anode 
capacity, thus allowing a large fraction of the 
capacitance of the resonant circuits to be pro- 
vided by low impedance straps. The resulting 
tight strapping gives a frequency ratio of 1:2 
between the operating mode and the next higher 
mode. The double ring straps are looped inside 
the anode U’s. This allows the straps to be 
attached near the centers of the anode faces, 
eliminating unbalanced r-f voltages on the 
anodes. Also, the anodes shield the cathode and 
the rest of the tube from the fundamentally 
unbalanced double ring straps. The ten resonant 
circuits formed by the anode loops extending 
radially inward from the tube shell are electrically 
similar to conventional vane resonators. 

The over-all size of the magnetron is deter- 
mined by the output seal, which must be large to 
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Fic. 4. Magnetron mounting and connections. The tube 
is slid downward into the cylindrical magnet, simul- 
taneously plugging in to the coaxial output transmission 
line which is rigidly mounted on the magnet. The iron 
magnetic circuit is shown in heavy cross-hatching. 
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handle the power. In this cylindrical structure the 
seal occupies the entire diameter of the tube shell. 
The coaxial arrangement of the seal with the 
cathode and magnetic pole pieces allows the tube 
to be plugged into its magnet and output 
transmission line simultaneously. Figure 4 shows 
the method of mounting the tube. An electro- 
magnet is used, consisting of a solenoidal coil 
with an iron shell whose ends line up with the 
iron pole pieces inside the tube. Air is blown 
through the center conductor of the output line 
to cool the seal, the inverted cup shape of the seal 
giving a good distribution of the cooling air as 
well as protecting the seal from breakage. 

The power output is obtained by attaching a 
flat strap from the center conductor of the seal to 
one of the anodes. The exact point of attachment 
is selected to provide proper loading of the 
oscillator when the 52.5-ohm transmission line on 
the output has no standing wave. 

The cathode is a simple spiral of pure tungsten 
wire with molybdenum end hats. It requires 600 
watts heating power. In a magnetron having 
considerable electron bombardment of the cath- 
ode, it is desirable to use a very inefficient high 
wattage cathode, so that the back heating be- 
comes a small fraction of the total heating power. 
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Fic. 5. Performance chart. Solid lines are contours of 
constant magnetic field, dashed lines are constant power, 
and dotted lines constant efficiency. Power, in watts, and 
efficiency measured are indicated on the experimental 
points. This chart is taken with a load matched to the 
transmission line. 
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Fic. 6. Rieke diagram. The effect of load impedance on 
the efficiency and wave-length of the oscillator is best 
shown by this chart. The data were taken at fixed condi- 
tions of 1.0 amp. plate current and 1200 gauss magnetic 
field. Upper figures are wave-length (add 29.000 cm) and 
lower figures percentage efficiency. Constant efficiency and 
constant wave-length contours are shown by solid and 
dotted lines, respectively. The radial coordinate is voltage 
reflection coefficient, with the scale marked in standing 
wave ratio in decibels. The azimuthal coordinate is the 
phase angle of the reflection coefficient, with the scale 
marked in wave-lengths displacement of the standing 
wave minimum from an arbitrary zero. 


MAGNETRON OPERATING CHARACTERISTICS 


The performance of the magnetron as an 
oscillator is graphically presented by several 
charts. Shown in Fig. 5 is the ‘performance 
chart” in which power output (upper numbers) 
and efficiency (lower numbers) are plotted 
against plate voltage and current. Lines of con- 
stant magnetic field are also shown. 

Another important oscillator characteristic is 
shown in Fig. 6, the ‘Rieke diagram,” in which 
efficiency and wave-length are plotted as a 
function of load impedance for constant current 
and magnetic field. The radial coordinate is 
standing wave ratio in the output line (in voltage) 
plotted on a scale proportional to reflection 
coefficient. The angular coordinate is the phase 
angle of the standing wave minimum. Thus the 
position of a point on the chart represents a 
complex load impedance, which might be trans- 
formed into resistance and reactance components 
by superposition of a Smith chart. 
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Fic. 7. Effect of plate voltage on output power. With a 
constant magnetic field as obtained from a permanent 
magnet, power output varies rapidly with plate voltage, 
following the dotted curve. An electromagnet excited by 
the magnetron plate current reduces this variation. 


It is apparent from Fig. 5 that the plate current 
changes very rapidly with plate voltage at 
constant magnetic field. That is, the variational 
resistance for plate modulation is around 250 
ohms. This is a good situation if one desires to use 
amplitude modulation, as the modulating power 
is low. However, in dielectric heating, a constant 
r-f output is desired, and account must be taken 
of variations in line voltage. It is obvious from 
the curves that a small variation in line voltage 
and hence in plate voltage will produce a large 
change in plate current and in output power. 
Figure 7 shows this relation, neglecting the 
regulation of the plate supply rectifier. 

A method of greatly reducing the effect of line 
voltage fluctuations and, incidentally, plate 
supply ripple, is to run the tube’s electromagnet 
in series with the low voltage end of the plate 
supply. This gives a magnetic field directly pro- 
portional to the plate current. Since the plate 
voltage is almost directly proportional to the 
“magnetic field, the plate current will in turn be 
proportional to plate voltage. The resultant 
ohmic plate resistance characteristic produces the 
variation of r-f power with plate supply voltage 
shown by the solid curve in Fig. 7. 

The relationship between d.c. plate current, 
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plate voltage, and magnetic field also affects the 
way in which the magnetron output depends on 
the load impedance. The data plotted in Fig. 6 
were taken under the easily reproducible condi- 
tions of constant current and magnetic field, 
letting the voltage find its own value. The 
voltage varies much less at constant current than 
the current would at constant voltage. 

In practical operation we do not have either a 
constant current or constant voltage power 
supply, and the regulation characteristics of the 
supply affect the relationship of magnetron 
power output and load impedance. It is obvious 
that this is a very important relationship for 
heating applications. 

With the magnetic field supplied by a series 
electromagnet, some of the effects of power 
supply regulation are simulated, and r-f power 
output tends to be stabilized. For a constant 
voltage power supply and series magnet, the 
variation of output power with load for a pure 
resistive load is plotted in Fig. 8. For comparison, 
there is also shown the relations obtained with 
constant voltage power supply, and with constant 
current, 
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Fic. 8. Effect on output power of varying pure resistance 
loads. For loads which present a pure resistive impedance to 
the magnetron, the power delivered is plotted as a function 
of the ratio of load resistance R to the characteristic 
impedance Zp» of the transmission line. At constant mag- 
netic field, the power increases with load resistance if the 
magnetron plate voltage is held constant, but decreases if 
the current is constant. Using a series electromagnet with 
constant voltage supply, the better characteristic shown 
by the solid line is obtained. 
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It has been found that if plate voltage is applied 
to the magnetron with no magnetic field but 
with the magnet in series, the tube will start in a 
very high frequency mode of oscillation, drawing 
plate current so that the voltage never rises 
enough to start the desired mode of circuit 
oscillation. To prevent this, it is necessary to 
have a field of at least 600 gauss present before 
plate voltage is applied. A rectifier across the 
magnet as shown in Fig. 9 supplies this starting 
current, and normally stops supplying current 
when the magnetron current builds up a voltage 
drop across the magnet coil greater than the 
output voltage of the rectifier. 

The magnet supply rectifier can be made to do 
another useful job as a fine control for plate 
current and output power. In the plate supply 
rectifier, continuous voltage control is expensive, 
but a few taps for step control are cheap. Now at 
any given plate voltage, the magnetron is ap- 
proximately a constant-magnetic-field device, 
independent of plate current. If the plate current 
supplies part of the magnetic field and the rest 
comes from the magnet rectifier, any increase in 
rectifier current will tend to increase the magnetic 
field and thus will cause a corresponding decrease 
in magnetron current. In this way continuous 
adjustment of output power is obtained by 
regulating the 100-watt magnet rectifier. Figure 9 
shows the Variac connection used for control. 


R-F LOADING EQUIPMENT 


The problem of carrying the 1050 mc power 
from the magnetron and coupling it into a load is 
one which has almost as many answers as there 
are varieties of material to be heated. The power 
is taken from the oscillator on a stub-supported 
coaxial line, the start of which is shown in Fig. 4. 
A line size of 3§’ O.D. with 13” center conductor 
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Fic. 9. Current supply for electromagnet. An auxiliary 
d.c. supply for the magnet is necessary to assure proper 
starting. This variable-voltage rectifier may be used to 
control the magnetron power output. 


has been adopted as adequate for this power 
level and frequency. For most loads, a large 
step-up in impedance is necessary to match the 
line to the load. This usually takes the form of a 
resonant cavity tuned to the operating frequency 
and coupled to the transmission line by a loop or 
antenna. 
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An Approximate Solution of the Differential Equation of the Ultracentrifuge 
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In determining molecular weights with the ultracentrifuge by the equilibrium method, it is 
necessary to drive the rotor at relatively high speeds for very long periods of time. It would be 
a real advantage to be able to deduce molecular weights from the distribution of concentration 
within a cell placed in a rotor without waiting for equilibrium. The solution of a certain dif- 
ferential equation gives the theoretical distribution within the cell at any time, and a com- 
parison of the experimentally determined concentrations with those deduced from theory 
should yield molecular weights without waiting for equilibrium. The exact solution of this 
differential equation is known, but it is not suitable for numerical work. The present paper 
presents an approximate solution from which numerical results can be obtained with ease. 





(y= of the most important methods of meas- 
uring high molecular weights is by the use 
of the ultracentrifuge as developed by Svedberg 
and Beams. A dilute solution of the substance 
whose molecular weight is to be determined is 
placed in a small cell, with radial sides and cylin- 
drical ends, near the periphery of the rotor, and 
under the influence of the centrifugal force the 
distribution of concentration changes from an 
initial uniform distribution to an exponential one. 
When this final equilibrium state has been 
reached the concentration “‘c,”” at any point in 
the cell, is given by the expression 


M (i — Vp)w?(r? — 10) 
=C9 exp| | 
2RT 





(1) 


where “ M”’ is the molecular weight of the dis- 
solved substance, ‘V”’ its partial specific vol- 
ume, p the density of the solvent, w the angular 
velocity of rotation, and ‘‘r’’ the distance of any 
point in the cell from the axis of rotation. “* C,”’ 
is the concentration at the point ‘79.’ “‘ R” is the 
gas constant per mole and “7” the absolute 
temperature. The equilibrium state is reached 
when the tendency for solute particles to settle 
to the outer part of the cell because of the cen- 
trifugal force is exactly equal to the tendency for 
_them to diffuse inwards because of the concen- 
tration gradient. It is obvious that Eq. (1) gives 
““M”’ in terms of quantities that can be meas- 
ured, although the experiment is one that can 
only be performed successfully with faultless 
optical and mechanical equipment. 

The time required for the attainment of the 
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equilibrium state is often very long—sometimes 
the centrifuge must run for many days, or even 
weeks, at constant speed and at a constant tem- 
perature. If the time of a run could be shortened 
this method of molecular weight determination 
would be much more useful. The experimental 
difficulties would be greatly reduced and the 
number of runs made with any centrifuge would 
be increased. This latter consideration is an im- 
portant one in view of the many substances that 
could be studied by this method. If the distribu- 
tion of concentration within the cell at any time 
after the start of an experiment could be deduced 
theoretically, a comparison of the measured and 
computed distribution might enable ‘ M”’ to be 
found without waiting for equilibrium. The 
function c(r, t), which gives the concentration at 
any point at any time, would be required for 
this purpose and this is obtained by solving the 
following differential equation : 


10 Oc Oc 
-— (d= —wrsc)r| =—. (2) 
ror or ot 


In (2), ‘D” is the diffusion constant of the solute 
and ‘“‘s’”’ the sedimentation constant, i.e., the 
velocity of settling of the dissolved substance in 
unit gravitational field. The quantities “‘D’’ and 
‘*s”” are related to each other. It is possible to 
show that 


D/s=RT/M(1i— Vp). (3) 


A proof of this relation is given in a previous 
paper by the author.! The differential Eq. (2) was 


!W. J. Archibald, Phys. Rev. 53, 746 (1938). 
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first derived by Lamm.’ A simple derivation is 
to be found in reference 1. 

A treatment of Eq. (2) which is free of any 
mathematical approximations has been given and 
is to be found in two places in the literature.* 4 
While the solution given is exact and can be 
employed for purposes of computation, yet the 
cost in time of getting numerical results is so 
great that little use has been found for it. The 
expression obtained for c(r, ¢) is simple in form 
and would be most serviceable if the mathe- 
matical functions which it contains could be 
looked up in tables. However such tables are 
not available and the task of compiling them 
would be a major project in applied mathe- 
matics. It is the purpose of this paper to derive 
another expression for c(r,¢) which, while not 
exact, is felt to be good enough for many pur- 
poses, and which contains functions that are 
tabulated. 

In dealing with Eq. (2) it is more convenient 
to use ‘“‘z”’ and “7” rather than “r” and ‘‘?” as 
the independent variables, where 


ws M(1i— Vp)w? 
s=—r’= r? and r=2wst. (4) 
2D 2RT 





As 7;=r=re2 where r; and rz are the inner and 
outer limits of the cell, respectively, then 
a=z=b where : 








“ M(i- en 
b 2RT 


9 


ro~ 


In terms of these variables we find 


b 
f Mle 1, z)dz 


a 


c(z, 7) b—a 











e+> 
Co e—er 3 n=1 
f Ute, 1, z) }*dz 


XM(an, 1, ze", (5) 


where M(a, y, 2) is a solution of the differential 
equation 





d?M v dM «a 
(=-1 —-—-M=0 (6) 
dz? z dz z 





2 Lamm, Ark. Mat. Astron. Fysik, 21B, No. 2 (1929). 
3W. J. Archibald, Phys. Rev. 54, 371 (1938). 
*W. J. Archibald, Ann. N. Y. Acad. Sc. 43, 211 (1942). 
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and the a,’s are the “eigenvalues” of the pa- 
rameter a. The a,’s are those particular values 
of a which give the functions M(a,, 1,2) the 
property that 


d 
—M (an, 1, 2) = M(an, 1, 2) (7) 
dz 


for z=a and z=). “co” is the initial uniform 
concentration. (The reader is asked to consult 
reference 3 or 4 for a detailed proof of (5) and 
a full explanation of the results stated here.) 
M (an, 1, 2) is an abbreviation for 


F(an, 2, 2) 


F(a,, 1, 2) -——————- 
aa: <r 


W(an, 1,2), (8) 


where F(an,, y,2z) and W(a,, 7, 2) are the two 
independent solutions of (6) which are to be 
found written out in full.** It is the difficulty of 
obtaining numerical values for the F and W 
functions which limits the usefulness of Eq. (5). 

It will be our purpose to see if it is possible to 
derive some expression for the_functions M(a,, 
1,2), other than that given in (8), and which 
contains only well-known functions. This al- 
ternative expression will of necessity be an ap- 
proximation and a careful study of the magni- 
tude of, the different quantities involved will be re- 
quired to make the approximation a good one. 
A further purpose will be to discover a simple 
method for determining the eigenvalues ap. 

In (5) the first term, in e*, gives c/co for the 
equilibrium state. In an actual experiment_it is 
not desirable to have the final concentration 
vary too widely throughout the cell. If one were 
to accept as a reasonable set of conditions that 
finally the concentration at rz be three times 
that at r, then 


e/e~3, or (b—a)—=1.1. 


This serves to show the narrow limits within 
which the variable ‘‘z” is confined. Furthermore, 
since 
b/a=r;?/r2?, 
we have 
ry r,? 


a=(b—a) ™1.1 


re? —r;" re? —r;" 





The cells used are usually small and one might 
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choose as typical values the following: r2=5.0 
cm and r;=4.5 cm. This gives a~4.7. As a re- 
sult of these considerations we are entitled to 
treat the variable ‘‘z’’ as being moderately large 
(of the order of 5.0) and as being confined within 
a short range of values. The results of a computa- 
tion on a typical experiment** show that 
5.0=2=6.2 and these figures are probably typi- 
cal of most cases in which the equilibrium method 
is used. 

A further point of interest has to do with the 
quantities a,. The exact theory gives a=1 as 
one of the eigenvalues. This a, of course, gives 
rise to the first term in (5). All the other a’s have 
large negative values. In the computations re- 
ferred to above it was found that a,;= —38.7. 
a2 was not determined, but the indications were 
that it was less than — 100. 

The first term in (5) presents no difficulty: 
but the terms in the summation part of the ex- 
pression are bothersome. Thus we will confine 
our attention to these terms. This means that 
we must study Eq. (6) (with y =1), remembering 
that @ is very large and negative and that “2” 
is moderately large and of restricted variation. 
In (6) let 

M =e?!*z-4V. (9) 


The equation for “ V” then becomes 
d?V 1 1 (a—}) 


——- — V=0. 
4 42? Zz 





(10) 
dz? 


The considerations given above suggest that it 
is permissible to drop the term 1/42’. This term 
is very small in comparison with —(a—4)/z and 
is even quite small in comparison with 3. Then 
(10) becomes 


avy 1 (@-%) 
—+| omni | v=o. 
dz 4 2Z 





This is still an inconvenient equation to deal 
with but it can be further simplified by replacing 
the term in braces by 


“ee 





(11) 


“oft 


This term is a very close approximation to 
—1—(a—})/z when 25, a is large and b—a is 
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TABLE I, 

















z F(z) f(z) approx. 
5.0 7.59 7.56 
$2 7.29 7.27 
5.4 7.01 7.00 
5.6 6.75 6.75 
5.8 6.51 6.52 
6.0 6.28 6.30 
6.2 6.07 6.10 
(a+b) 1 
f(z) approx. = { oo =——(a —}) \ > 


With a = —38.7, a =5.0, and b =6.2. 


of the order of unity. As an example we may use 
the numerical data from reference 3 or 4. Putting 
a= —38.7,a=5, and b=6.2 with 5.0=:s=6.2 we 
obtain the results given in Table I. It will be 
observed that the maximum difference between 
the two quantities is less than 0.5 percent. 

Thus if one uses 


@V a? 


dz? 3 





(12) 


as the differential equation for the function ‘ V”’ 
it seems reasonable to suppose that (9) will give 
a close representation of the function ‘ M”’ for 
“*s”’ between “a” and ‘‘b.” 

Now Eq. (12) has for its most general solution 


V(z) =3!{ AJ,(202')+BN,(202')}, (13) 


where J; and JN, are Bessel functions of the first 
and second kind, respectively, and A and B are 
arbitrary constants to be determined by the con- 
ditions of the problem. Combining (9) and (13) 
we have 


M(a,z) =e7!?{AJ,(202!)+BN,(202')}. (14) 


It must be emphasized that this expression can 
only be made to approximate M(a, 1, z) (by the 
proper choice of A and B) for large values of a 
and for z between a and b. For any other range 
of variables it will undoubtedly be wide of the 
mark. Even subject to these restrictions there is 
no reason to suppose that e*/?J,(20z') is an ap- 
proximation to F(a,1,z) or e*/2N,(2ez') to 
W(a, 1, z). All we are entitled to expect is that a 
linear combination of F(a,1,z) and W(a, 1, z) 
can be replaced by some linear combination of 
the functions e*/2J,(20z4) and e*/2N,(20z?). 
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We will now adopt for the function M(ap, 1, 2), 
which occurs in Eq. (5), the form given in (14). 
However it is still necessary to discover what 
values of a, to use. (It must be remembered that 
@ is contained in the quantity o.) It will be re- 
called that the function M(a,z) must be such 
that 


adM/dz=M (15) 
for z=a and z=b. This equation will be found to 
hold only when the a’s are properly chosen. Let 

Z,(202!) = AJ, (202) + BN, (202?) 
and in order to be specific let N,(20z2') be the 


Neumann function whose properties are well 
known. Then 





M(q, 2) =e*!2Z;(20z?). (16) 
Substituting (16) in (15) one gets 
dZ, zi 
—-—Z, =(. (17) 


d(2ez') 2¢ 


(In the following the argument of Z, J, and N 
will be omitted.) But 





dZ, 1 
=Z)—-——Z. 
d(202') 2o23 
Therefore, (17) becomes 
Zo—Z, tan¢=0, (18) 
where 
1si 
tang =—(—+2!). (19) 
2a\ 23 


In terms of the J and N functions (18) becomes 
A(Jo—Ji tant)+B(No— JN, tant) =0. (20) 


At this point it is useful to make a further ap- 
proximation. We have carried the Bessel and 
Neumann functions as far as it is convenient— 
to carry them further leads to complicated ex- 
pressions. However the argument 2¢:2? is so large 
that they may be replaced by their asymptotic 
representations. The quantity 202+ will usually 
have values in excess of 20 and in this range we 
may use the following: 


Jo= (moz')—! cos(202!— 2/4), 
No=(xoz')—? sin(2c2'— 7/4), 
J, = (mo2')— cos(20z' — 32/4), 
Ni =(mo24)- sin(20z'—3/4). 


- (21) 
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Substituting these expressions in (20), and sim- 
plifying, one obtains 


A(cose—tanf sine) + B(sine+tanf cose) =0, 


where ¢=202!'— 2/4. Dividing through by cose 
this becomes 


A(i—tanf¢ tane)+B(tan¢é+tane) =0 
or 


A+B tan({+e) =0. (22) 
Now ¢ and ¢ are both functions of z and equa- 


tion (22) must hold for z=a and z=b. There- 
fore, we get the two equations 


A+B tan{{(a)+«(a)} =0| | 
A+B tan{{(b)+e()} =0 


These equations can only have roots other than 
the trivial ones A = B=0 if 


(23) 


tan{{(a)+e(a)} =tan{f(b)+e(b)}. (24) 
Equation (24) can be satisfied if 
{¢() +e(b)} — {F(@) +€(@)} 
=nr,n=0,1,2,--- (25) 


Introducing the functions ¢ and ¢ explicitly (25) 


becomes 
(=) | 
x a} 


(b}—a}) 146 
—tan | : (— )}. (26) 


where x=20(b!—a'). The permissible values of 
a are determined by finding the values of x for 
which this equation is satisfied (x is a function 
of a). Before looking at (26) in detail it is obvious 
that when 2=0 its solution is x=0. This means 
o=0or 


x=nr+tan—! 











a=}—2(a+)d). 


This will give a small negative value for ‘a. 
Remembering that out treatment cannot be ex- 
pected to be correct when a@ has a small absolute 
value it will be realized that this is the approxi- 
mation to the eigenvalue a=1. However the 
term in (5) arising from this eigenvalue (viz., the 
term in e*) is simple and there is no point in 
doing otherwise than taking it as it stands. 
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Fic. 1. The full line is the function M(—38.7, 1, zs) and 
the dotted line the function —0.4205 M(—38.77, z). 





Since b—a~™~1, we have 








1 
(b}—a!)= ~—. 
bi+ai 2b! 
Therefore, 
(b4—a!) (1+5) (-)}— 0.5 
x b ON x bx 
When n=1, 2, --- in’ (26) an actual computa- 


tion using representative figures shows that 
xcvnx. Thus the arguments of the inverse tan 
functions will be small, and these terms will be 
very small in comparison with nz. This being 
the case we will replace x by mz in the inverse 
tan functions and expand them retaining only 
the first term of the expansions. This gives 


: am =| 


x=nt + 
a‘ bi 


nT 1 fita 1+0 
[ite 141 oy 


at bi 





nT 
or 





 2(b'—a}) ne 


on 


It is interesting to see how accurately this for- 
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mula gives a; (i.e., »=1) for the problem re- 
ferred to above in which a=5.0 and 6=6.2. A 
simple computation gives o,;=6.154 or a, 
= —38.77. The correct value is a,= —38.70. 
Thus it is seen that for large absolute values of 
a the approximation is a very good one. Using 
the exact theory, several weeks of computing 
were required to get a, and the labor involved 
in getting a, seemed prohibitive. Formula (27) 
gives a result, correct to within 0.2 percent in 
this case, with only a few minutes required for 
the computation. When m= 2, we have o2 = 12.36 
and a,=—154.6. This checks with the exact 
theory as to order of magnitude and is probably 
quite accurate because for such large values of a 
the approximations are especially valid. This 
eigenvalue and the ones which follow it are so 
large that rapid convergence of the summa- 
tion in (5) is assured even for small values of r. 

Returning to a consideration of Eq. (14) we 
see that it should now be written 


M (ain, 2) =A ne?!? { Ji(2onz*) +K,Ni(20,2%)}, (28) 
where 
K,=B,/An=—cot{fn(a)+en(a)}. (29) 


Since it is not usual to find the J; and N, func- 
tions tabulated for such large values of the argu- 
ment as occur in this problem it seems best to 
replace them by their asymptotic representations 
in M(a,, 2). Then (28) becomes 


M (an, 2) =e7!22-* {cos(2on2' — 32/4) 
+K,, sin(2e,2'—32/4)}. (30) 


It will be noted that the A, as well as some other 
constants which occur as multiplying factors have 
been dropped. That this is permissible is at once 
evident from Eq. (5). M(as, 1, 2), or its approxi- 
mation M(ap, 2), is found twice in the numerator 
and twice in the denominator and any constant 
multiplying factor will cancel. Hence M(an, 2) 
as written in (30) will replace M(an,1,2) as 
given by the exact theory, and the a,’s will be 
given by (27) and (11). 

A noteworthy feature of the function M(ap, 2) 
is its simplicity and the ease with which nu- 
merical results can be obtained from it. Its 
validity is best illustrated by comparing the 
curve obtained from it with that obtained from 
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the exact function M(ap,, 1, 2). This latter func- 
tion is tabulated for 5.0=z=6.2 and a= —38.7.! 
In Fig. (1) the two functions M(—38.7, 1, z) and 
— 0.4205 M(—38.77, 2) are plotted for purposes 
of comparison. (It is to be observed that the 
two M functions are not identical—one is a 
multiple of the other.) No attempt has been 
made to find the best multiplying factor: the 
factor —0.4205 is the one that makes the two 
functions identical at z=5.0, the extreme end of 
the range. The usefulness of the approximation 
seems well attested. within the range of values 
for which it was devised. 

The evaluation of the integrals which occur in 
Eq. (5) does not seem of importance. Since 
M(an,2) must be computed in any case’ the 
simplest procedure is to evaluate these integrals 
by numerical integration. 

Two sets of values of ¢/co for several values of 
7 are to be found in Table II. In the first column 
of each group are the results given by the ap- 
proximate functions and in the second column 
those given by the exact functions. These latter 
values are obtained from reference 4, for the 
problem considered there. In computing the 
values of c/co from the exact solution the terms 
etc., were not used. This did not 
seem to introduce a serious error into the results 
because it was felt that these higher terms were 
insignificant even for such a small value of 7 as 
.015. However, no proof of this could be given 
and it is reassuring to discover from the approxi- 
mate solution that this omission is justified. If 
the term in a is included the numerical values 
most seriously affected are only changed by a 


in ao, a3, *** 
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TABLE II. c(z, r)/co. 











7 =.015 7 =.03 7=.05 r=.10 
2 A A B A B A B 
5.0 .726 .717 .632 .627 .569 .567 .524 .524 
5.1 .795 .784 .695 .689 627 .625 
5.2 .853 .841 .754 .747 .687 .684 .640 .639 
5.3 .900 .889 .809 .804 .748 .746 
54 .938 .926 .863 .857 .813 .810 .778 .777 
5.5 .968 .956 .916 .910 .882 .879 
5.6 .994 .983 971 .965 .955 .953 .944 .944 
5.7 1.021 1.013 1.030 1.026 1.037 1.035 
5.8 1.051 1.050 1.096 1.095 1.126 1.126 1.148 1.148 
5.9 1.094 1.099 1.174 1.177 1.228 1.229 
6.0 1.157 1.165 1.269 1.273 1.344 1.346 1.398 1.398 
6.1 1.241 1.254 1.382 1.389 1.476 1.479 
6.2 1.358 1.374 1.519 1.528 1.628 1.632 1.705 1.705 








Columns A contain the approximate values and columns B the exact 
values of the function c(z, r) /co. 


few units in the third figure beyond the decimal. 
(It should be stated that the values given in 
Table II for the approximate solution do not 
include the terms in ae,a3,---. They were 
omitted purposely in order to have the same 
conditions in both sets of calculations.) The 
greatest differences between the two sets of 
figures occur in the section headed 7r=.015. 
(This corresponds to a time equal to about 1/7 
of that required to reach full equilibrium.) But 
even here the approximation is within 1.2 per- 
cent of being correct over most of the range of z. 
For larger values of r the differences are much 
less. 

To obtain a set of curves from the exact solu- 
tion would take an experienced computer several 
weeks: using the functions derived in this paper 
the same curves could be obtained in six or eight 
hours. 
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The Motion of a Conical Coil Spring 


I. EpsTern* 
Department of Physics, University of Cincinnati, Cincinnati, Ohio 
(Received November 22, 1946) 


The elongations and natural frequencies of conically-shaped coil springs are derived theo- 
retically for several dynamic boundary conditions. Frequencies are verified experimentally on 
nickel and piano wire springs, attesting the reliability of the computations. 


1. EQUATION OF MOTION 


HE dissimilarities in the analytic expres- 

sions for the oscillatory motions of coni- 
cally-shaped springs on one hand and cylindrical 
ones on the other, arise as results of the vari- 
ability of mass per unit length in the direction 
of motion and of the inconstancy of the spring 
“constant” in the case of the conical type. If a 
constant mass per unit length of wire, m, is 
assumed, the total mass of a single complete 
turn of wire can be stated as a function of the 
diameter of the turn, or as a function of the 
diameter of the largest turn and the linear dis- 
tance from it. Let it be required to find the mass 
of a complete turn of wire at an average distance 
x from a coil of diameter Do (Fig. 1). The cir- 
cumference of the turn in question is r(D»—2ax), 
where a is the tangent of the half-angle at the 
vertex of the cone. Hence, its mass is m times 
this value. The “‘inertial force’’ (mass times ac- 
celeration) of a single turn of wire is then, 
m(Dy—2ax)m@E/dt?, where E is the displace- 
ment of the turn (elongation of the spring), 
measured from its unstretched position. The 
inertial force of an element dx must be equated to 
the net force necessary to give the element dis- 
placement E. The formula for the deflection of a 
single turn of wire (as measured by one 360° 
revolution of the coil around the longitudinal 








a=tang 








hem X—or} 


Fic. 1. Schematic diagram of spring 
(longitudinal section). 


* Graduate student. 
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axis of the cone) of loop diameter D and wire 
thickness d, under a load of P pounds, is 
f=8D*P/Gd‘, where f is the deflection in inches, 
D and d are measured in inches, and G is the 
torsional modulus of elasticity of the wire used.! 
Since the spring constant K is defined as the 
force per unit deflection, the conical spring “‘con- 
stant” is K=Gd‘*/8D*, and the force on an 
element of wire due to spring tension is: 


0 OE 0[Gd* dE 
F=u —|x Jes = -|- -- — lax 
OxL Ox OxL8D* dx 





It is advantageous to use as independent 
variable the diameter of a turn of wire, defined 
by the equation Z=D=D,—2ax. Thus, equat- 
ing the inertial force to the tension force, the 
equation takes the form: 


2amZ'@?E ®E 3 OE 1) 
a’Gd* at? aZ2? ZaZ 








A solution of the form E(Z, t)=w(Z)-T(t) is 
assumed, and substitutions are made back into 
Eq. (1). The variables Z and ¢ may be separated 
and two ordinary differential equations. arise. 
For the equation in ¢ we get: 


aT /dt?+T =0. (2) 


The general solution of this equation is T=T7o 
cos(At+ 6), where \ = 2zf, where f; are the natural 
frequencies of the oscillatory system to be de- 
termined by boundary conditions operating on 
the ordinary differential equation in Z. For the 
equation in Z we get: 

— 2m) i dw 3 dw (3) 

a°Gd'  ZwdZ? Z'wdZ 


1 National Metals Handbook (American Society for Steel 
Training, Cleveland, 1933), p. 990. 
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For ease in computation to be done later, the 
following changes are made at this point. Z is 
changed so that z=Z/Dpo. Thus, the range of the 
independent variable, rather than being from 0 
to Do, will now be from 0 to 1. Further, let 
a’Gd*/2xmD,° be denoted by a single bulk con- 
stant c. Equation (3) then reduces to: 





+ 0. (4) 
dz* zdz c 
Series solutions of Eq. (4) may be identified 
with solutions in Bessel functions of orders 2 
and — of the form: 


W(y) = (7/2)! {A J_2)3(y) + BJo,3(y) ie (5) 


where y=)z*/3yc. Finally, the general solution 
of Eq. (1), in terms of the independent variable 
z, is: 


E(z, t) = (A/2) $2? {A J_2/3(Az?) 
+ BJo;3(Az*)} cos(At+6), (6) 


where \/3Vc has been replaced by A. An in- 
vestigation will now be conducted into the par- 
ticular solutions of this equation to be used to 
satisfy specific boundary conditions. 


2. SPECIAL BOUNDARY CONDITIONS 
I. Both Extremities Fixed 


The particular solutions of Eq. (6) which sat- 
isfy the boundary conditions E(0, ¢) = E(1, t) =0, 
are: 


E(z, t) = (A/2) 4s? BJe)3(Azs) cos(Ai+ 4). (7) 


By taking 6=0, we assume that the motion 
starts from rest, that is, dE/dt=0 at t=0. In 
order to satisfy the boundary condition at the 
base, it is necessary that E(1,¢) in Eq. (7) be 
zero for all ¢. This can be accomplished by limit- 
ing the values of A to those which are roots of 
the equation J2;3(A)=0. The first two roots of 
this equation are A,=3.38 and A,=6.53, while 
other roots may be deduced by adding 7 suc- 
cessively to each former root.? Since A;=),/3c}, 
and \;=2zf;, then f; (the ‘natural frequencies of 
the system under assigned boundary conditions) 


2 Dinnik, Archiv d. Math. und Phys. [3], 18, 337-338 
(1911). 
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TABLE I. 
Si(experi- 
Root Aj Fi(c/sec.) mentally) 
Ni 1st 3.38 2.59 2.58 
2nd 6.53 5.21 
Piano wire 1st 3.38 2.53 2.55 
2nd 6.53 5.08 





will be f;=3c!A;/27, where the A; are roots of 
J23(A) =0. The constant ¢ may_be determined in 
one of three ways. First, by the defining formula 
for c: c=a’Gd‘/2rmD,*, using quantities which 
are deemed reliable for a, G, d, and m. Another 
method is to analyze the static case for an ex- 
perimental determination of Gd‘ and using values 
for the other quantities to the best of one’s 
knowledge. A third method is to solve for yc 
using one experimentally determined natural 
frequency. Once determined in this manner, c 
will retain the same value regardless of boundary 
conditions or initial distortions. It is best to 
arrive at a suitable value of c by using all three 
methods simultaneously, and arriving at con- 
sistent values for all of the constants included 
in ¢. 

In experiments performed in conjunction with 
this problem, two springs (one of nickel and one 
of piano wire) were employed. By the methods 
just outlined, the value of c adopted for the nickel 
coil was 2.75 and for the piano wire coil, 2.65. 
Thus, for the nickel coil, 


3(2.75)8A; A; 
628 1.26 





and for the piano wire coil, 
3(2.65)'A; =A; 
6.28 «1.29 





j 


Table I shows the results of this case. 


IIA. Vertex Fixed, Base Free 
In this case, the boundary conditions to be met 
are that E(0, ¢)=0, and dE/dx be zero at z=1. 
As in the previous case, in order to fulfill the 
condition at the vertex, A, (Eq. (6)), must be 
made identically zero. Thus, for a particular 
solution, we choose: 


E(z, t) = (A/2)!s?BJe;3(Az*) cos(AZ). (8) 


369 








+ Jy (a) 














LA," .92— 


i 








Ig) \ 


Fic. 2. Graphical solution for 
first ‘two roots, A, and Ag of 
Eq. (9). 
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dE/dx must be zero for z=1, that is, 
(Az*) Jo)’ (Az*) = —$Jo)3(Az*) for z=1. 


Hence, the equivalent expression to be satisfied 
is: 





4 Jo/3(A) 
J 5/3(A) =— oa, 
3 (A) 


The roots of this equation will determine the 
values of A which will satisfy the boundary con- 
ditions at the base of the spring. Equation (9) 
may be solved graphically from tabulated values 
of Jsys(A) and of 4/3[J2;3(A)/(A) ]. By inspection 
of Fig. 2, Eq. (9) is satisfied when A=1.92, 4.9, 

_+++,. Using the formulas developed in Section I, 
the natural frequencies may be computed. Table 
[I shows the results of this case. 


(9) 


IIB. Base Fixed, Vertex Free 


In this case, the boundary conditions to be 
met are that E(1, ¢)=0, and that dE/dx be zero 
at =0. It is advantageous at this point to write 
the general solution of Eq. (4) in the form: 


W(y) = (y/2)!{ A J_2)3(y) +BJas(y) }, 
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and noting that y=0 when z=0, y=A when 
z=1. Now, from the definition of J_23(y), it is 
evident that 


1 
lim[(v/2) 4 J_9/3(y) ] ee oe 
y—0 a(—2/ 3) 


where m(m) =I'(n+1). Thus, when y=0, E(z, ¢) 
= E(0, t) ¥0, as is required by boundary condi- 
tions, and the solution is not trivial. Inclusion of 
Jz3(y) in the solution is unnecessary to meet 
boundary conditions and B can be taken as 
zero, thus reducing our solution to: 


E(y, t) =(y/2)§A J_23(y) cos(At), — (10) 














TABLE Ii. 
fi(experi- 
Root Aj Fi(c/sec.) mentally) 
Ni Ist 1.92 1.53 1.53 
2nd 4.9 3.9 
3rd 8.1 6.4 
Piano wire Ist 1.92 1.49 1.51 
2nd 49 ° 3.8 
3rd 8.1 6.3 
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or, in terms of z: 
E(z, t) = [(A/2) iA ‘]2? J—2/3(Az*) cos(AZ). (1 1) 


We now aim to satisfy the condition that dE/dx 
be zero at z=0, using Eq. (11) as the’particular 
solution of Eq. (1). This requires the following 
relationship to be satisfied : 


32(Az*) J_2/3'(Az*) = —2z2J_2)3(Az*) for z=0, 
or the equivalent form: 


3z[ ( on 2/3) J—2)3(Az*) = Az®Jy/3(Az*) | 
Ss=— 22J—2/3(Az*) . 


Terms in J_2;3(Az*) are mutually annihilated and 
the following relationship remains to be satisfied : 


— 3Az4*J1/3(Az*) =0 for z=0. (12) 


But this is true for z=0, hence, Eq. (11) satisfies 
the boundary condition at the vertex. 

To satisfy the boundary condition at the base, 
E(1, t)=0, the choice of A must be restricted to 
roots of the equation J_23(A)=0. Hence, the 
boundary conditions of this case are satisfied by 
the particular solution (11) where A; are roots of 
J_2;3(A) =0. The first three roots of this equa- 
tion are A; =1.25, A,=4.43, and A3;=7.58, while 
other roots may be deduced by adding 7 suc- 
cessively to each former root.? After converting 
these roots to fundamental frequencies, the re- 
sults of this case may be stated in tabular form 


(Table III). 


IIIA. Weight Attached to Base, Vertex Fixed 


The boundary conditions in this case are that 
E(z, t) be zero at z=0 and that KdE/dx = M@E/ 
df? at z=1, where K =Gd‘*/82', the spring ‘‘con- 
stant’”’ in this problem, and M the mass of the 
attached weight. ? 

In order to satisfy the boundary condition at 
the vertex, A in the general solution of Eq. (1) 
must be taken as zero. Hence, 


(A/2)§Bz? Jo/3(Az*) cos(AZ) (13) 


satisfies the physical conditions imposed at the 
point of support. To fulfill the boundary condi- 
tion at z=1 (point of attachment of weight), the 
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TaBLe III. 
fi(experi- 
Root Aj fi(c/sec.) mentally) 
Ni 1st 1.25 1.00 1.04 
2nd 4.43 3.52 
3rd 7.58 6.04 
Piano wire Ist 1.25 .96 .92 
2nd 4.43 3.44 
3rd 7.58 5.88 








following equality must hold: 
aGd* 





{4 Jo)3(A) — 3AJ5y3(A) } = 9 McA*Jo/3(A), 


or, writing this equation in a form more easily 
accessible to computation, 


9Mc 4 Joj(A) 
{1-| |}. =Jsys(A). (14) 
aGd* 3 (A) 








Thus, the particular solution of Eq. (1) which 
satisfies the conditions of this case is Eq. (13) 
wherein A; are roots of Eq. (14). The defining 
equation for A;, (14), and hence, for the funda- 
mental frequencies of this system, may be solved 
by graphical means. The right member in Eq. 
(9) tabulated for Case IIA, is present here again. 
The function to the left of the equals sign in 
Eq. (14) is then built up as a tabulated function 
of A and is plotted on the same sheet as J5/3(A). 
The intersection of these curves gives the roots 
A;. Three different weights (60, 80, and 110 
grams) are used in each case. For the coils used 
here, it was found that aGd‘=5.22 for the nickel 
coil, and for the piano wire coil, a corresponding 
value of aGd‘=3.72. Figure 3 gives, by inspec- 
tion, the roots of Eq. (14) necessary for the ful- 
fillment of the boundary conditions of this case. 
These roots are then converted into natural 


frequencies, yielding the following summary 
(Table IV). 


IIIB. Weight Attached to Vertex, Base Fixed 


The boundary conditions in this case are that 
E(z,t) be zero at z=1, and that —KdE/dx 
=M2#E/d at z=0, where K=Gd‘/8z', the 
spring ‘‘constant”’ in this problem, and M the 
mass of the attached weights. 
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Fic. 3. Graphical solu- 
tion of Eq. (14) for six 
- values of 9.Mc/aGd'. 
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0.6 4.0 14 48 2z 


To fulfill the boundary conditions of this 


problem, it is necessary to employ Bessel func- 


tions of orders 3? and —?. The solution taken is 


the general one: 
E(z, t) = (A/2)!s? 
x {A J~2/3(Az*) + BJ»/3(Az*) } cos(AZ). (15) 


The following two equalities must hold for the 
solution of this case: 


aGd‘ 





(A/2)3{3A (Az*) J_o/s'(Az®) +3B(Az3) Joys’ (Az?) 


+2A J_2)3(Az*) +2BJo/3(Az’) } 
= 9 McA?(A/2)!z? {A J_o)3(Az*) + BJo3(Az*)}, (16) 
‘for z=0, and 
(A/2)§2? {A J_2/3(Az*) + BJo/3(Az*) } 
Xcos(M)=0 (17) 
for z=1. The condition of Eq. (17) reduces to: 
A J_2/3(A) + BJo/3(A) =0. (18) 


In order to reduce Eq. (16) it is advisable to con- 
vert the derivatives of Bessel functions into re- 
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26 20 3¢)—COtéC<CSN 


lated Bessel functions, and then to take limits 
of each term as z—0. Doing this, the condition 
to be met at the vertex becomes: 


aGd‘At 9 McA%4} 
e- -| ———_|1 (19) 
(0.90275) 44 (2.67893) A! 


Applying also condition (18), A; and B; can be 
eliminated, and the final ratio obtains: 


aGd* Js 3(A) 
| -| oo) js ma ie 
Mc J _2)3(A) 











where M has been converted into grams. The 
values of A;-which satisfy Eq. (20) are the ones 
which, when put into solution (15), will have 
both boundary conditions satisfied. Equation 


(20) may be solved graphically. J2;3(A) and 














TABLE IV. 
Mass attached 1st frequency 
Wire in grams Ai fi filexperimentally) 

Nickel 60 1.06 .84 84 
80 .96 .76 14 

110 .84 .66 .64 

Piano wire 60 .96 74 73 
80 85 .66 .64 


110 74 57 56 
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Mass attached ist frequency 














Wire in grams Ai fi filexperimentally) 
Nickel 60 .89 as 71 
80 82 .66 65 
110 74 .59 58 
Piano wire 60 81 .62 61 
80 mi 58 56 
110 .66 51 .50 





J_»;3(A) have been tabulated against A as inde- 
pendent variable. aGd*/ Mc will have a constant 
value, depending on the spring used and the 
mass suspended. Figure 4 shows the final de- 
termination of the roots of Eq. (20). It gives, by 
inspection, the roots of Eq. (20) necessary for 
the fulfillment of the boundary conditions of 


TABLE VI. 








. . . . : : 
nickel coil piano wire coil 





Ratio of ist root fi fi fi fi 
diameters Ai theory exp. theory exp. 
Case I (1: ©) 3.38 (see results of Case I) 


3.3 
1:10!, (1:2.154) 3.6 2.86 2.86 2.80 2.81 
1:83, (1:2.000) 3.7 2.94 2.94 2.88 2.89 
1:63, (1:1.817) 3.85 3.05 3.02 2.99 2.97 


this case. These roots may be converted into 
natural frequencies in the usual way. The fol- 
lowing is a tabulation of the results of the theory 
and their experimental verification (Table V). 


IV. Two Intermediate Points Fixed 


Next considered is the general case of the spring 
fastened at any two intermediate points. The 
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Fic. 4. Graphical solution of 





Eq. (20) for six values of gg 
aGd‘/ Mc. 
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TY) +4 
Jag (#v) 
” f , a 
.20 
i Fic. 5. Graphical solution of 








Eq. (24) when ratios of largest 
diameter of coil to smallest are 
as 1:10!, 1:2, and 1:63. 











eee (ks10) b 
—— (k-@) . 
—-—4 (4-6) ‘ 
































boundary conditions, stated mathematically, are 
that E(2’, )=0 for 2’ =a and for z’=b. There is 
no loss in generality if the widest upper boundary 
is taken to have the value 1. This can be done 
since in Eq. (1), ¢ is defined as a’Gd‘*/2xmD,', 
where Do is the diameter of the widest turn of 
wire. This is, in effect, thinking of the upper 
bound as the base of the spring. Thus, the 
boundary conditions to be used are: E(1, ¢) =0 
= E(b, t), where 0<b<1. We take as a solution 
of Eq. (4) the general solution : 


E(z, t) = (A /2) 172 {A J ~2)3(Az*) + BJo)3(Az*) } 
Xcos(A/). (21) 


To satisfy boundary conditions at z=1 and at 
z=b, the following equalities must be met 
simultaneously : 


A J_2)3(A) + BJoj3(A) ™ 0, 
A J ~2)3(Ab*) + BJo;3(Ab*) =(. 


(22) 
(23) 


The solution of this pair of simultaneous equa- 
tions depends upon the satisfaction of the fol- 
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lowing relationship : 
J-2j3(A) = Jnya(A) 
F-a(Ab*)— Jaja(Ab*) 





(24) 


Equation (24) may be solved graphically for 
specific values of b*. Values of J,(kA,) for n=3 
and —3, and for k=6, 8, and 10 are derived from 
tabulated values of J_2/3(x) and Ja3(x) by a 
shift of axis of the independent variable. It will 
be noted that this case degenerates into Case I 
when k becomes infinite (that is, when the lower 
diameter is taken as zero). Figure 5 shows the 
solution of Eq. (24) for values of b?=6, 8, and 
10. It gives, by inspection, the values of A; which 
will satisfy Eq. (21) and both boundary condi- 
tions. A; may be translated into f; and in the 
usual way. Table VI below gives the results of 
this case. 
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An Introduction to the Mathematical Theory of Plasticity 


W. PRAGER 
Brown University, Providence, Rhode Island 


(Received December 12, 1946) 


The paper aims at introducing the reader to the principal theories of plasticity. Since a pres- 
entation of the general stress-strain relations used in these theories would require too much 
space, the discussion is restricted to the mechanical behavior of plastic materials under shear. 
Theories of plastic deformation (Hencky, Nadai) and theories of plastic flow (Saint Venant- 
Lévy-Mises, Prandtl-Reuss, Prager) are illustrated by the example of a prismatic bar under 


torsion. 





1. INTRODUCTION 


N designing a structural member or a machine 
part the engineer, as a rule, uses formulas 
which are based on the theory of elasticity in a 
more or less rigorous manner. The highest stress 
predicted by these formulas is then compared 
with the maximum stress of the same type which 
the material can stand without failure. A so- 
called factor of safety is derived from this com- 
parison. Now, in all but the very simplest, cases 
the factor of safety obtained in this manner will 
differ considerably from the factor of safety de- 
fined as the ratio of the ultimate load to the de- 
sign load. The reason for this is, of course, the 
fact that the stresses set up by the ultimate load 
will exceed the elastic limit of the material and 
that, consequently, the theory of elasticity is 
unable to predict these stresses. In order to ob- 
tain the second factor of safety, the only one 
with a physical meaning, the mechanical be- 
havior of the material beyond the elastic limit 
must be taken into account. This behavior must 
also be taken into consideration in the study of 
manufacturing processes which aim at obtaining 
a favorable stress distribution under service con- 
ditions by setting up suitable initial stresses 
through previous plastic deformations. Another 
important field of application of the theory of 
plasticity is furnished by the phenomenon of 
creep. Under a constant load acting over a long 
period of time and at a high temperature a slow 
but steady increase of the deformation may take 
place. The problem of predicting the total de- 
formation due to creep which may occur during 
the lifetime of a structure is of an obvious im- 
portance. In the cases mentioned so far, the 
plastic deformations must be rather small, as 
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otherwise they would endanger the safety of the 
structure. This assumption of small deformations 
will no longer be justified in the study of tech- 
nological processes like rolling or drawing. How- 
ever, in an introduction to the theory of plas- 
ticity, it would certainly be unwise to try to 
cope at the same time with the difficulties arising 
from the consideration of finite deformations 
and from the unelastic behavior of the material. 
In the following it will, therefore, be assumed 
that the deformations, though greater than those 
possible within the elastic range, can still be 
treated as small. In order to show that such an 
assumption is reasonable, let us consider a bar 
of mild steel in simple tension. At the elastic 
limit the unit extension is of the order of one 
per mille. Now, we certainly are entitled to treat 
a unit extension of one percent still as small, but 
it is already far greater than any unit extension 
possible within the elastic range. 


2. ELASTIC TORSION 


The essential characteristic of the mechanical 
behavior of an elastic material is adequately ex- 
pressed in Hooke’s famous statement ‘‘ ut tensio 
sic vis." No similarly concise description of the 
mechanical behavior of plastic materials is pos- 
sible, and a mathematical theory attempting to 
take into consideration all mechanical phe- 
nomena observed in the plastic range would not 
be practical. A workable mathematical theory 
must idealize the mechanical behavior of struc- 
tural materials to a certain extent, and different 
idealizations are being used in the various fields 
of application. The fact that there exist almost 
as many theories of plasticity as writers on this 
subject is often a source of confusion to the stu- 
dent beginning to explore this field. 
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In the following, some of the more important 
theories of plasticity will be discussed briefly. 
Since a presentation of the general stress-strain 
relations used in the mathematical theory of 
plasticity would require too much space, the fol- 
lowing discussion is restricted to the mechanical 
behavior of plastic materials under shear such 
as occurs at a generic point of a cylindrical or 
prismatic bar in torsion. Using a system of rec- 
tangular coordinates x, y, z, with the z axis 
parallel to the generators of the cylindrical or 
prismatic surface of the bar, we let 7, and 7, 
denote the components of the shearing stress 
acting at a generic point of a cross section. We 
shall assume that all other stress components 
vanish and that 7, and 7, are independent of z. 

The equation of equilibrium,' 


07,/dx+07,/dy =0 (1) 


is fulfilled if the stress components are derived 
from a stress function (x, y) in accordance with 


T2=O0Y/dy, Ty= —Oyp/dx. (2) 


The condition that the cylindrical or prismatic 
surface of the bar is free from stresses requires 
that the stress function ¥ has a constant value, 
for instance zero, along the boundary of the 
cross section of the bar.! 

In order to obtain a differential equation for 
the stress function, we must consider the de- 
formations of the bar and make use of the stress- 
strain relations. As in elastic torsion, the dis- 
placement components are assumed to be of 
the form 


u=—yz0, v=x2z0, w=o¢(x, y, 6), (3) 


where 6 denotes the angle of twist per unit length, 
and the function ¢ represents the warping of an 
originally plane cross section of the bar.? The 
shear strains corresponding to the displacements 
(3) are’ 


V¥z2=9¢/dx—YO, yy=A¢/dy+x. (4) 


Up to this point, our discussion of the torsion 

1 See, for instance, S. Timoshenko, Theory of Elasticity 
(McGraw-Hill Book Company, Inc., New York, 1934) p. 
230. 

2J. Mandel, Ann. Ponts Chaussées 116, 1-33 (1946) 
assumes w=O6(x, y). Such proportionality between warp- 
ing and angle of twist can hardly be expected in the plastic 
range. 

+See, for instance, S. Timoshenko, reference 1, p. 233. 
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problem applies equally well to elastic and 
plastic materials. Further progress, however, de- 
pends on the introduction of a stress-strain rela- 
tion, and different materials require different 
stress-strain relations. For the purpose of the 
present discussion, the stress-strain relation of an 
isotropic material is best expressed as a relation 
between the vector « with the components 1;, 7, 
and the vector y with the components yz, Yy,, 
both vectors being considered as functions of 
time, as is the angle @. In addition to the vectors 
< and y, a stress-strain relation may therefore 
contain their derivatives with respect to time, 
d2/dt, dy /dt, 0°/dt, etc. In expressing the stress- 
strain relation in vector form, we’make certain 
that it is independent of the choice of the co- 
ordinate axes x, y, as it must be for an isotropic 
material. 

For an elastic material, the stress-strain rela- 
tion (Hooke’s law) has the following particularly 
simple form 

<«=Goy, (5) 


where Go is the shear modulus, a characteristic 
constant of the material. Introducing (2) and (4) 
into the scalar equations equivalent to (5), we 
obtain 
dy /dy =Go(dye/dx— ye), 
— dy /dx =Go(d¢/dy+x6). (6) 


Elimination of the warping function ¢ from Eqs. 


(6) furnishes the following differential equation | 


for the stress function: 
ay /dx?+0*y/dy? = — 2G, (7) 


where Gp and @ may be considered as given con- 
stants. Together with the condition that y=0 
on the boundary of the cross section, Eq. (7) 
completely defines the stress function y. 

As L. Prandtl‘ has pointed out, the function y 
satisfying (7) and vanishing along the boundary 
can be obtained by experiment. Let us consider a 
thin membrane, for instance a soap film, which is 
uniformly stretched over an opening in a rigid 
plane plate. The edge of the opening is given the 
shape of the contour of the cross section for 
which we desire to solve the boundary value 
problem formulated ‘above. Furthermore, the 
membrane is fastened to the edge of the opening. 


*L. Prandtl, Jahresber. Deutsch. Math.-Ver. 13, 31-36 
(1904). 
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If we now have a small excess pressure p on one 
side of the plate, the membrane will be deformed 
and it can easily be shown that the deflections 
w of the membrane out of its original plane will 
satisfy the differential equation 


0w/dx2+0*w/dy? = — p/T, (8) 


where JT is the uniform surface tension in the 
membrane. Along the edge of the opening we 
have w=0. Choosing the value of p/T in a suit- 
able manner, we can therefore obtain the surface 
z=y(x, y) by experiment. Equations (2) show 
that the magnitude of the vector ¢ is represented 
by the slope of the deformed membrane, and 
that the vector « is tangential to the contour 
lines of the deformed membrane.5 


3. STRESS-STRAIN RELATIONS FOR PLASTIC 
MATERIALS UNDER SHEAR 


In its initial phase, engineering research in 
plasticity was concerned with the limits of the 
elastic range rather than the phenomena occur- 
ring beyond this range. The condition that a 
given state of stress and strain is on the boundary 
of the elastic domain is expressed mathematically 
by the vanishing of a certain expression in the 
components of stress and strain (yield condition). 
For isotropic materials, this expression must be 
invariant with respect to rotation of the co- 
ordinate axes. Since Hooke’s law holds up to the 
boundary of the elastic domain, all strain com- 
ponents which might appear in the yield condi- 
tion can be expressed in terms of stress com- 
ponents. The yield conditions for an isotropic 
material then stipulates the vanishing of a cer- 
tain stress invariant. In the torsion problem, 
where the state of stress at a point is defined by 
the vector «, the only yield condition of this 
kind is 

|<| —k=0, (9) 


where || denotes the intensity of stress, i.e., the 
magnitude of the vector +, and k’the yield stress 
in pure shear. 

The yield condition (9) establishes the limit 
up to which the stress distribution can be ob- 
tained from the theory of elasticity. To carry 
the stress analysis beyond this point, we must 


5For a comprehensive survey of the elastic torsion 
oma see T. J. Higgins, J. App. Phys. 10, 248-259 
(1942). 
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explore the stress-strain relations in the plastic 
range. It was already pointed out above that, 
for the torsion problem, a stress-strain relation 
reduces to a relation between the vectors + and 
y and their derivatives with respect to time. We 
shall restrict the following discussion to stress- 
strain relations which do not involve time de- 
rivatives of an order higher than the first. More- 
over, we shall neglect all viscosity effects. This 
means that we shall only admit stress-strain 
relations which are homogeneous in the time 
rates of stress and strain. In spite of these re- 
strictions which are suggested by the results of 
experiments on metals such as steel, copper, and 
aluminum, we still retain considerable freedom 
in fitting stress-strain relations to experimental 
evidence. 

With each homogenous state of stress and strain 
in a plastic body there is associated a certain 
permanent strain, i.e., the residual strain after 
the stress is reduced to zero by unloading the 
body. That part of the strain which disappears 
upon unloading is called the elastic strain. In 
the case of the torsion problem, we shall denote 
the elastic and permanent strains by y’ and y”, 
respectively. Experience shows that the elastic 
strain is related to the stress by Hooke’s law: 


2=Gor’. (10) 


In establishing stress-strain relations for plastic 
materials, we must distinguish between loading 
and unloading. An infinitesimal change of the 
state of stress and strain in a plastic material is 
said to constitute loading if it involves a change 
of the permanent strain. On the other hand, 
during an infinitesimal change of stress and 
strain which constitutes unloading the per- 
manent strain remains constant. In the case of 
the torsion problem, we have dy/dt=0dy'/dt for 
unloading, and Eq. (10) suggests 


Ax /dt=Grdy/at (11) 


as the stress-strain relation for unloading. Equa- 
tion (11) is used in all theories of plasticity as 
the stress-strain relation for unloading; these 
theories therefore differ only in the stress-strain 
relations for loading. Moreover, even for loading 
the elastic strain is supposed to be related to the 
stress by Hooke’s law, so that the task of es- 
tablishing a stress-strain relation for a plastic 
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material consists essentially in defining the man- 
ner in which the permanent strain, or its rate of 
change, depend on the stress. For the torsion 
problem, in particular, we must define the man- 
ner in which y” or dy’’/dt depend on ¢. 

Let us first consider the case where y”’ is made 
to depend on ¢ (theories of plastic deformation). 
The isotropy of the material requires then that 
these two vectors have the same direction. This 
may be expressed by the equation 


Goy’”’ =(A—1), (12) 


where A is a scalar factor of proportionality 
which, as a rule, will be a function of the co- 
ordinates x, y, and the time. To obtain a proper 
relation between y’’ and « we must supplement 
Eq. (12) by a relation between the intensities 
ly’’| and |+|: 

|*| =Gof(|y’"|). (13) 


In the following, Eq. (13) will be called the flow 
condition. The flow condition enables us to 
eliminate the factor of proportionality, \, from 
Eq. (12). In the following, we shall consider two 
kinds of flow conditions: 

(a) for perfectly plastic materials: Gof=k 

=const., 

(b) for materials with strain-hardening: for the 
small permanent strains to which this dis- 
cussion is restricted, f is a monotonically 
increasing function of its argument; it is 
then convenient to write the inverse of 
(13) in the form 


Gol x’ | =|2|Le(14l /Go) —1]. (14) 


The elastic strain is given by Eq. (10), the 
permanent strain by (12) and (13). Combining 
Eqs. (10) and (12), we obtain the following rela- 
tion for the total strain y=y'+y”: 


In the case of a perfectly plastic material Eq. 
(15) together with the flow condition |+| =k 
constitute the stress-strain relations for loading. 
.This stress-strain relation must, of course, be 
supplemented by Hooke’s law, Eq. (5), which is 
supposed to be valid as long as |*! <k. Stress- 
strain relations of this type were first discussed 
by H. Hencky.® In the case of a material with 


*H. Hencky, Zeits. f. angew. Math. Mech. 4, 323-334 
(1924). 
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strain-hardening the form (14) of the flow con- 
dition enables us to eliminate \ from Eq. (15). 
Since both y’ and y” have the direction of ¢, 
we have 


Gol” | =GoL |x| — |v’| ]=Go| x| — |]. 


i 


Equating the last member of this continued 
equation to the right-hand side of (14), we find 


Go|! =|! (|| /Go). 


Comparing this to the equation obtained by 
taking absolute values on both sides of (15), 
we find 


A=g(|¢| /Go). 


The stress-strain relation for the loading of a 
material with strain-hardening therefore takes 
the form 


Goy =tg(|2| /Go). (16) 


Setting Go/g=G, we may write (16) in the 
form *=Gy, where G is a function of |*|/Go. 
This stress-strain relation which, of course, ap- 
plies to loading only, has the form of Hooke’s 
law, but the shear modulus G, the so-called 
secant modulus, depends on the intensity of 
stress. The relation «=Gy covers the elastic as 
well as the plastic range, if G is made to equal 
Go for |*|<k. The practical stress analysis is 
simplified, however, if a single analytic function 
G(||/Go) is used for both the elastic and the plas- 
tic range. To obtain good results one must, of 
course, choose a function which departs but little 
from the initial value Gp as long as || <k. On ac- 
count of its mathematical simplicity, this type of 
stress-strain relation which was introduced by 
A. Nadai’ is widely used in the modern litera- 
ture ;° it gives a continuous transition from the 
elastic to the plastic state. 

The stress-strain relations of Hencky and 
Nadai are sometimes criticized because they let 
the stress depend on the instantaneous strain 


only and not on the history of straining as might 


7A. Nadai, Plasticity, a Mechanics of the Plastic State of 
Matter (McGraw-Hill Book Company, Inc., New York, 
1931), p. 75. 

8 See, for instance, B. Finzi, Atti Ac. Sci. Torino, Cl. 
Sci. Fiz. Mat. Nat. 76, 222-238 (1941); F. P. Cozzone, 
J. Aero. Sci. 10, 137-151 (1943); A. A. Ilyushin, Prikladnaia 
Mat. Mekh. 7, 245-272 (1943); L. M. Kachanov, ibid. 6, 
187-196 (1942); E. Melan, Oesterreich. Ing.-Arch. 1, 14-21 
(1946). W. R. Osgood, J. Aero. Sci. 11, 213-226 (1944); 
W. W. Sokolovsky, J. Appl. Mech. 13, Ai-A10 (1946). 
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be expected for plastic materials. We shall re- 
turn to this criticism after discussing some rival 
theories of plasticity. 

Let us now consider the case where dy’’/dt 
is made to depend on ¢ (theories of plastic flow). 
The isotropy of the material requires that these 
two vectors have the same direction, and we may 
write 


Gud” /at =e, (17) 


where yu is a factor or proportionality similar to 
the factor \ appearing in (12). Combining Eqs. 
(10) and (17), we obtain the following relation 
for the total strain y=y’+y’: 


Goody /dt = de/dt+ pr. (18) 


| For a perfectly plastic material with the flow 
condition || =k or -<=k?, we have ¢-02/dt=0. 
Scalar multiplication of each side of Eq. (18) by 
« vields therefore 


u=Gor-dy/k?dt=GoW /k?, 


where W=+<-dy/dt denotes the rate at which 
work is done. Substituting (19) into (18) we 
obtain 


(19) 


Gody/dt =e /dt+GoWs/k? with |z| =k. (20) 


This stress-strain relation for loading must, of 
course, be supplemented by Hooke’s law, Eq. 
(5), which is supposed to be valid as long as 
«| <k. Stress-strain relations of this type were 
first discussed by L. Prandtl® and E. Reuss.?° 
As has already been indicated above, the 
sharp distinction between the elastic and plastic 
ranges with separate stress-strain relations com- 
plicates any stress analysis considerably. It 
therefore seems desirable to modify Eq. (20) so 
as to obtain a gradual transition from the elastic 
to the plastic state and take account of strain- 
hardening at the same time. As the present 
author" has shown, this can be accomplished 
without difficulty. In fact, let us drop the condi- 
tion |<| =k in (20). Scalar multiplication of each 
side of this equation by ¢ then furnishes 


GoW = 1S+G WS, ‘R?, 


(21) 
*L. Prandtl, Proc. ist Internat. Congr. Appl. Mech., 
Delft, 1924, pp. 43-54. 


1 E. Reuss, Zeits. f. angew. Math. Mech. 10, 266-274 
(1930). 


1 W. Prager, Proc. 5th Internat. Congr. Appl. Mech., 
Cambridge, Massachusetts, 1938, pp. 234-237. 
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where S=*4-t= ||? is a measure of the intensity 
of stress and the dot denotes differentiation with 
respect to time. Solving (21) for W and substi- 
tuting into (20), we obtain 


dy Ot S/k? 


Ginna msl cerimnanitninly, (22) 
at at 2(1—S/k*) 


Let us apply this stress-strain relation to the 
case of pure shear, where 7,=7,=0. Equation 
(22) then takes the form 





Tz7,/k? ts 
Go¥z2=t2+ T= ° (23) 
1—7,?/k? 1—7,7/k? 
With the initial condition 7r,=0 for y,=0, in- 
tegration of (23) furnishes 
=k tanh(Goy-/k). 


The graph of 7, vs. yz has the initial slope Go 
and the asymptote r,=k. The stress-strain rela- 
tion (23) thus provides a continuous transition 
from elastic behavior governed by Hooke’s law 
to the flow of a perfectly plastic material. 
Equation (23) is easily modified so as to in- 
clude strain-hardening. To this end we write 


r 

4 -( 1), 
S 

where G* denotes the so-called tangent modulus 


which depends on the stress intensity S. Indeed, 
specializing as in (23), we obtain 


Trtz 
Goy:= tz2+ (=- ‘> 


dr,/dy; = 


OY =. 


a =—+ 


(24) 





oh 


or 
G.* 


From a graph of 7, vs. yz the tangent modulus 
G*'can be found as a function of 7,7; the quantity 
G* appearing in (24) is the same function of S. 
The stress-strain relation (24) was suggested in 
1942 by H. J. Laning in an unpublished manu- 
script ; more general stress-strain relations of this 
type are discussed by G. H. Handelman, C. C. 
Lin, and W. Prager.” 

The first term on the right-hand side of Eq. 
(18) represents the rate of change of the elastic 


2 G. H. Handelman, C. C. Lin, and W. Prager, Quart. 
App. Math. 4, 397-407 (1947). 
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strain. For perfectly plastic materials, it has 
sometimes been suggested to neglect the elastic 
strains at the side of the often much larger per- 
manent strains. Equation (18) then takes the 
form 

dy/dt=y's with |<| =k. (25) 


Stress-strain relations of this type were the first 
to be used in the theory of plasticity, in the two- 
dimensional case by B. de Saint Venant," and 
in the three-dimensional case by M. Lévy™ and 
R. v. Mises.“ As will be shown below, serious 
objections can be raised against the neglection of 
elastic strains.'® In view of this fact, the amount 
of effort which has been spent on developing this 
particular theory of plasticity is surprising." 


4. PLASTIC TORSION 


In applying the stress-strain relations of the 
preceding section to the torsion problem, let us 
first consider perfectly plastic materials. These 
possess a clearly defined yield point below which 
they obey Hooke’s law (Hencky, Prandtl-Reuss) 
or are rigid (Saint Venant-Lévy-Mises) ; in the 
plastic range they flow under constant ||. Once 
the applied torque is sufficiently large to produce 
yielding, the cross’sectioncan be divided into 
elastic (or rigid) and plastic regions. In the latter 
we have || =: or, in terms of the stress function, 


(dy /dx)?+ (dp /dy)? = k?. (26) 


As the applied torque increases the plastic re- 
gions expand until, finally, they cover the entire 
cross section. 

Let us first consider this final stage of the 
torsion test. Here Eq. (26) is satisfied throughout 
the cross section. The left-hand side of this 
equation is the square of the slope of the tangent 
plane of the surface z=y(x, y) with respect to 
the plane z=0. According to (26) the surface 
z=y(x, y) therefore is a surface of constant slope. 


( a70) de Saint Venant, Comptes rendus 70, 473-480 
1 , 

4M. Lévy, Comptes rendus 70, 1323-1325 (1870). 

16 R. v. Mises, Goettinger Nachrichten 582-592 (1913). 

16 See also W. Prager, J. App. Phys. 15, 65-71 (1944). 

17 Excellent surveys of the Saint Venant-Mises theory of 
plasticity have been given by H. Geiringer, Fondements 
mathématiques de la théorie des corps plastiques isotropes 
(Gauthier-Villars, Paris, 1937),.and S. G. Mikhlin in the 
chapter entitled ‘“The mathematical theory of plasticity’’ 
of the book Some New Methods in the Mechanics of Con- 
—— Media (Russian), (Izd. A. N. USSR, 1938), pp. 
157-216. 
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We can produce such a surface experimentally 
by heaping sand on a horizontal tray which has 
the shape of the cross section. Indeed, the slope 
of such a sand hill has a constant value corre- 
sponding to the angle of friction, and the bound- 
ary condition y=0 is enforced by the shape of 
the tray. For a bar of quadratic cross section, 
for instance, the surface z=y(x, y) is a pyramid, 
for a rectangular cross section it has the form of 
a roof and for a circular cross section the form of 
a cone. 

In the case of a perfectly plastic material, and 
for a torque which produces plastic flow through- 
out the entire test specimen, the stress distribu- 
tion can thus be obtained without reference to 
any specific relation between stresses and strains. 
This fact is usually expressed by stating that, 
under these circumstances, the stress distribu- 
tion in torsion is statically determinate.'* 

So far, we have assumed that the torque is 
sufficiently large to produce plastic flow through- 
out the entire test specimen. For smaller torques 
only part of the material will become plastic. 
According to the theories of Hencky and Prandtl- 
Reuss, the stress function y will satisfy Eq. (7) 
[with (dy/dx)?+ (dy/dy)? <k*] in the elastic re- 
gion of the cross section and Eq. (26) in the plas- 
tic region. Continuity of the stress vector + at 
the boundary between elastic and plastic regions 
requires continuity of the first derivatives of y. 
The complete problem for the stress function has 
therefore the following form: to determine a func- 
tion W(x,y) with continuous first derivatives, which 
vanishes along the contour of the cross section and 
satisfies Eq. (7) with (dp/dx)?+(dp/dy)? <k? in 
part of the cross section and Eq. (26) in the rest.'® 
As Nadai'® has pointed out, Prandtl’s membrane 
analogy can easily be modified so as to enable 
us to solve this problem by experiment. Let us 
consider the square cross section as an example. 
A membrane is then stretched uniformly over a 
square opening in a rigid plate and fastened to 
the edge of this opening. Over the opening we 
erect a pyramidal surface of the slope which 
corresponds to the yield stress k. If now the 
membrane is deformed by slowly increasing pres- 

18 This wag first pointed out by A. Nadai, Zeits. f. 
angew. Math. Mech. 3, 442-454 (1923). 

'9For an ingenious manner of constructing particular 


solutions of this problem by an inverse method see V. 
Sokolovsky, Prikladnaia Mat. Mekh. 6, 241-246 (1942). 
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Fic. 1. Elastic and plastic regions in the cross section of a 
prismatic bar under torsion. 


sure on one side of the plate, it will at first 
assume a shape which corresponds to the stress 
function in the elastic case. As the pressure in- 
creases, however, the membrane will partly be 
pressed against the pyramidal surface. Where 
this is the case, the membrane represents the 
plastic stress function satisfying Eq. (26). Where 
the membrane is not in contact with the pyra- 
midal surface it represents the elastic stress func- 
tion satisfying Eq. (7). 

Having determined the stresses by means of 
this so-called soap film—sand hill analogy, we 
can determine the strains. Let us first apply the 
theory of Hencky, again considering a square 
cross section as an example. In Fig. 1, let the 
shaded areas be the plastic regions as obtained 
from the soap film—sand hill analogy and the 
unshaded cruciform area the elastic region. Since 
the stress function, and hence the stresses, are 
determined by the analogy, the strains in the 
elastic region can be found by applying Hooke’s 
law. As Eqs. (4) show, this means that, in the 
elastic region, the warping function ¢ can be 
determined to within an immaterial constant. 
To show how the warping function in the plastic 
region is found, we consider the shaded area 
ABC in Fig. 1. Here 7,=k, 7,=0, and the y 
component of Hencky’s stress-strain relation (15) 
gives 


vy =99/dy+xd=0, (27) 
while the x component does not furnish any in- 
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formation on account of the presence of the un- 
known factor X. Integration of (27) furnishes 


= —xyb+f(x). (28) 


Along the arc ABC the warping function ¢ is 
known from the elastic solution, and this en- 
ables us to determine the function f(x) in (28). 

Next, let us apply the theory of Prandtl and 
Reuss to the same problem. The soap film—sand 
hill analogy shows that, for any given torque, 
this theory furnishes the same stress distribution 
as that of Hencky. Moreover, it can be shown 
that the warping function, too, is the same for 
both theories.”° Indeed, let us consider a generic 
point P of the cross section. As the torque is 
increased, the boundary between the elastic and 
plastic regions sweeps across the section of the 
bar. Once P has been engulfed in the plastic 
region, the stress ¢ at this point remains con- 
stant. For instance, if P is engulfed in the shaded 
area ABC of Fig. 1, the stress at P will have the 
constant components 7r,=k, r,=0. Accordingly, 
02/dt is zero, and the y component of the stress- 
strain relation (18) of the Prandtl-Reuss theory 
furnishes y,=0, while the x component does 
not give any information on account of the pres- 
ence of the unknown factor uw. The strain com- 
ponent y, therefore retains whatever value it has 
at the instant when the boundary between the 
elastic and plastic regions sweeps across P. On 
account of the continuity of the stress vector, 
however, the stress components at P just before 
this point is engulfed into the plastic region must 
be r,=k, tr, =0. At this instant we therefore have 
¥,=0 at P, by Hooke’s law. Since from this 
instant on y,=0 at P, we have y,=0 at a generic 
point P inside the shaded area ABC of Fig. 1. 
From this and the identity of the stress dis- 
tributions for the two theories follows the iden- 
tity of the warping functions. 

The theories of Hencky and Prandtl-Reuss are 
thus seen to furnish two identical predictions in 
the case of the torsion problem. This result must 
not be unduly generalized, however. As a rule, 
the predictions furnished by the two theories 








2% This was first pointed out by the present author in 
H. Geiringer and W. Prager, ‘“Mechanik isotroper Koerper 
im plastischen Zustand,”” Ergebnisse d. exakt. Nat. wiss. 
13, 310-363 (1934), p. 331. 
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differ; conditions under which they agree were 
investigated by A. A. Ilyushin.” 

As regards the theory of Saint Venant-Lévy- 
Mises, it does not give a solution of the torsion 
problem except in the case where the torque is 
sufficiently large to produce plastic flow through- 
out the entire test specimen. Indeed, the stress 
distribution in the plastic part being the same as 
was studied above, Eq. (25) gives 


Vy=9¢/dy+xb=0 
or 


¢= —xyb+f(x) (29) 


for a point in the shaded region ABC in Fig. 1. 
Since the cruciform unshaded region now is sup- 
posed to be rigid, we have ¢=0 along the arc 
ABC. This condition enables us to determine 
f(x), and hence ¢, inside the shaded area ABC 
provided that its boundary, the arc ABC, is 
known. The theory of Saint Venant-Lévy-Mises 
does not furnish any condition, however, from which 
this boundary between the rigid and plastic regions 
could be determined. Moreover, the stress distribu- 
tion in the rigid part is indeterminate. 

Let us review the results obtained so far. The 


Saint-Venant-Lévy-Mises theory is applicable 


only if the torque is large enough to produce 
plastic flow throughout the entire bar. For 
smaller torques the theories of Hencky and 
Prandtl-Reuss furnish identical results. In all 
but the very simplest cases, however, the theo- 
retical determination of the boundary between 
the elastic and plastic regions encounters prac- 
tically unsurmountable mathematical difficul- 
ties.” These difficulties are avoided by the theo- 
ries of Nadai and Prager which stipulate a con- 
tinuous transition from the elastic to the plastic 
state. 

Let us first apply Nadai’s theory to the torsion 
- problem. In terms of the stress function the 
intensity of stress, on which the shear modulus 
depends, is given by 


dy\? ay\ 27) 
| | -|(—) +(~) = |grady|. 
Ox oy 


#1 A, A. Ilyushin, Prikladnaia Mat. Mekh. 9, 207-218 
(1945). 

2 For a numerical treatment of this problem by: relaxa- 
tion methods see D. G. Christopherson, Trans. A.S.M.E. 
62, Al-A4, and F. S. Shaw, Austral. Counc. Aeronaut. 
Rep. ACA-11 (1944). See also R. V. Southwell, Relaxation 
methods in theoretical physics (Oxford Press, 1946). 
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The stress-strain relation «=Gy therefore gives 


1 dy d¢ 


——=—-y), ——-—=—+x0 


ye, 
G dy Ox G dx dy 


Elimination of gy between these equations 


furnishes 
1 dy 1 dy 
—)+ —(— —)= 26 (30) 
G OxJ/ dy\G dy 


which equation replaces Eq. (7). Here G is a 
function of |grad y|, and the stress function y 
vanishes on the contour of the cross section. 
Equation (30) is non-linear, and somewhat 
reminiscent of equations which are encountered 
in connection with steady two-dimensional flows 
of compressible fluids. Integration methods de- 
veloped in fluid dynamics might therefore prove 
useful; so far, however, the boundary value 
problem formulated in connection with Eq. (30) 
has been solved only in a few special cases.** 
Finally, let us apply the stress-strain relation 
(24) to the torsion problem. To this end, we set 


Sj Go ds 
{; = 1) = log (5). 
9 2\G*(s) $ 


This enables us to write Eq. (24) in the form 


OY a 
Go— = A—(*/A). (31) 
ot ot 


Introducing the warping function and the stress 
function into the two scalar equations equivalent 
to (31) and eliminating the warping function be- 
tween them we obtain** 


0 0/1 dv\7 
156%) 
OxL Ot\A dx/ J 


ar 1 dy a0 
+— Pid -—)|- —2Ge—. (32) 
dyL at\A ay at 





This equation contains the time rates of change 


* The as yet unpublished Ph.D. thesis of G. H. Handei- 
man (Brown University, 1946) is concerned with the 
application of this theory to thin-walled sections. See also 
L. A. Galin, Prikladnaia Mat. Mekh. 8, 307-322 (1944). 

* This equation was first established in W. Prager, 
“Theory of plasticity,” (mimeographed lecture notes) 


Brown University, Providence, Rhode Island, 1942, p. 112. 
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of y and @ in a homogeneous manner. Accord- 
ingly, the time scale is irrelevant and the angle 
of twist 6 may be taken as measure for the time. 
Thus, 


0 0/1 dy 
“p5(02) 
OxL O0\A Ox 


af a1 oy 
+—| 4 (— ~)]- —2Go. (33) 
dyL a0\A dy 


Considered as an equation for y(x, y, 6), Eq. 
(33) is non-linear. It is worth noting, however, 
that this equation is linear in dY/0@ with co- 
efficients which depend on y. This remark sug- 
gests the possibility of integrating (33) approxi- 
mately by solving a number of linear problems 
for dy 00 with coefficients which depend on the 


function y built up from the previously obtained 
values of dy/00. Physically speaking, this pro- 
cedure amounts to producing the desired angle 
of twist in a series of small steps, and assuming 
that during each step the tangent modulus at 
any point of the cross section retains the value 
which corresponds to the stress intensity exist- 
ing at the beginning of this step at the point 
under consideration. 

For perfectly plastic materials, the theories of 
plastic deformation and of plastic flow furnished 
the same solution of the torsion problem; except 
for the circular cross section these theories will 
yield different results, however, when strain- 
hardening is taken into account. Torsion of a 
cylindrical bar of non-circular cross section may 
therefore afford an opportunity of deciding ex- 
perimentally between these two theories. 





On the Theory of Noise in Radio Receivers with Square Law Detectors* 
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Department of Mathematics, Cornell University, Ithaca, New York 


AND 
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For the video output V of a receiver, consisting of an i-f stage, a quadratic detector, and a 
video amplifier, the probability density P(V) has been obtained for noise alone and for noise 
and signal. The results are expressed in terms of eigenvalues and eigenfunctions of the integral 


equation 


[°K ®p(s—b)f (dt =V(s), 


where p(r) is the i-f correlation function (i.e., the Fourier transform of the i-f power spectrum) 
and K(t) is the response function of the video amplifier (i.e., the Fourier transform of the video 


amplitude spectrum). Two special cases are discussed in which the integral equation can be 
solved explicitly. Approximations for general amplifiers are given in the limiting cases of wide 
and narrow videos. Some applications of the method to other problems are shown in Sections 


7B and 9. 


I. THE RECEIVER MODEL 


HE response of a linear circuit to noise has 
been investigated by several authors. The 


* The major part of this paper is based on work done 
for the Office of Scientific Research and Development under 
contract OEMsr-429 with Cornell University and contract 
OEMsr-262 with Massachusetts Institute of Technology. 

** Now at the Physics Department, Syracuse University, 
Syracuse, New York. 
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probability densities and correlation coefficients 
of the noise outputs of some types of detectors 
are also known.' To find the probability density 
of the noise voltage output of the next stage 


! References to the work in this field can be found in the 
articles by S. O. Rice, Bell Sys. Tech. J. 23, 282 (1944); 
24, 46 (1945), and M. C. Wang and G, E. Uhlenbeck, 
Rev. Mod. Phys. 17, 323 (1946). 
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Fic. 1. The receiver model. 


(video or audio amplifier) is a problem of a dif- 
ferent type, since after the detector, the noise 
voltage is no longer Gaussianly distributed. 

In this paper we have, therefore, calculated 
the probability density P( V) of the noise voltage 
V at the end terminals of a receiver consisting 
of an intermediate-frequency amplifier, a quad- 
ratic detector and a video (or audio) amplifier 
(Fig. 1). The response of both amplifiers is as- 
sumed to be linear and characterized by the 
steady-state response function y(f) and C(f), 
respectively, or the response functions Q(t) and 
K(t), such that the output voltages of the ampli- 
fiers are y(f)e"“* and C(f)é**, respectively, if 
the input voltage is e*'. If, therefore, the video 
input voltage is v(t), the video output voltage V 
is given by 


vio=f agciper f dt'e?*'st’y(t’) 


—z 


; (1.11) 
= f dt'K(t—t’)v(t’) 


where 


K(x)= fi dscipyer (1.12) 


Since V(t) can depend only on the values which 
v(t’) assumes at times ¢’ =, it follows that K(r) 
=0 for r<0, which means that all poles of C(f) 
lie above the real axis. Since K must be real, 
we have C(—f)=C*(f). Both conditions are ful- 
filled for passive networks. We shall normalize 


‘K by 
f K(t)dt=1. 


In writing i-f voltages Vi.", and V1,r. for 
input and output it is generally convenient to 
consider them as modulated carriers, i.e., we 
write 

Vir.Y= Re {u(t)e?*#o*} 


= x(t) cos2rfot+y(t) sin2rfot, 
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where 
u(t) =x(t)—iy(t), 


and fo is the carrier frequency. The i-f output 
voltage is then 


Vir.O= Rel f dfe**'/'y(f) 


x f dtetir (terre | 


= Re| ei f dt'u(t’) 


xf dferris-re-eny(9)| 
= Re} U(t)e***/o"} 


where 


vo= f dt'Q(t—?’)u(t’), 


and 


an)= f dfe?*i/ myify= f df'e**'l"*( f’ + fo). 


We shall normalize y and Q by demanding that 
f QG(r)dr=1; 


then we have 


f Ornddr= ff farapasrerorsr ry s+fo 
7 fn Xy(f" +fo) =1, 


or 


J df'y(fotf y(fo—f’) =1. 


The calculations are appreciably simplified by 
assuming a symmetrical pass band* for the i-f 
amplifier, i.e., 


V(fotf) =7*(fo—f’). (1.20) 


We also put |y(fo+f)| =B(f). Inserting this in 
the normalization equation we have 


firoar=1, 


i.e., the power spectrum is normalized. The sym- 


* See, however, footnote 3, Section IV. 
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metry of y implies that Q is real and Vj.r. can 
now be written as 


Vir. =X (t) cos2arfot+ Y(t) sin2wfot 


with 


(1.21) 


x= f dt'Q(t—t’)x(t’) 
_ (1.22) 
v= f dt'Q(t—t’) y(t’). 


Q(r) must of course vanish for + <0. 
The detector is assumed to yield the output 


X*(t)+ Yd), (1.31) 
if the input voltage is 
X(t) cos2afot+ Y(t) sin2afot. (1.32) 


In order that this be a reasonable model of an 
actual detector it is necessary to assume 


|x’ (t) ly’ (t) | 
i—|<fo <fo, 
) | y(t) | 





and 
| x(t 

which implies that the width of the pass band of 

the i-f amplifier is small compared to the inter- 

mediate frequency fo. 


II. REPRESENTATION OF NOISE 


The i-f output noise voltage (omitting the ir- 
relevant d.c. component) is represented? as 


- /N me , et 
Vier. =— Dlrv(he')| 
VT kar 


x (X;,’ cos2xf,’t+ Y,’ sin2mf;,’t) (2.1 1) 


with 


fi’ =k/T, 


where X,’, Yi’, X2’, Yo’: --are independent nor- 


mally distributed random variables each having 
the probability density 


nm texp(—u’). 


The representation holds only for the time in- 
terval 0<t<T, but this is adequate if T is chosen 
large compared to the time during which ob- 
servations are made. We shall later let T-@. 
The properties of y(f) were discussed above and 
N denotes the input noise power per unit 
frequency. 


2S. O. Rice, ‘Mathematical analysis of random noise,” 
Bell Sys. Tech. J. 23, 282 (1944), see p. 328. 


VOLUME 18, APRIL, 1947 


According to (1.31) and (1.32), the detector 
output voltage arising from this i-f input voltage 
is given by 


X")+Y2)= Vo, 





where 
F i 5 
X(t) oe p> l-y(fe’) | (Xu! cos2a( fe’ — fot 
a +Y,' sin2«(f,’ — fo)t] (2.12) 
¥(t)=— S| y(fr’) |[—Xe’ sin2a( fi’ —fo) 
VT km 


+Y,’ cos2x( fi! —fo)t]. 


Making use of (1.20) and the fact that the band 
width of the i-f amplifier is small compared to 
the intermediate frequency fo, we may write 


V/N 2 
X(t)= VT > Bf.) (X, cos2afit+ Yi; sin2axf;it), 
ve . (2.13) 
Y(t) “VT = BC fi) 
xX (—X;, sin2rfi.t+ VY, cos2rf,t), 
where 
fr=k/T, B( fx) = | v(fet+fo) |, (2.14) 


Xp=X'k+ko VY. = Y'k+ko 


and 


with ko= 7 fo, assumed to be an integer. 
The video output voltage V(t) is given by 
(1.11) as 


V(t) -f K(t—’)(X7(/) + V2") Jd’ 
~ (2.15) 
-f K(u)[X?(t—u) + ¥?(t—u) du. 


In the derivation of Chapter 3 we shall consider 
the probability distribution of the slightly modi- 
fied quantity 


T 
ve)= f K(u)[_X2(t—u)+ Y*(t—u) |du, (2.16) 


and then pass to the limit as T> ~. 


Ill. “WHITE NOISE” 


The i-f output X(t) cos2afot+ Y(t) sin2afot 
with X(t), Y(t) given by (2.13) can be described 
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as caused by the “white noise” 


x(t) cos2afot+ y(t) sin2rfot, 


where x(t) and y(t) are given by 


‘N 
x(t) —. = (Xz cos2afut+ ¥, sin2a frt) 
k=—L 
VN os 
y(t) =- VT za (— 2% sin2rf,.¢+ VY, cos2rf;,t), 


where L is a number so large that the frequencies 
fox L/T lie outside of the pass band of the i-f 
amplifier. This‘ white noise”’ is the input thermal 
and shot noise which actually is known to have 
a constant power spectrum up to very high 
frequencies. 

In the derivation of the next section it will be 
advantageous to represent x(t) and y(t) in a dif- 
ferent orthogonal system, and we will now show 
that the expansion coefficients in any other com- 
plete orthogonal system are Gaussianly dis- 
tributed and independent of each other, if L is 
sufficiently large, i.e., if the power spectrum is 
constant up to sufficiently high frequencies. 

Let ¢;(¢) be an orthonormal, complete set in 
the interval (0, 7). Expanding x(f) in this system 
we have 

. 


x(t)=>0 uj39,(0), (3.1) 


1 
where 


N\i L T 
1;= (=) ps (xf o,(t) cos2rf,tdt 
7 bom, 0 


T 
vif ¢;(t) sind fut) (3.2) 
0 


The coefficients u; are thus linear combinations 
of the Gaussian variables X;, Y; and are, there- 

‘fore, Gaussianly distributed. To prove inde- 
pendence of linear combinations of independent 
Gaussian variables, it is necessary only to show 
that (uu;)=0 for 1Aj. Using the statistical 
properties of X; and Y;, we have 


uu j/m 5 th t rs) j (S 
l A l 


L 
X>> cos2rfi(s—t). (3.3) 


—L 
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For sufficiently large L 
1 L 
— > cos2rf,(s —t) 
F wk 
can be replaced by 6(s—?#) and, therefore, 
| N 
(Hitt jw jl. (3.4) 


In the same way we can show, for the co- 
efficients v; of y(t), that 


uel 3.5 
(vm jn j (3.5) 
7 l 
and 

(401) yO. (3.6) 


IV. DERIVATION OF PROBABILITY DENSITY 
FOR NOISE ALONE 


To derive the probability density P(V) we 
use the above representation for the incoming 
‘white noise’ writing the i-f input voltage* as 

Vy.r. =x(t) cos2rfot+y(t) sin2drfot, (4.11) 
with 
x(t) => ; u;¢;(t) 
y(t) => 5 230,(0). (4.12) 


and 


Choosing N as unit of power we have 
(1,1) 0 = (UV) wv = 45. (4.13) 
The i-f output voltate Vi.r. is given by 
Vir. =X (t) cos2rfott+ Y(t) sin2wrfot, (4.14) 


where 


X(t) = f . Q(t—1)x(1r)dr, 
| = (4.15) 
Y(t) -{ O(t—7)y(r)dr. 


’ The authors are aware that the results of this section 
can be obtained in a more elegant way by using the 
representation (4.12) for Vi.r.“ rather than for the 
components x(t) and y(t), and by defining the quadratic 


detector by 
Vo =(Vi.r. P 


rather than by Eq. (4.31), provided that the I.F. and 
video band widths are small compared with fo [cf. Phys. 
Rev. 70, 449, C6 (1946)|. The simplifying assumption of 
symmetrical I.F. pass band can then be easily dispensed 
with. The method presented here was chosen to keep the 
continuity with the existing literature. 
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The correlation coefficient 
p(t) =(X(t)X(t+7))w/(X2(t))y (4.21) 


can be expressed in terms of the function Q as 
follows: 


(XOXt+ w= f foU-—mMee+r—9) 


XX (X (71) X (72) ad rid r2 


1 x 
--{ O(t— 71) Q(t+7—71)dr1 


—@® 


1 DL 
--{ Q(8)Q(0+ 7)dé, (4.22) 


1 0 
(2) w= 5 f Q?(0)d0=3, (4.23) 


and therefore 


n)=f | Q(0)Q0(0+ 7)dé. (4.24) 


We note that 


[ eronae 
-{ J O(A)Q(6+ r)e?*""dédr 


=f ff ferarurtinver+so 


K err ils'O+s"' O+r) +i d adr 
=f fararrir tims +i0 
Ko +f’ 6f" +f) 
which is a special case of the Wiener-Khintchine 
relation. 
The detector output voltage Vp is, therefore, 


Vp =X*(t)+ Y*(0) 


= f fow- ri) Q(t— T) 
X [x (41) x( 72) +y(71)y(72) |dridtz, (4.31) 
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and the output of the video amplifier 


V(t)= J "K(t—6)[X*(0)-+ ¥(0) do 


- f “K(t—0d6 J fow-re0-r 


b 4 [x( 71) x( T2) +y(71)9(72) ‘|drid7z. (4.32) 
Changing variables to * 
6 =t—0, 1/=t—1m, Te’ =t—T2, (4.33) 


we have 


vo=f K(o' do’ ff Or 0) QC! -0") 


X [x(t— 71')x(t—y2’) 
+y(t—71')y(t—T2’) \dri'dro’. (4.34) 


If the i-f input noise voltage is a homogeneous 
random function x(t—7) and y(t—7) are given by 


x(t—7) = 20; uj'o,(7), 
y(t— T) = > 5 2/o;(7), 


(4.35) 


where u,;’ and 9,’ are random variables with the 
same properties as u; and v;, and we have— 
omitting the primes— 


XD i(uajtva)deiri)dst2)dndrz (4.36) 
= Diwan; +o.) 


xf fac T2)Oi(71)Oj( T2)d r1d T2 (4.37) 
with 
A(11, 72) -f K(6)d0Q(r1 —O0)Q(12—8). (4.38) 


We note that ¢ does not appear anymore on the 
r.h.s., that means that V(t) is a homogeneous 
random function. This is caused by the fact that 
the input noise is itself a homogeneous random 
function. 
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The essential step in our calculation is now to 


choose as the function ¢;({7) the normalized eigen- - 


functions of the integral equation 
f A(1, t2)o(T2)dt2=do(71). (4.41) 


Then V(¢) is expressed as 
V=D i(u?+v,*)ri, (4.42) 


where the quantities \; are the eigenvalues of the 


integral equation. The quantities 
u?Z+v7=Z; (4.43) 


are independent and their probability densities 
W(Z;) are given by 


e~7i if Z;>0, 


WZ)=\9 ig Z:<0, 


(4.44) 
because u; and v; are independent Gaussian vari- 
ables of spread 3. The probability distribution 
of the sum representing V is therefore 


pvydvaavf +f av-x.nz) Il. W(Z,)dZ, 
0 0 


1 a x 2 
=d V— f dee" f . f L]e e77*#-™*dZ, 
2r —© 0 0 





1 as 1 
=dV— f dée~‘¥§_____., (4.45) 
2a J_« TIe(1 —7Ac€) 
For V>0O the path of integration is closed in the lower half plane and contracted around all poles 
¢,= —1/X, with \,>0 and vice versa. We thus obtain 


e~ ive eV As) 





k#s 


P(V)= an a i = y » 
=e ae 1 “1i—it) “ 
. (—iny (=. ) 





for V>0 (4.46) 


wI'(1-~) 
” ote As 


where }>*+ means summation over all values of s for which \,>0. For V <0 we obtain 


1 dé 
P(V) =— 2 





aes : " (: 1 
(—4 n ee 
iXn 


where >_,,“~ indicates summation over all for 
which \, <0. The product is extended always 
over all k with exception of k=n. 

It is sometimes more convenient to have the 
integral equation in a different form, involving 
_ the functions K and p rather than K and Q. This 
form is obtained by introducing the new unknown 


f(0)= { Qr2-6(r2)dra, (4.51) 


so that 


Multiplying both sides with Q(7,:—6’) and in- 
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e~ive e(—ViAn) 
eel => - (4.47) 
‘(1-74 &) n XE 
AE Dak An (1 -~) 
kn n 








tegrating over 7; we obtain 


f K(o)d0 f d7,Q(71—@’) 


X O(71— 49) f(0) =Af(0’). (4.53) 


We can now write 


f Q(r1-0/)Q( 71-8) 74, 


-{ drQ(r)Q(r+0’—0) = p(6’—8), (4.54) 


and have the integral equation in the form 


J K(6)p(0 — 0) f(0)d0=rf(0’), 
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or 


br 
de 


Wi 


ru 
gi 


Ww 


Ww 


ar 








or, since K(@) =0 for 6<0, 
[ K@o(e'-9 s040=r400) (4.55) 


Finally, by writing (K(@))!f(@)=(6) we can 
bring the equation into the form which will be 
derived directly in Section 6: 


K (6’) }'p(0’ — 0)[_K(@) }'b(0)d0 = r¥(0’). 
J, Kore OK OP Od= 90). 
By means of the expression 

V=d diz; (4.61) 


we can easily verify formulas for the moments. 
For the first two moments we have for instance 





(Vyw=E MZw=E = f " K(s)ds=1, (4.62) 


(V2) ,=> AP(Z 7?) wt dD AAdu(Z Zi) 


j#k 
=> A?v+(Z Aj)? (4.63) 
1 
=] +> d;. 
For the fluctuation 
2= (V2 )w—(V) a. (4.64) 


We thus obtain 


e=> nt= ff K()KWo%(s—vasas, (4.65) 


this being the trace of the first iterated kernel. 


V. DERIVATION OF THE PROBABILITY DENSITY FOR SIGNAL AND NOISE 


For the case of signal and noise only one derivation will be given here, since the direct derivation 
runs parallel to the analogous derivation for noise alone. The i-f input voltage V1.r.“" is in this case 


given by 


Vier.“ =(x(t) +a(t) | cos2rfot+[y(t)+8(f) | sin2xfot, (5.11) 


where x(t), y(¢) are the input white noise amplitudes, and a(t), B(t) are the input signal amplitudes. 


The i-f output voltage is then given by 


Vier. =(X()+ p(t) | cos2afot +L Y(t) + q(t) | sin2rfot, (5.12) 


where X and Y are defined as above and 


p= f Q(t—r)a(r)dr and at)= f QO(t—1)B(r)dr (5.13) 


are the signal amplitudes after the i-f amplifier. The detector output voltage is then 


Vo =([X()+p() P+LVYO+¢) ? (5.14) 


and the output of the video amplifier is 


Vi)= f K (t—0)d0 {CX (0) +p() P+LYV(0+9(0) 7} 


=-{ Ku-oaef [ Q(8—71)Q(8— r2) {LX (71) Fa(r1) Jx(72) +a( 72) J 


—2 


Changing variables to 


+[ (71) +8(71) Ly(72) +8( 72) J} dridrs. (5.15) 


6’ =t—8@, 7) =t—n7, T2’ =t—Te, (5.16) 
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we obtain 


vio= Kode f f O(71' — 8’) O( 12’ — 0’) { [x(t — 71’) tat — 71’) x(t — 12’) +eal(t— 72’) J 


+[y(t—71') +8(t— 11’) [v(t — 12’) + B(t— 12’) ] }d ri/d 72! 


- 
zx 


= i) A(ri, t2)d rd ro | [x(t— 71) +a(t— 7) J[-x(t— 72) +a(t— 72) | 
, +[y(t— 71) +8(t— 11) v(t—72) +8(t—12)]}. (5.17) 
We now expand noise and signal amplitudes in terms of the eigen functions ¢; of A, writing 
x(t—7r) =D; ujo(7), a(t—7) =D; a,()o<(7), ; 
y(t—r) => ;v)¢/(7), Blt—7r) =D; B()o 7), ta 
and obtain 


V(t) =>; Aj Las tea,(t) ]?+[0j;+8,) }*}. (5.19) 


The probability density of the variables u; and v; are r~) exp(—a;*) and 27 exp( —2;), respectively, 
and the characteristic function for V(t) is, therefore, given by 

; exp[ — (u;?7+2;7) ] 
(ce) y= f- f exp {ié > AL(ujt+a,;)?+(@;+8,)7]} [] ————————adu dv; 


] T 


1 ® 
=I du dv; exp { —[(u?+0,°)(1 id) — 248 (u jj +08;) ]+itdj(a? +87) } 
- jdV j CX] 


1 ¥ | 1£X 50; ° 1EX 5B; *| 
== ff duqac exp| (=u in)! - ) = (24 -ieay'—— ; ) 
ind JS "i : " (1-itd))! sO * (1 ieay 7 J 


- 7 £*(a;7+8;*) ; . J ‘ 
exp] ———————— + it (a? +87) |, (5.21) 
1—7£); 








exp| 7£A,;(a;?+6;7) /(1—7&A;) 
(et¥), =T] expLié Aas Ps J (5.22) 
j 1—2td; 





The Fourier inversion integral of this expression which represents the probability density cannot be 
evaluated in closed form. 

We shall check this expression by deriving the known * probability density for the voltage after a 
quadratic detector. This distribution is obtained in the limit of infinitely wide video band width, 
which implies that 


f _K(6)¥(0)d0=(0) (5.31) 


for any function (8). Then A(r172) becomes simply the product Q(71)Q(72) and the integral equation 
reduces to 


Or) f Q(12)6(72)dt2=dG(71). (5.32) 
The only eigenfunction of this integral equation is of the form CQ(71) where C is a constant, and from 


*S. A. Goudsmit, RL 43-21; K. A. Norton and V. D. Landon, Proc. I.R.E. 30, 425, Sept. 42; J. C. Slater, RL V-23; 
D. G. Fink, Report T8; D. O. North, RCA Tech. Rep. PTR6C. 
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JS07(r)dr =1 we have C=1 and \=1. The characteristic function, therefore, reduces to 





exptt(a? +B”) /(1 —7€) 
1—it 


where 


a= f alti—)Q(r)de= f 


e- | B(t— r)Q( ndr= f B(7)Q(t—7')dr’ =q(0). 


a(r’)O(r—7')dr’ = p(t), 


Calling the signal after the detector s?=s*(t) we have 


and 





P(V)=— 
2r 


6) 





1 s-i- 
Oe 


for V>0 and zero otherwise. With u=1—7é this becomes 


1 0, 
P(V) =exp(— V—s?)— f 
2m1 


4 (5.33) 
S(t) =p()+97(t) =a°+8", (5.34) 
1 * exp(—7&V +72ks?/(t—7 
[> é a 8) ae 
- (1—12€) 

xp(—1&V+12£s?/(1 —7)) 

Se (5.35) 
1—1é& 
(+ 
exp[ Vu+s?/u |du/u 

=exp(— V—s*)Jo(2isy/ V) for V>0 and (5.36) 


zero otherwise, where Jy is the Bessel function of order 0. This result is in agreement with the earlier 


results. 


VI. DIRECT DERIVATION FOR NOISE ALONE 


In order to give a direct derivation of the re- 
sults obtained in Section 4, we investigate the 
probability density of the function 


T 
vio=f K(u)[x*(t—u)+y2(t—u) ldu (6.11)* 
0 


and let, at an appropriate instant, 7 approach 
infinity. For the sake of simplicity, we shall re- 
strict ourselves to the case K(u)>0. We divide 
the interval (0,7) into a large number 2 of 
equal subintervals and consider the expression 


7 n 
— 2D [*(t—uj) +9°(t—u)) JK (u)); 
Nn j- 
= (6.12) 
ee Ra 
T<u;<-—T. 
n n 





* For convenience, x and y are used in this section to 
denote the quantities formerly denoted by X and Y. 
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The characteristic function of the joint prob- 
ability density of 





X:=x(t—u1), X2=x(t—ue2), *+°, Xn=x(t—Uu,) 


defined as the average 


(exo > Ex(t— «)|) 


is easily found to be 


pr(u)—m) (6.13) 


i k=l 


1 nm 
exp} -- > 
4 Zh 
where the correlation coefficients pr for the finite 
interval T are defined as 
is 
pr(t) “T x. A(fx) cos2n f;t, (6.14) 


and 


A(f)=B*(f) = |y*(fafo)|. (6.15) 
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Thus the joint probability density of 
Nip Le, XB, ° °°, Xn 


is _ by the formula 


xe ol - 3 pr(uj;—uy,)&; sie. dé, 
| 44 
aL | 6.16) 
— —— exp ux jXee, (6. 
~ (y/m)" /D a a “eae 


where the matrix [a, | is the inverse of the cor- 
relation matrix [pr(u;—u,.)] and D the deter- 
minant of the correlation matrix. 

We next find the characteristic function of 


: i n 
— > K(u)x*(t—u,), 
nm 
i.e., the average 
eT * 
(exp| = ~ K(u,)x?(t— 
n 


), 


We have 


tT n 
(ols > K(u)x*(t—u) ' 
n Ay 
—f- ff exp; — 
mses VD —» 


xexp| - ails -+dx,. (6.21) 
1 


1¢T n 
ye K(uss | 





Since K(u;)>0 we can introduce the new 
variables 
xj =_K(u;) ]'x;. (6.22) 
Obtaining 
Cex i > K(u;)x*?(t—u;) ). 
“7 - = 4 TICK )} 
* " 1éT 
xf f exp|~ Eas) 
Xexp; — > Baris’ ban’ --dx,', (6.23) 
1 





where ; 
aires (6.24) 
” [K(uj) }LK(u)}} 


Let wi™, pe , un™ be the eigenvalues of 
the matrix [8 ]; then changing the coordinates 





(n) 
, 


in the last integral to ‘principal axes’ of the 
“ellipsoid”’ 
p> B jex jx), = | (6.25) 
1 


and noting that }°x/? remains invariant under 
the transformation, we get easily 


1¢T n | 
(exp| = > K(u,)x?(t—u,;) ) 
n 1 Ay 


1 
—— ee Or (6.26) 


VD TICK (u)} (ui) 


n 





Since we have assumed that A(f) is an even 
function we can show without much difficulty 
that 


> K(ujx*(t—u;) and > K(uj)y?(t—u,;) 
1 1 


are independent in the statistical sense of the 
word. Thus the characteristic function of 


(6.27) 
n 


— > K(u,) {x*(t—u) +y2(t—u,)} 
1 


1 
- : itT 
D]I K(u;) (4 -=) 
1 1 n 


we now observe that the 
verse of the matrix 


[(K(u;))}pr(uj— uy) (K(ux))*] 
and that the determinant of (6.28) is 





matrix [jx] is the in= 


(6.28) 


n 


D II K(u;). 
1 
Denoting by Ai, A2™, , \,™, the eigen- 
values of matrix (6.28) we see that the y’s are 


the inverses of the \’s and that 


D II K (uj) =), “Ao ae ‘> 
1 
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These facts imply that the characteristic func- 
tion of (6.27) is 


1/T1(! -—1,08). (6.29) 


It follows from Hilbert’s treatment of integral 
equations that the eigenvalues of matrix (6.28) 
when multiplied by 7/n approach, as m ap- 
proaches infinity, the eigenvalues of the integral 
equation 


J [K(s) Por(s—)[K() Vf (dt=rf(s), (6.31) 


provided these eigenvalues are simple. Since the 
distribution of (6.27) approaches that of V(t) we 
deduce that the characteristic function of V(t) 
is given by the formula 


1 T 1—72A;(78), 6.32 
/ 1 inj) (6.32) 


where 47, A", , are the eigenvalues of 
the integral equation (6.31). For large T we can 
replace the integral equation (6.31) by 


f | LK(s) }'e(s—)LK(t) }'f()dt= f(s), (6.33) 


and the characteristic function of V(t) by 
1/[] (1 —7),£), (6.34) 


where i, Ae, As, «+: are now the eigenvalues of 
the integral equation (6.34) which we again as- 
sume to be simple. The probability density P( V) 
of V can now be expressed by means of the in- 
version formula 


1 . oe 
P(V)=— { ———az. 6.35 
aa J. nase C' 


At the end of Section 4, the integration has been 
carried out by the method of residues. The justi- 
fication of the formal integration is not quite 
simple since one must investigate the behavior 
of the entire function 


¢(Z) =[][(1—7A,Z) 
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for large |Z|.5 However, it can be considerably 
simplified if one assumes that for some »(0<» <1) 
the series 

LAs 

1 
converges. On the other hand it is easy to see 
that the exact probability density P(V) is the 
limit, as m approaches infinity, of P,(V) where 

— ® exp(— V/s) 
P,(V)=2 


s=] 





; V>0. (6.36) 


if 


Xs []’(1—Aj/d,) 
1 
P,(V)=0; V<0O since with the assumption 
K(u)>0 negative eigenvalues cannot occur. In 
fact, the characteristic function corresponding 
to P,,( V) which is seen to be 


- 1 
i i —1£; 


approaches the characteristic function associated 
with P(V) thus P,(V) approaches P(V). 


VII. EXPLICIT SOLUTIONS IN SOME 
SPECIAL CASES 


1. Single-Tuned I.F. with Simple Low Pass Video 


For the special case of a single tuned i-f 
amplifier with the video amplifier acting as a 
simple low pass filter the integral equation can 
be solved explicitly and the eigenvalues can be 
expressed in terms of roots of Bessel functions. 
For one special ratio of video band width to i-f 
band width it is even possible to obtain P(V) 
in closed form. For other ratios it is not too 
difficult to sum the series for P(V) numerically. 

The correlation function p(t) is in this case 


p(t) =e-@'#! (7.11) 
which is equivalent to writing 


A(f)= Const. (7.12) 
1+(2nf/a)? 


5 We note that ¢(Z) =D(iZ) where D(z) is the Fredholm 
Determinant of the above integral equation. It is note- 
worthy that this relationship between the characteristic 
function of P(V) and the Fredholm Determinant persists, 
even in the general case in which K(t) is not assumed 
positive (in which case the simple product representation 
need.not hold, since 2A need not converge). In the general 
case one must take the Fredholm Determinant of the 
kernel K(t)p(s—t). 
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For the video we assume 


,u>O0 


a™)= 8 exp(—Bu), «<0 


1 
C(f) =—————. 
1+2mif/B 


The integral equation (6.33) now assumes the 
form 


(7.1 


wm 
—_ 


af e~ale—tle~ (8/2) +8) F(thdt=Xf(s). 


y(t) ~exp(F) f00 


* N 

f exp(—a|s—t)|) exp(—Bt)y(t)dt =—-(s). (7.17) 
3 
f 


Putting 


(7.16) 


we have 


This integral equation is solved by® 


2a\} 
¥(s) =Joais 2( ) e~ Bs | (7.18) 
Br 
and the eigenvalues are determined by 
J (2a 8) 4[.2(2a Bd)*]=0. (7.19) 


Thus the éigenvalues of our integral equation 
are expressible in terms of roots of Bessel func- 
tions and are seen to be simple. It is interesting 
to note that if 2a/8=% we can express P(V) in 
terms of one of the theta functions. In fact, if 
2a/8=%, Eq. (7.19) assumes the form 





Ji[(6/d)*]=0 (7.21) 
and since 
2\! sinx 
a} /Xx 
‘we see that 
6 1 
An=—— (n=1, 2, ---) (7.23) 
rn? 
_ and formula (6.36) for V>0 becomes 
© ©@ s* exp(—r’s?V/6) 
P(V)=— ft — —— (7.24) 
6.1 []’(1—(s?/j7*)) 
I 


* For the solution and discussion of the integral equation 
see M. L. Juncosa “An integral equation related to Bessel 
functions,’ Duke Math. J. 12, 465-471 (1945). 


394 





It can be shown quite easily that 


s? cosms a 
Il’{ 1-— j= -— (7.25) 
j - 2 
and therefore, 
P(V)= - > (—1)**'s? exp(—2*s?V/6). (7.26) 
2 ont 
Introducing the theta-function notation 
6(q¢; 0) =1—29g+2q'—2q°+--- (7.27) 
we have 
g0’(q; 0) = —2 & (—1)**'s2q" (7.28) 
s=1 
and finally 
: r rV als 
P(V) = —— exp Jere rV/6-()). (7.29) 
6 6 


For other values of the ratio a/8 we can cal- 
culate values of P(V’) and plot the corresponding 
graphs. 


2. The Noise Power, Averaged over a 
Finite Time Interval 


In some problems it is of interest to know the 
distribution of the average 


1 T 
M= i) I*(t)dt. 
T d/o 


The derivation of Section 4 is clearly applicable 
to this case if we put 


(7.31) 


a \1/7, 0<u<r 
~ |0 elsewhere. 


K(u) (7.32) 
The probability density is given by (4.46) re- 
placing V by A/ except that the \’s are now eigen- 
values of the integral equation 


1 T 
i) p(s—t) f(t)dt=Af(s). 


T</0 


(7.33) 


In the simple case 


p(t) =e-a'#! (7.34) 


the eigenvalues can be shown to be given by 
the formula 


a 


Xm = - | 
g 1+7,.” 





(7.35) 
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where y,, is the mth (positive) root of the equation 


tan gy= —2y/(1—y’) (7.36) 


and g=ar. 


VIII. APPROXIMATIONS FOR P(V) 


At the end of Section 5, it was shown that, for 
infinitely wide video, P(V) becomes identical 
with the well-known probability density of the 
output voltage of a quadratic detector. We will 
derive now an approximation for the case that 
the video band width is very large compared 
with the i-f band width. 

It can be shown that in this case one of the 
eigenvalues say \;, tends to unity while all others 
approach zero. In order to obtain the dependence 
of the eigenvalues on the band width ratio we 
shall write 


aa i 
A(f’)= (5) (8.11) 
B \B 
and f-f 
C(f) =c —), (8.12) 
where ‘ . 
f a(v)dvy=1 (8.13) 
and “ae 
c(fo/b) =1 (8.14) 


in conformity with our original normalization for 
A(f) and K(t). The limiting process is carried 
through for B/b=e-—0, while the functions a(v) 
and c(v) remain unchanged. In terms of A and C 
our integral equation becomes 


| A(fiCUfe—fig(fadfo=delfi) (8.21) 


and traces S; are, therefore, 


S,=1, 


Si= ff dndra(eryates) 


— 2 


x c(€[ 2 — |— fo/b)c( el — v2 | — fo/d), 


(8.22) 
Si= fff ardrsdra(esyatv:)a(rr) 


Xc(elve—v1 |—fo/b)c( el v1 —v3 | —fo/d) 
Xc(eLvs—ve]—fo/b), 
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etc. Expanding c(ev—f)/b) in a power series 





around fo/b we shall write 


c(ev— fo/b) =1+Cre+Cre?v?+---. 


Assuming that 
f a(v)v'dy 


is finite we obtain the following approximations: 
So>1+ae?+0(e'), 
S31+3ae?/2+0(), 
Ss1+2ae?+0(e'), 


(8.23) 


(8.24) 


where 


a= 2(2C.—C3') f a(v)v*dv (8.25) 


and essential use has been made of the fact that 
a(v) =a(—v). Noticing that 


AS2> Sy (8.26) 


(since S,=>oy dx”, and dy is by defintion the 
largest eigenvalue) we obtain 


rAr> Ss So = 1 +ae’+0(e*). (8.27) 


Therefore, )>7d,2 is of order ¢«* or higher, and 
> yA.‘ is of order & or higher. Now we notice 
that 


(S,)# =i + 0(é) (8.28) 
so that 
hi, =1+ae?/2+0(e). (8.29) 
Making use of the expression 
on NZx (8.31) 
1 
(cf. Eqs. (4.42), (4.43) ], we write 
"=Vi+0 (8.32) 
with 
Vi:=1—Ai+A\Zi (8.33) 
and 
v=) Ax(Z.—1). (8.34) 


If Pi(Vi) is the probability density of Vi, 
given by 


P,(Vi)=0 for Vi<i—vAy (8.35) 


and . 
Vi-—(1—Ay) 


P,(V,) =A, exp— for Vi>(1—A,), 


1 


and w(v) is the probability density of v, we have 





395 





for P( V) the exact expression 


Pv)= dvP,( V —v)w(v) 


V—I1+A, 
-{ dvP,(V—v)w(v). (8.36) 


Since the width of w(v) is of order € we can use 
P,(V) as the first approximation of P(V). To 
discuss the reliability of this appproximation we 
expand P(V) in the form 
V—-14 

vdvw(v) 


P(V)=P,(V) -Py(v) f 


« 


V—1+A\1 
+3P,''( vf v*dvw(v) 


1 V—14+\1 


4+ — 


T vdyw(v) P,’""(V —&) 


(8.37) 


where, in the last term, 
0<é< |p. 


For any fixed value of V>0, the integrals in the 
first and second term can be replaced by 1=0 
and (v’)=>°>A,2 with an error of order ¢ or 
smaller, because 


f vdvw(v) 
vV—(1—dy) 


<f- : w(v)dv=O0(e), (8.38) 
o LV-(1-)} 


f vdvw(v) 
V—(1—A)) 


rs) v' 
<f [Va yy erm ). 


The last term in (8.37) can be shown to be of 
order ¢® or higher. Thus we have 


P(V) =Pi(V) + 3P 1" (V) vw +0(€*) (8.39) 








-for any fixed V>0. Substituting 
hi =1+ae?/2+0(e) 


and 


(0° )m= 2 ae? =O(e'), 


k=2 


we see that for fixed V>0, P,(V) is a good ap- 
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proximation except for terms of order ¢* and 
higher. 
If, in addition, the existence of 


f a(v)v®*dy 


is assumed the approximation can be improved 
by including the term 


2P1"(V)(0")m 


(8.41) 


where 


(0* wv = S2— (Se)'+-0(e'?) 
and A, [which occurs in P;( V) ] is calculated from 


Ai = (S4)'+0(€). (8.43) 


(8.42) 


Substituting these values we obtain an approxi- 
mation for P(V) valid up to terms of order e* 
inclusive. In the neighborhood of V =0, no ap- 
proximation which utilizes only the traces of a 
limited number of iterated kernels, can be ex- 
pected to hold, since there—as can be seen from 
the general solution—the eigenvalues _indi- 
vidually determine the shape. 

The considerations of this section cannot be 
expected to be applicable to the case treated 
exactly in Section 7 since there not even 


[core 


exists. In that case one sees directly that the 
eigenvalues are such that 


1=1-0(-.), ¥d2=0(e). 


In order to be able to obtain a general approxi- 
mation for wide videos one must demand that 
the i-f spectrum falls off sufficiently rapidly at 
infinity. 

In the opposite limit, b/B-—0, similar con- 
siderations show that P(V) becomes Gaussian 
with mean (V),=1 and spread proportional to 


(b/B)}. 
IX. APPLICATION TO A DISCRETE CASE 


For some applications of noise theory it is 
desirable to have the characteristic function 


exp id fhe) =$(Ey-*+&) (9.11) 


k=1 
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of the joint probability of m successive values 
I, In-++, In; TL=x?(t) +y°(t.) of the detector 
output, from which the averages of products 
(U1 eli"*)w (m, integer) can be derived. $(£1- + - En) 
can be obtained by carrying through the idea of 
transformation to principal axes for the discrete 
case, that is, by a slight modification of Section 
6. However, we shall derive $(£,---£,) from the 
results of Section 4 since this seems the quicker 
way. 
We simply chose for K(t) the function 


1 n 
K(t) “ LD &5(t—t) (9.12) 
1 


where &=)(7& for the sake of normalization, 
and the time variable ¢ is chosen in such a way 
that all ¢, <0. For this choice of K(t) we have 


a) 1 n 
V= f KOLMO + Om E tele (9.13) 


and, provided that the J, are normalized so that 
(Tw =1, we have 


b(E1° + En) Seimtele = (eit V),, = (9.14) 


IT i(1 —2a;€) 
The integral equation, whose eigenvalues are \,; 
degenerates in this case into a system of linear 
equations for m discrete variables f(t.) =fi: 


1 n 
: DL ExPemfe=Afm, (9:15) 
kool 
where 
Pkm = p(t, — head 
with 


A= £X. 


We write these equations as 


> ExPimbk = Afn. (9. 1 6) 
k=1 
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To obtain ¢, it is not necessary actually to de- 
termine the eigenvalues since the product 
TJ ,(1—7A,) can be expressed in terms of the de- 
terminant det | £.pim-+t5im| as follows: 

We note that 


n 


I] (tu —iA,) is a polynomial in u which has zeros 
1 
at «=A; and starts with i"u". The determinant 


det | E:Pim —USim| is a polynomial in u, which has 
the same roots and starts with (—)"w". We there- 
fore have 





I] (iu —1d,) =4-* det | Ek Pim — USkm| (9.17) 
1 


and, specially, for «= —i 


II(1 — iA;) =1-" det | i Pim +15 km | A 
I 
The characteristic function is thus obtained as 
4” 


(ei=kelk),, = > ° 
det | Expim+75km | 





(9.18) 


While it has not been possible to evaluate the 
joint probability 


1 
a) a nd Pore ee 
(i++ Ts) rood | fas : 
Xo(Er + EeHM, (9.19) 


explicitly (except for n=1 and n=2, in which 
cases they are well known) one can obtain aver- 
ages of the form 


(1x Zk) for integer values of m,, 


by expanding both sides of (9.18) in power series 
in the variables £, and by comparing coefficients. 








7 For the general kernel a corresponding relation can be 
derived between the above product and the Fredholm 
determinant of the kernel. 
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The analogy between thermal and electrical networks is used to obtain solutions to problems 


relating to the heating and cooling of air-backed and solid-backed bolometer elements exposed 
to radiation of periodically fluctuating intensity. From the equations derived, thermal response- 
frequency characteristics are computed for “typical” bolometer elements. Features of these 


characteristics are discussed, and the physical influences of the various thermal properties of the 
bolometer material and of the backing material are pointed out. 


INTRODUCTION 


N some of the present-day applications of radi- 

ation bolometry the heat-sensitive receiver 
element is placed in series with another resist- 
ance, and a d.c. biasing voltage is applied across 
this series. The radiation incident upon the 
element is mechanically interrupted, as by a 
rotating sector disk, so that the element alter- 
nately heats and cools. Under these conditions 
an a.c. voltage is generated across the element, 
which, when amplified and measured, serves as 
an indication of the intensity of the incoming 
radiation. 

The amplitude of the a.c. signal thus generated 
is approximately proportional to the biasing 
voltage across the receiver element. This, in 
turn, is limited by the rate at which the element 
can dissipate the Joule heat produced by the 
bias. If the element is placed in thermal contact 
with a solid backing through which this heat 
can be more rapidly dissipated, the bias voltage 
can be increased. The amplitude of the tempera- 
ture fluctuation of the element due to the 
absorption of a given chopped radiation signal 
is reduced by the presence of such a backing, 
but this reduction, for certain ranges of chopping 


frequency, can be more than compensated for by 


the use of the higher biasing voltage made 
possible by the backing. In addition, the use of 
a backing brings about a faster approach to 
temperature equilibrium at each interruption of 
the radiation beam, so that the time constant 
of the element is correspondingly reduced, a con- 
sideration of importance, for example, in the 
rapid scanning of an infra-red spectrum. 

There are a number of questions in connection 
with the backing of bolometer elements for which 
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both experimental and theoretical answers have 
been sought. It is of importance to know, for 
example, what thermal properties and what 
geometry the backing should have in order to 
secure maximum dissipation of the bias heating 
with minimum reduction of the amplitude of 
the temperature fluctuations of the element. It is 
equally important to know the effect of such a 
backing on the frequency dependence of the 
element. 

It is the purpose of this paper to demonstrate 
a method for obtaining answers to such questions 
by applying the mathematics of electrical net- 
work theory. The electrical analogy has long been 
useful in the solution of problems in mechanics 
and acoustics, and some applications have already 
been made in the field of heat flow.'~* Analytical 
expressions governing the behavior of bolometers 


STRIP BACKING 




















u(i+e™) RZ WS 


Fic. 1. Physical arrangement of bolometer strip, 
backing, and holder. 


1V. Paschkis and M. P. Heisler, J. App. Phys. 17, 246 
(1946). 

2C. D. McCann and H. E. Criner, Machine Design 17, 
137 (1945). 

3 V. Paschkis and H. D. Baker, Trans. A.S.M.E. 64, 105 
(1942). 

4M. Avrami and V. Paschkis, Trans. A.I.C.E. 38, 631 
(1942). 
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Fic. 2. Electrical analog of the bolometer 
assembly of Fig. 1. 


will be developed from this analogy, and these 
will be employed in presenting graphical illustra- 
of the temperature fluctuation of the 
element in its dependence on the thermal proper- 
ties of the bolometer and the backing, and on 
the chopping frequency. 


tions 


General Equations for Backed Bolometers 


Consider a thin strip of temperature-sensitive 
bolometer material of 1 sq. cm area and Ax 
thickness. Consider this strip to be fastened to a 
solid slab of backing material, also of 1 sq. cm 
area, and extending in thickness from x=Ax to 
x=L. At x=L let the backing slab be fastened 
in turn to some kind of holder or housing. This 
physical arrangement is sketched in Fig. 1. 

Let the exposed front surface of the bolometer 
strip absorb an incident radiation signal of 


Uo(it+e*)* (1) 


watts, where Up is half the radiant power per 
unit area in the unchopped beam, w is 27 times 
the chopping frequency, and ¢ is the time. Let 
heat be dissipated by reradiation and convection 
from the front face of the strip, according to 
Newton’s law of cooling, at a rate of o(T,~0—T,) 
watts, where o is the dissipation constant of the 
front face in watts per sq. cm per degree, T,<o is 
the instantaneous temperature of the front face, 
and 7, is the ambient temperature. Similarly let 
heat be lost from the rear surface of the backing 
slab to the holder at a rate o’(T,-1—T,) watts, 
where o’ is the backing-to-holder dissipation 
constant (here considered to be thermally non- 
reactive) and T,., is the instantaneous tempera- 
ture of the rear face of the backing. These various 
heat flows are indicated in Fig. 1. Consider the 


*In (1) there is considered, for simplicity’s sake, only 
the fundamental component of the multiply-harmonic 
radiant power input obtained by the usual mechanical 
chopping of the incident beam. 


VOLUME 18, APRIL, 1947 





lateral surfaces of the bolometer and backing to 
be thermally insulated so that heat can flow 
only in the x coordinate. 

Let k and k’ be the thermal conductivities of 
the bolometer material and backing material, 
respectively, in watts per sq. cm per degree per 
cm, p and p’ the respective densities in g per cc, 
and C, and C,’ the respective specific heats in 
Joules per g degree. Then the respective thermal 
diffusivities K and K’ are constants given by 


K=k/pC, K’=k'/p'C,’. 


and 


The electrical analog of the backed bolometer 
strip is a short transmission line (Fig. 2) having 
distributed series resistance Ro and shunt capacity 
Cy per unit length, looking into a longer trans- 
mission line representing the backing and having 
distributed series resistance Ro’ and shunt ca- 
pacity Co’ per unit length. The front-surface 
dissipation of the bolometer strip corresponds to 
a conductance Go hung across the input end of 
the short line, while the backing-to-holder dissi- 
pation corresponds to a conductance G, termi- 
nating the long line. The following table of 
equivalents may be found useful in visualizing 
the correspondence of analogous thermal and 
electrical entities (Table I). 

The problem is now that of finding the input 
voltage amplitude Vo|2-0 as a function of fre- 
quency when a current 


o(1+€?**) 
is fed into the network. The a.c. solution of this 
network may be obtained by considering the 
problem as that of a short line of input impedance 


Z (taken just inside the input dissipation con- 
ductance) and terminated at x = Ax by the input 











TABLE I. 
Entity Thermal Electrical 
Input flow, a.c. part Ue?! ine?** 
Input amplitude To| 2=0 Vo|z=0 
Series resistance per unit length 
1 
bolometer z Ro 
r 1 
backing Pp Ro’ 
Shunt capacity per unit length 
bolometer pC. Co 
backing p’C,’ Co’ 
Front-surface dissipation o(Tz-0—Ts) GoV\|z-0 


Backing-to-holder dissipation  o/(Tz.1—T.) GriV|2or 
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impedance Z’ of the long line, which in turn is terminated at x=L by the dissipation G,. Thus, 
the solution is 


Vo! zm0 = 10 (Go+1 Z) ) 
where Z is given® by 
| 


_ ZosinhyAx+Z’ coshyAx 





4 


0 ’ 
Zy coshyAx+Z’ sinhyAx 
and where, in turn, Z’ is given by 


1 
Zo sinhy’(L— Ax) +— coshy’(L— Ax) | 


Gr 


-_ =Zy' 





1 
Zo coshy’(L— te sinhy’(L— a) | 
TL 


In these equations Zp and y are the characteristic impedance and propagation constant, respectively, 
of the short line representing the bolometer strip, and Zo’ and y’ are the corresponding constants for 
the long line representing the backing. For series-resistance shunt-capacity lines these quantities 
have the values: 
Zo = | Ry ‘jCw)} 
7 (jRoCw)! 
Zo - (Ro’ jCo'w)' 
= (jRo'Co’w)! 


) 
F (3) 


With the quantities (3) replaced into (2), and with the corresponding thermal equivalents sub- 
stituted into the resulting expressions, the amplitude of the temperature fluctuation at the front 
surface of the bolometer strip becomes 


To| <0 = Uo/(o+1/Zu), ) 


1 , w\? jor! 
-( ) sinh(~*) Ax+Z' th cosh() Ax 
jw K K 


where 


| x 
~. 





9 


5 jw } 
Ax+Z’ th sinh( =) Ax 
K 


1/K\'k 
Zu=-(—) 

k\ jw 1/K\3 jw 
-(—) cosh{ — 
k\ jw K 


mits —) sinn(2*) L- as) + cosh (2 “) (L—Ax) 
Z'un= ~) + . 
R’ K’ jw jw} 
<8 cosh( = "\ (L—- ax)-+5 sinh (= ) (L—Ax) | 
R jw K’ K’ 


_ *A derivation of this expression for the input impedance of a distributed-constant line terminated by another 
impedance can be found, for instance, in Slater’s Microwave Transmission, Chap. I 


(4) 


and where 
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er 





The subscript th indicates that the impedances 
designated are thermal impedances. 


Application to Air-Backed Bolometers 


The above equations can equally well be 
applied to air-backed bolometers. For an air- 





To | r=0 = 


backed element the only heat loss is by air 
dissipation at the front and rear surfaces. The 
thermal impedance of the backing in this case is 
replaced by the appropriate air-dissipation factor, 
so that Z's, in (4) becomes simply Z’,=1/c. 
With this substitution the first two equations of 
(4) combine to give 


Uo 





1/K\} joy! 1 jo\) 
(—) cosh(=) Ax+- sinh(=) Ax 
jw K o K 





jw tk 
mG 
K/ 1 


Two approximations to (5) can be made, each 
being valid over a particular range of parameters, 
which simplify the task of interpreting physically 
the results obtained from (5). For small argu- 
ments the sinh is nearly equal to the argument 
and the cosh is nearly equal to unity. Thus, for 


| (jw/K)'Ax| $0.3 ) 

Eq. (5) becomes 
Uo(oAx+k) , (5a) 
jk?wAx 


To | r=0— 








a(oAx+2k)+ 





For large arguments the sinh and cosh are 
nearly equal. Hence, if 


| (jw/K)'Ax| 31.8 
Eq. (5) becomes 
, (5b) 
Uo 


To| z-0™ . 
o+k(jw/K) 





Equations (5), (5a), and (5b) can be used to 
calculate theoretical response-frequency curves 
for various bolometers. By thermal response is 
meant the ratio T|,-0/Uo of the amplitude of 
the front-face temperature fluctuation to the 
amplitude of the incoming chopped radiation 
signal. Such a response-frequency curve is given 
by the solid curve of Fig. 3 for a typical bolometer 
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K\3} jw? 1 jor} 
-(—) sinh( =) Ax+- cosh(=) Ax 
RN jw K o K 





strip made of thermistor material* and having: 


Ax = 0.0018 cm 
k=0.01 watt/cm deg. 
p=5 g/cc 

C,=0.75 joule/g deg. 
o =0.005 watt/sq. cm deg. 


This curve is composed of the following three 
distinct regions, each being characteristic of a 
particular frequency range. 

I. A frequency-independent region for very 
low frequencies in which, physically, the tem- 
perature of the entire flake is always in equi- 
librium with the incident radiation. From Eq. 


WITH CONSTANTS AS GIVEN BELOW 








WITH @ MADE “| REGION I 
5 TIMES LARGER 


SOLID CURVE : 
4k 0.0018 CM 
@ = 0.005 WATT/EM"' DEG 


p= 375 JOULES /CM’ DEG 
k= 0.01 WATT /CM DEG 


RESPONSE, DEGREES PER WATT PER SQCM 
i] 


' 
w, 22M . ad 


Fic. 3. Response-frequency characteristic for 
air-backed bolometer. 
* Thermistor material is a mixture of certain semi- 
conductor oxides. It has a large negative temperature 
coefficient of resistance. 
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(5a), appropriate to this region, it is apparent 
that the temperature fluctuations of the strip 
are in phase with the fluctuations of the incoming 
radiation. 


Il. An intermediate region in which the re- 
sponse is inversely proportional to the frequency. 
Here the strip temperature is no longer in 
equilibrium with the incoming radiation, and the 
phase of the former lags that of the latter by a 
quarter cycle. As will be shown below, the 
instantaneous temperature of the whole strip is 
effectively equal to that of the front surface over 
the entire frequency region for which this rela- 
tionship holds. The thermal time constant of the 
strip is the reciprocal of the frequency at which 
the transition from region I to region II occurs. 

Ill. A high frequency region in which the 
response obeys a 1//w relationship, as given 
by (Sb). Here the phase of the front-face tem- 
perature lags that of the incoming radiation by 
an eighth cycle. 

The physical meaning of the transition from 
region II to region III can be understood from a 
consideration of the speed v and wave-length \ 
with which waves of temperature traverse the 
thickness of the strip in consequence of the 
periodic heating and cooling of the front face. 

The propagation constant y of the electrical 
line representing the bolometer, as given by 
Eq. (3) is: y=(jRoCw)'. The imaginary part 
of this quantity, called the phase constant 8, is: 
B=(RoCw/2)'. And the velocity v of electrical 
waves along this network is given by v=w/8 
=(2w/RoCo)'. The thermal analog of this 
velocity expression is Utherm=(2Kw)!, while the 
wave-length \ of the temperature waves is 

2rv 
\=— =22(2K/w)}.* (6) 

Ww 
Comparing (6) with the condition in (5b) for 
which the transition from region II to region II] 
becomes practically complete, it is evident that 
Aum ~4Ax. Thus, the transition corresponds, with 
increasing frequency, to the physical situation 
in which the temperature waves have become so 
short that their wave-length is comparable to 


* These expressions for wave-length and velocity of 
propagation of temperature waves in a_ semi-infinite 
medium were derived (otherwise) and used by Kelvin in 
his studies of the penetration of diurnal and annual waves 
of temperature into the earth. 
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the thickness of the strip, and the strip is no 
longer at the same instantaneous temperature 
throughout its entire thickness. Hence for chop- 
ping frequencies lying in region III various 
layers within the bolometer strip have tempera- 
tures different from 7'|\,~0, and the electrical 
output response of the bolometer bridge is no 
longer determined by To| 2-0, but rather by the 
average throughout the thickness of the strip 
of the local temperature fluctuation amplitude, 
with due regard for phase. It can be shown* 
that this electrical response continues to decrease 
like 1/w with increasing frequency, even though 
the front-face temperature fluctuation amplitude 
falls off like 1/+/w. It appears from these con- 
siderations that for chopping frequencies in 
regions I and II this solution could equally well 
have been obtained from a lumped-constant 
analogy as from a transmission line analogy. 

In Fig. 3 are plotted two dashed curves show- 
ing how the response-frequency characteristic is 
altered if a change is made in the thickness of 
the strip or in the value of the dissipation con- 
stant. The dashed curve characteristic of the 
5-times-greater dissipation constant becomes 
identical with the solid curve for #30. It can 
be seen from this illustration that if, by backing 
or otherwise, the dissipation from the element 
can be increased, an advantage will be gained in 
the operation of the bolometer. For, with this 
5-fold greater dissipation a \/5-fold greater d.c. 
bias can be placed across the bolometer bridge, 
and for w= 30 a «/5-fold increase in the electrical 
signal to the amplifier will be obtained for a 
given radiation signal. Furthermore, the time 
constant has been decreased to 1/1/5 of its 
former value. 

Unfortunately, there is no solid backing ma- 
terial whose dissipation is entirely conductive, 
and some thermal reactance due to the volume 
heat capacity of the material is inescapable. The 
effect of this reactance, as will be shown in the 
next section, is to modify the advantage gained 
by the use of a solid backing. 


Application to Solid-Backed Bolometers 


As an example of the kind of response-fre- 
quency curve to be expected for a solid-backed 


* The author is particularly indebted to W. H. Brattain 
for helpful discussion of this subject. 
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bolometer there is shown in the solid curve of 
Fig. 4 the characteristic of the same strip 
bolometer as was used in the calculations for 
Fig. 3, but now fastened intimately to a rocksalt 
backing having the following dimensions and 
properties. 


L=0.1 cm: area=1 sq. cm 

k’ =0.07 watt/cm degree 
p’C,’ =2 joules/cm* degree 

o’ =0.1 watt/cm? degree. 


In the same figure appear a number of dashed 
curves showing the effects produced by varying, 
one or two at a time, some of the parameters of 
the backing. The caption adjacent to each 
dashed curve explains which parameter has been 
changed, and by how much. All of these curves 
have been computed from (4), using, wherever 
appropriate, the type of simplifications indicated 
in (5a) and (5b). 

These curves exhibit features which are capable 
of simple physical explanation. The peculiar 
shape of the solid-line characteristic A is caused 
by the effect of the thermal reactance associated 
with the volume heat capacity of the backing. 
If this reactance could be eliminated, as for 
curve B which, to illustrate an extreme case, 
has been computed for p’C,’=0, the character- 
istic would resemble that of an air-backed unit 
with unachievably high surface dissipation. This 
circumstance is of importance in the practical 
selection of a backing material. It can be con- 
cluded generally that for backings having low 
p’C,’ the response can be expected to remain 
frequency-independent out to higher chopping 
frequencies than for backings having comparable 
k’ but higher p’C,’. 

Curve C, which has been computed for the 
extreme case k’=0, appears to coincide exactly 
with the solid curve of Fig. 3 except in the 
frequency-independent region, where the limiting 
response is twice that of the air-backed unit. 
With k’=0, the rear surface of the bolometer 
strip is thermally insulated, so that it now has 
only half the total dissipation of the air-backed 
unit. Of course, such a low-conductivity backing 
would be of no practical value in dissipating the 
bias heat of the element. 

In curve D, to illustrate still another extreme 
case, the backing has been given an infinite 
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Fic. 4. Response-frequency characteristic for 
rocksalt-backed bolometer. 


thermal conductivity. This is equivalent to 
placing the entire thermal capacity of the back- 
ing directly in contact with the rear surface of the 
bolometer strip. This unbuffered reactance brings 
about the rapid decrease in response in the 
range from w=10 to w=100 and results in the 
nearly frequency-independent response in the 
range from w=500 to w=10,000. In the latter 
range the frequency is so high that the heat 
capacity of the backing presents practically a 
‘‘short-circuiting’’ reactance to the rear surface 
of the bolometer strip. Thus, while the tempera- 
ture at the front surface of the strip is being 
driven by the fluctuation of the incoming radia- 
tion, the temperature of the rear surface is being 
held nearly constant by this reactance. And, 
while the whole strip is at the same phase 
throughout, it is not at the same temperature 
throughout. This state of affairs does not change 
much with increasing frequency until so high a 
frequency is reached that the wave-length of the 
thermal waves in the strip becomes shorter than 
the thickness of the strip. Then, of course, the 
front-surface temperature of the strip becomes 
independent of the nature and properties of the 
backing. 

It is to be noted that in the solid curve A 
there appears some tendency toward a levelling- 
off in the neighborhood of w~ 100. The reason is 
as explained in the last paragraph, except that 
here the full short-circuiting effect of the backing 
reactance is buffered by the distributed thermal 
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resistance of the backing, and the levelling-off is 
not as pronounced. 

Curve E shows the modification in the re- 
sponse-frequency curve brought about by making 
the backing five times thicker than for the 
solid curve A. Increasing the thickness of the 
backing increases the heat capacity, and the 
parallel reactance of this extra capacity makes 
itself felt at a lower chopping frequency than 
for curve A. Curve E tends toward a higher low 
frequency limiting response than curve A_ be- 
cause of the decrease in dissipation, brought 
about by the increased thermal resistance of the 
thicker backing. 

Finally, curve F shows the effect of making 
the backing-to-holder dissipation ten times larger 
than for curve A. This increased dissipation 
lowers the low frequency limiting response and 
makes it necessary to go to higher chopping 
frequencies before coming to the fall-off due to 
the heat capacity of the backing. 

It is often stated that an ideal backing material 
is one having a large value of the diffusivity K’. 
In making such a statement a distinction must 
be made whether the diffusivity is desired to be 
large by virtue of a large thermal conductivity 
k’ or by virtue of a small volume heat capacity 
p’C,’. A comparison of curves B and D shows 
why this distinction is necessary. In a bolometer 
assembly of such geometry and thermal proper- 
ties as are considered here the highest permissible 
d.c. bias across the bolometer strip is deter- 
mined chiefly by the backing-to-holder dissipa- 
tion rather than by the thermal conductivity of 
the backing. As long as the latter is large enough 
(as it is for rocksalt) so that the backing offers a 
comparatively small thermal resistance between 
the element and the holder, there is little virtue 
-in having k’ still higher, as, for example, for 
quartz. However, any decrease in the value of 
p’C,’ improves the response in the frequency 
region from w~1 to w~1000 and also decreases 
the time constant, both without affecting the 
permissible d.c. bias. 


DISCUSSION 


In the interest of simplifying the development 


of this presentation a number of short cuts and . 


omissions have been employed which would have 
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to be taken into account in the application of 
the method here outlined to any physical problem 
where accuracy of theoretical prediction is de- 
sired. There will be pointed out, in the next few 
paragraphs, the nature of some of these simplifi- 
cations, their consequences, and, so far as pos- 
sible, their specific remedies. 

All the equations appearing in this paper were 
derived on the assumption of a single sinusoidal 
amplitude fluctuation of input radiant power. If 
the fluctuation is obtained with a sector disk 
producing a nearly square-wave input of radiant 
power, a more satisfactory solution would be 
obtained by solving the appropriate equation for 
the fundamental chopping frequency and two or 
three of its harmonics as well, assigning to each 
harmonic solution a weight corresponding to the 
Fourier coefficient of the particular harmonic in 
the square-wave expansion, and summing over 
the expansion. 

While backed bolometers can be made for 
which the thermal contact between the element 
and its backing is intimate, it is usually necessary 
in practice to fasten the two together with a 
thin layer of cement whose thermal properties 
may contribute to the behavior of the whole. 
The backed bolometers discussed in connection 
with the curves of Fig. 4 have been idealized to 
the extent of neglecting this layer and its effects. 
A paper is now in preparation by W. van 
Roosbroeck in which some of these details are 
treated. Van Roosbroeck’s treatment will also 
include a consideration of cases in which the 
incoming radiation, instead of being absorbed, 
as here assumed, at the front surface of the 
bolometer strip, is absorbed throughout the 
thickness of the strip. 

In usual practice a backed bolometer is made 
of a narrow strip of sensitive material cemented 
to a backing slab much larger in area than the 
strip. Hence the physical problem is not one 
dimensional except in a small volume of the 
backing immediately adjacent to the strip. At 
larger distances away from the strip in the 
backing the heat flows along nearly radial lines 
and the isothermals are nearly semi-cylinders 
coaxial with the center line of the strip. The 
analogous electrical problem is one in which the 
series conductance and parallel capacitance of 
the transmission line representing the backing 
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increase nearly linearly with distance along the 
line. Except for special cases exact solutions for 
varying-parameter lines are not easily obtainable. 
However, approximate electrical analogs can be 
set up which are capable of solution. Such 
approximations, over certain frequency ranges 
are assisted by the fact that temperature waves 
in a semi-infinite solid are attenuated in ampli- 
tude by a factor e°* per wave-length. Hence, 
for higher frequencies the temperature waves, 
which become shorter in wave-length with in- 
creasing frequency, do not penetrate far enough 
into the backing to get into the region where 
the fanning-out of the heat flow lines occurs, 
and the constant-parameter solutions already 
derived can be used as good approximations. 

In deriving the equation for the solid-backed 





element the assumption was made that the 
holder presented a non-reactive dissipation to 
the rear surface of the backing slab. This assump- 
tion could lead to error only in frequency region I 
where the chopping rate is so low that the tem- 
perature waves penetrate clear through the 
backing. For higher frequencies the exact nature 
of the holder’s impedance would not be of any 
consequence. 
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The Speed of Propagation of Brittle Cracks in Steel 


G. Hupson* AND M. GREENFIELD** 
David Taylor Model Basin, USN, Washington, D. C. 
(Received September 30, 1946) 


Brittle failure was produced in notched steel tensile specimens. The speed of propagation of 
the brittle crack was measured and found to be approximately 40,400 in./sec. 


HERE have been recent reports of the 

brittle failure of steel plate used in a variety 

of engineering structures. This rupture occurred 

when the steel was subjected to certain severe 

conditions. In the main these brittle failures may 
be characterized by the following: 


(1) A crack is propagated with extreme rapidity through 
the steel plate, 

(2) The surface of rupture created by this crack is 
orthogonal to the plane of the plate, 

(3) This rupture surface shows a characteristic “herring- 
bone” pattern, 

(4) There exist small permanent strains, often of less 
than 2 percent, in the neighborhood of the fractured 
surface. It was thought that a study of the speed of 
propagation of brittle cracks in steel might contribute to an 
understanding of the mechanisms of these phenomena. 


A preliminary note on the beginning of this 
study has been given in an earlier paper.' This 
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work was subsequently extended with a more 
refined technique, to include more variation in 
the path length over which the average velocity 
of the crack was measured, and a variation in the 
distance of this path from the point of initiation 
of the crack. In the present paper, in addition to 
reporting the later results which agree with the 
previous ones within the limit of experimental 
error, more details are given of the technique 
used. 

It has been found possible to produce brittle 
failure in a steel tensile specimen at room tem- 
peratures in a. controlled experiment of the 
following type. The specimens consisted of rect- 
angular plates of a plain carbon structural steel 
(also known as medium steel) with a carbon 
content not exceeding 0.31. This steel had a 
minimum average ultimate strength of 60,000 
lb. /in.2, and a minimum average yield strength of 
34,000 Ib./in.2. This steel had not been heat 
treated. The specimens were approximately 2 
feet long, 53 inches wide, and 1 inch thick. To fix 
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Fic. 1. Schematic diagram showing the positions of 
one-mil Nichrome break wires 1, 2, and 3 on the tensile 
specimen. 


the point of initiation of the failure, a notch 
consisting of a milling saw cut of the order of 35 
inch in width and } inch deep was machined. into 
one edge of the specimen midway between its 
ends. See Fig. 1. 

The specimen was placed in the grips of a 
tensile machinef and pulled to failure. A brittle 
crack initiated at the milling cut, as indicated by 
Fig. 1, was propagated across the specimen, at 
right angles to the edge. 

An attempt was made to photograph the crack 
as it spread through the plate. The camera used 
was a Fastax Western Electric, operating at 
about 3000 frames per second. Even this speed 
was insufficient to catch the phenomenon, as the 
crack invariably spread across the specimen be- 
tween two successive frames. It was possible, 
however, to infer by this method a lower limit of 
the order of 5000 in./sec. for the velocity. 

The following method using electronic recording 
was finally adopted for determining the speed of 
propagation of the crack in these specimens. The 
specimens were first prepared by heating them to 
a temperature sufficient to melt but not boil a 








t These specimens were actually tested in two different 
types of machines; one is a Baldwin Southwark hydraulic 
type with a 600,000 Ib. capacity located at the David 
Taylor Model Basin, USN, Washington, D. C., and the 
other, an Olsen screw type with a 600,000 Ib. capacity, 
located at the National Bureau of Standards, Washing- 
ton, D. C. 
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Fic. 2. A typical record obtained from the breaking of a 
tensile specimen like that shown in Fig. 1. The numbers 
1, 2, and 3 designate the positions of the electron beam at 
the instants of breaking of wires 1, 2, and 3, respectively. 


layer of medium de Khotinsky cement. This 
cement was spread smoothly on the surface of the 
metal in the region through which the crack was 
to be propagated. While the de Khotinsky was 
still fluid, several lengths of one-mil Nichrome 
wire were embedded in the cement, as shown 
schematically in Fig. 1. These wires were arranged 
to be approximately parallel, about one inch 
apart, and perpendicular to the expected path of 
the crack. The de Khotinsky thus served not only 
to cement the wire firmly to the tensile specimen, 


TABLE I. Speed of propagation of brittle crack. 


Testing Specimen Wire Distance Average velocity 
machine number numbers inches inches/second 
TMB E-1 2-3 2.23 42,700 
TMB E-3 2-3 4.41 37,900 
NBS F-3 2-3 1.01 36,500 
F-3 3-4 0.98 43,000 
F-3 4-5 0.96 39,500 
F-3 5-6 0.96 41,200 
NBS G-2 3-4 0.96 42,400 
G-2 4-6 1.99 41,600 
NBS G-3 2-4 1.97 44,200 
G-3 4-5 0.96 37,400 
G-3 5-6 0.98 43,700 
NBS G-5 3-4 1.00 41,200 
C-5 4-5 0.99 33,000 
C-5 5-6 0.98 41,000 


Average of above velocities = 40,400 +2500 in./sec. 
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but also insulated the wire from the specimen. 
The ends of the wires which protruded from the 
cement were insulated from the specimen by 
strips of paper stuck to the specimen partly by 
the de Khotinsky cement and partly by Scotch 
Tape. After cooling the specimen, the ends of the 
wires were then soldered to binding posts on 
Bakelite strips. These were held by rubber bands 
to the specimen about two inches on either side 
of the path of the crack. 

When a crack is propagated from the base of 
the saw cut, it breaks the fine wires in succession. 
The breaking of the first wire is used to start a 
single horizontal sweep of the electron beam on 
the screen of a cathode-ray oscillograph. Each of 
the succeeding wires, on breaking, introduces a 
voltage step on the vertical deflection plates by 
means of a simple voltage divider arrangement. 
The resulting staircase-like record may be photo- 
graphed, e.g., see Fig. 2. The rounding off of the 
successive steps is caused by the impedance of the 
circuit used. A calibration wave for timing is 
provided by a standard oscillator. See Fig. 2. 

The average velocity of the crack between any 
two successive wires is calculated by dividing the 





measured distance between them by the corre- 
sponding time interval measured from the 
oscillograph record. 

Table I gives the final results. Column 1 states 
where the specimen was tested, i.e., either at the 
Taylor Model Basin or at the National Bureau of 
Standards. Column 2 is the specimen number. 
Columns 3 and 4 give, respectively, the numbers 
of the 1-mil Nichrome wires between which the 
speed was measured, and the distance between 
those wires. For example, Specimen F-3 had six 
(6) 1-mil Nichrome wires cemented to it. The 
breaking of the first wire started the single sweep 
across the c.r.o. screen. Then there were signal 
wires 2 through 6 spaced as indicated in column 4. 
Column 5 gives the average velocity for the 
corresponding distance listed in column 4. 

The average of all the velocities listed in 
column 5 is 40,400 in./sec. with a mean absolute 
deviation of 2500 in./sec. 

It is intended to continue this work to learn 
whether there is a significant variation of the 
velocity of the brittle crack with the thickness 
of the steel plate, its temperature, and its 
composition. 





Use of the Broadcast Band in Geologic Mapping* 


LARKIN KERWIN** 
Massachusetts Institute of Technology, Cambridge, Massachusetts 


(Received December 4, 1946) 


A review of the literature indicates that a method of geologic mapping by means of observing 
the effect of geologic anomalies on electromagnetic field intensity should be practical. Field 
equipment is designed and assembled for that purpose and proves accurate and simple to use. 
Preliminary investigations with it successfully locate several anomalies, 


INTRODUCTION 


HE practicability of using high frequency 
electromagnetic waves in sub-surface geo- 
logic investigation is generally viewed as dubious. 
Joyce's theoretical conclusions from Maxwell's 


* This paper is based on a portion of a thesis entitled 
“The Development of Field Equipment Suitable for 
Measuring the Significant Properties of the Earth and 
Their Effect on Electromagnetic Wave Propagation,”’ 
submitted for credits towards the degree of Master of 
Science at M.I.T. Financial support for the work was 
provided through a grant from The Geological Society of 
America, which is hereby gratefully acknowledged. 

** Now with the Department of Physics, Laval Uni- 
versity, Quebec. 
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work! promise great attenuation under the as- 
sumed conditions. Eve,? Terman,’ and Ollendorf 
support these considerations. However King 
predicted less attenuation, while Fritsch‘ declares 
that theoretical assumptions are never justified 
in practice and suggests that practical procedures 
are quite possible. 

Relatively few experimental investigations 


1 J. W. Joyce, Bureau of Mines Tech. Pap. No. 497 (1931). 

2A. S. Eve and D. A. Keys, Applied Geophysics (Cam- 
bridge University Press, New York, 1929). 

*F. E. Terman, Radio Engineers’ Handbook (McGraw- 
Hill Book Company, Inc., New York, 1943), p. 698. 

4\V. Fritsch, Beit. Ang. Geophysik [4] 5, pp. 375-395. 
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have been undertaken and these are inconclusive. 
Potter and Friis® determined that ordinary 
topographic variations did not have great effects 
on high frequency field intensities. On the other 
hand, Cloos,* Howell,’ and Fritsch* found that 
field intensities did vary with the sub-surface 
geology. Joyce, Eve,’ Allen, and Silverman and 
Sheffet'® worked qualitatively in mines or tunnels 
and found in general that broadcast frequencies 
often penetrate the ground for considerable 
distances, but that higher frequencies were 
quickly absorbed. Joyce also did some quantita- 
tive work at 500 cycles, but it is highly doubtful 
if this frequency justified the use of electro- 
magnetic wave theory in considering the results. 
The effect was undoubtedly a simple inductive 
one in view of the tremendous wave-length. 
Fritsch": lists many interesting results of suc- 
cessful underground propagation. 

These results and others indicated that the 
effect of geology on the field intensity of electro- 
magnetic waves might provide a useful method 
of geologic mapping. The Geologic Society of 
America provided a grant for the purpose of 
developing equipment, and preliminary results 
obtained with it last summer indicate that the 
band from 200 ke to 1000 ke will prove useful 
in future work. 

It was decided that the investigation of the 
field patterns of various frequencies over the 
same area would provide useful information. 
Geologic anomalies would cause irregularities in 
the pattern which would vary with frequency 
according to the penetrating power of the various 
wave-lengths. Equipment was purchased and 
adapted or constructed to-produce, monitor, and 
measure suitable electromagnetic fields. 


EQUIPMENT 


The field intensity recorder was assembled 
inside a wooden-bodied station wagon supplied 


6 R. K, Potter and H. T. Friis, Proc. I.R.E. [4] 20 (April, 
1932). 

®E. Cloos, Am. J. Sci. [166] 28 (October, 1934). 

7 Howell, Geophys. [2] 8 (April, 1943). 

sV. Fritsch, Beit. Ang. Geophysik [3] 5, pp. 315-364. 

® Eve, A.I.M. E. Tech. Pub. No. 316 

10 p- Silverman and D. Sheffet, ote [4] 7 (October, 
1942 

“ V, Fritsch, Beit. Ang. Geophysik [1] 6, pp. 100-119; 
(36, pp. 277-306; [4] 6, pp. 407-412; [417, pp. 449-461. 

2 VY. Fritsch, Neues Jahrb. f. Geol. B., Vol. 84, Series 1, 
pp. 90-116. 
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by the Geotechnical Corporation. Its basic com- 
ponent was a modified Model DAE radio direc- 
tion finder. The antenna consisted of a 14-inch 
shielded loop, mounted in a stationary cylinder 
provided with azimuth scale and commutator 
rings connected to a shielded lead. This loop was 
mounted on a spring tripod which was bolted to 
the floor of the station wagon. This mount was 
firm, but protected the loop from the shocks 
produced by the car’s motion. 

The associated receiver was mounted on a 
small oak table, directly in front of the operator's 
seat. This unit was of the superheterodyne type, 
covering the frequency range 240 kc to 2000 kc 
in three bands. It was provided with volume 
control, antenna tuning, and beat frequency 
oscillator, which made the receiving of unmodu- 
lated signals considerably easier. The case was 
provided with shock mounts screwed to the 
reinforced table top. In this unit the field pro- 
duced a proportionate I.F. voltage, and that of 
the 2nd I.F. stage was selected as the most 
satisfactory for recording purposes. It was 
tapped to a jack on the panel, the signal from 
which was then proportional to the field strength, 
but could be controlled by the volume control 
circuit. 

This voltage was measured by a vacuum-tube 
voltmeter of the ““Vomax”’ type, installed im- 
mediately above the receiver on shock mounts 
attached to the latter’s case. Sponge rubber pads 
were glued to the roof of the wagon and the 
voltmeter as added protection. For continuous 
recording, an Esterline-Angus recording milli- 
ammeter was used, having speeds up to six 
inches per minute. It was connected into the 
voltmeter circuit via a jack mounted on the 
meter face. -By using this system throughout, 
any component could be removed for alteration 
or replacement without necessitating a rewiring 
job each time. For transportation and moderate 
speeds of recording the meter was mounted in 
sponge rubber padded clamps on a shelf under 
the table. Here it was easily available for inking, 
marking the chart, etc. The needle was. weighed 
somewhat more than usual, and it was found 
that it would not jump from the paper if the 
car was driven at low speeds. For continuous 
operation at high speeds, a low frequency spring 
suspension should be used. 
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All of these units were operated with 110- 
volt, 60-cycle power. In the mobile unit this was 
provided by a Janette rotary converter operated 
by storage batteries. This produced a steady 
voltage at about 200 watts and operated silently. 
When well charged, the batteries provided 24 
volts to operate the equipment for about 8 
hours, more than a normal day’s field work. 
A battery charger was mounted underneath the 
operator’s seat and connected so that after the 
day’s operations a long power cable could be 
connected to any convenient 110-volt outlet and 
the batteries charged overnight. This system 
proved completely satisfactory. A constant-volt- 
age transformer proved useful for last-minute 
work at the day’s end, in case the batteries began 
weakening. A general view of the equipment in 
the station wagon is given in Fig. 1. 

It was intended that frequency ranges from 
200 ke and up should be investigated. For pre- 
liminary work the broadcast band seemed suit- 
able, and first field runs were made by recording 
the fields of various commercial transmitters as 
described later. For easily controllable work and 
for greater frequency range, a compact 75-watt 
transmitter was obtained. Known as the BC- 
375-E in the U. S. Army Signal Corps for whom 
it was designed, this transmitter provided con- 
tinuously variable frequencies from 200 kc to 
500 ke and from 1500 ke to 12,500 kc by means of 
six plug-in tuning units. This transmitter was 
designed to operate from 24-volt d.c. supply; 
plate voltages, etc., being supplied by a dyna- 
motor. This supply provided approximately 750 
watts. For portable operation battery operation 
was indicated, but was impractical because of 
the heavy drain of about 30 amperes. Accord- 
ingly, after some search a small gasoline driven 
generator was obtained which provided the 
required power, and delivered it through a 
number of ‘“‘floating’’ storage batteries. This 
provided a smooth d.c. supply of about 28 volts 
for which the transmitter was easily modified. 

During the process of taking field intensity 
measurements Over various ground conditions, 
it was possible that the electromagnetic field 
would undergo certain variations independent of 
the geology. Fluctuations at the transmitter, 
atmospheric changes, the time of day at which 
readings were taken, etc. all contributed to 
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Fic. 1. Field intensity recording equipment mounted 
in station wagon. 


fluctuation. It was therefore necessary to provide 
a stationary monitoring field intensity meter 
which would keep a record of these changes. The 
fact that relative field strengths and not absolute 
measurements were required simplified matters 
somewhat. A Hallicrafters Model SX-28 all-band 
receiver was obtained and provided with a 
simple fixed antenna. A vacuum-tube voltmeter 
was connected into the a.v.c. circuit of the re-. 
ceiver. Thus by the adjustment of the r-f amplifier 
amplitude and the use of a suitable voltmeter 
scale, a convenient deflection was obtained. This 
followed accurately the variations in field in- 
tensity. In order to obtain a continuous record 
of these changes, an Esterline Angus recording 
milliammeter was placed in the voltmeter circuit. 
This meter was adjusted to record at the rate of 
1 foot per hour, which provided a satisfactory 
record of field strength changes without noting 
any minor fluctuations. 

All of this equipment operated from 110-volt, 
60-cycle current which was usually available 
near the scene of operations, or could be provided 
by a small motor generator. 

In the course of the field work which was 
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Fic. 2. Field intensity recordings (taken from recording meter chart). 


extended over several weeks, the station wagon testing, and after a few adjustments following 
was driven over all kinds of roads, and a com-_ the first couple of days’ experience, the appa- 
plete gamut of speeds. The equipment and ratus proved entirely field-worthy. The only 
mountings were thus subjected toa most thorough precaution taken was to remove the recording 
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milliammeter needle if it was intended to proceed 
over rough ground at:high speed. At low speeds 
it recorded satisfactorily. On this score, it is 
believed that future equipment patterned after 
this will provide a simple and rugged tool for 
the investigator. 


FIELD PROCEDURE 


For preliminary work, a dike situated in 
Newton, Massachusetts was selected as a suitable 
testing ground. This dike had been roughly 
located during a survey by Mr. Blumberg of 
Harvard, who kindly lent some of his field maps 
to this project. An extension of basic melophyre 
was found protruding through the roxbury con- 
glomerate which was characteristic of the area. 
The contrast in the physical properties of these 
rocks is considerable, making them quite suit- 
able for the purpose in hand. 

The district was suburban, with considerable 
vegetation and the houses widely separated. 
Some time was spent in examining the geology 
in order to locate the dike factually if possible. 
In Fig. 3a the dike location based on surface 
geology is indicated. 

Field intensity measurements were then made 
under a variety of conditions. First, detailed runs 
were made by taking field intensity records of 
about five minutes duration at points about 
twenty-five feet apart along the road. In each 
case the location of the loop and its orientation 
were established by taping and bearings taken 
with a Brunton-type pocket transit. The latter 
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Fic. 3b. Field intensity contours. 
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readings were taken from two directions and at 
some distance from the station wagon, elimi- 
nating the latter’s effect on the compass. The 
bearings so taken were all consistent with those 
established from surveyed maps. 

The duration of the five-minute reading was 
important in order to achieve an averaging effect. 
An example of the readings obtained over several 
stations as taken from the field chart is given in 
Fig. 2a. One of the disadvantages of using the 
broadcast band with the present equipment here 
became apparent. The modulating voltage had 
some effect on the instantaneous readings, al- 
though not the average. Therefore readings 
extending over some period of time were neces- 
sary to obtain accurate values. The type of 
program being broadcast also has an effect, as 
the modulation of musical programs is much 
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Fic. 3c. Resistivity contours. 
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smoother than those of the “soap-opera’’ type, 
and gave records which were much easier to 
interpret. An example is shown in Fig. 2b. 
These difficulties are avoided when using the 
unmodulated wave of the local transmitter. 

The procedure was repeated several times, for 
several frequencies. Records were also made in 
the sunshine, rain, under dull clouds, and at 
night. 

In addition to this, several types of continuous 
recording were investigated. By keeping the 
loop in a constant position of initial maximum 
strength and driving the station wagon slowly 
over the area, a figure such as that shown in 
Fig. 2c was recorded. This shows the prominent 
features which are marked. Other runs were 
made by keeping the loop manually oriented as 
much as possible in the direction of maximum 
intensity. The result is shown in Fig. 2d, where 
it is seen to bear close resemblance to the former 
record. Note that runs were made in varying 
directions and must be properly interpreted as 
to orientation. Continuously rotating the loop 
so as to produce an average effect gave a record 
such as shown in Fig. 2e, which follows the 
same general pattern as before. 

A series of resistivity measurements was also 
made in the area in order to determine whether 
any relation existed between resistivity and field 
intensity. These measurements were made by 
inserting four metal stakes into the ground about 
twenty-five feet apart along the line of measure- 
ment. 150 volts, 60 cycle, as supplied by the 
small motor generator in the station wagon were 
applied between the two outer stakes. By meas- 
urements taken of the resistance between them 
the current put into the ground could be calcu- 
lated. The voltage drop between the two inner 
stakes was then measured with a vacuum-tube 
voltmeter, and from these readings the resistivity 
was calculated." By advancing the rear stake 
and changing terminals, successive readings were 
quickly taken. In this manner the resistivity was 
measured along the same direction as that in 
which the field intensity was recorded. 

From the accumulated data, field intensity 
contours were drawn on the field map, and are 
shown in Fig. 3b. We see that the contours 


4% C. A. Heiland, Geophysical Exploration (Prentice-Hall, 
Inc., New York, 1940), p. 709. 
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follow closely the direction of the dike as deter- 
mined by the surface geology. Resistivity con- 
tours were also determined and are shown in 
Fig. 3c, where they assume the same general 
pattern as the other figures. The resistivity of 
the dike is substantially lower than the con- 
glomerate, and as would be expected, the field 
intensity above the dike was also lower. 

Azimuths of maximum field intensity are given 
in Fig. 3d, but they do not appear to convey 
any conclusive information. Over the dike, on 
Oak Hill Road, they, seem to orient more 
towards the North than at other places. 

Certain characteristics of the graphs drawn 
from the accumulated data are of some interest. 
It was observed that data taken at night showed 
considerably less variation than that taken by 
day. Measurements of low frequencies showed 
greater disturbance over the dike than did those 
of higher frequencies. Data taken during rainy 
weather showed the fields to be in general 
stronger and less affected by geology. Low in- 
tensities over the dike pointed to its higher 
conductivity, as was verified by resistivity meas- 
urements. The low frequency irregularities were 
more widespread than the others, indicating that 
the ground anomaly was more than a surface 
effect. To some extent here, and more so in 
urban areas later investigated, the effects of 
manholes, steel rails, etc. were quite obvious and 
easily detected. 

Besides the detailed investigations in this 
area, other districts were investigated, using both 
stationary and rotating loop techniques. In the 
more urban area of Hammond Street, Brookline, 
a similar dibase dike protruding through the 
conglomerate was located by this method, and 
its presence verified by an examination of the 
surface geology. The boundaries of the dike were 
fairly sharply recorded as the field weakened 
over it. In Medford there is a noted extrusion of 
Medford dibase through Lynn volcanic complex. 
In this instance the resistivity of the dike is 
higher than the surrounding rocks, and the 
field intensity increased over it. 

It was also observed that topography appeared 
to have relatively little effect. In driving over 
stretches of country in Arlington where the 
geology was uniform, the field remained essen- 
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tially constant except for manhole effects, etc., 
whether the receiver was passing over rolling 


-hills, proceeding through cuttings, or on level 


ground. A notable exception was observed on 
the top of Pine Hill, a prominence beside the 
Medford Dike and very close to the transmitting 
antenna of the broadcast station being recorded. 
At the top of the hill the signal was extremely 
strong, but it weakened considerably as the bulk 
of the mountain was interposed sharply between 
the receiver and transmitter. 

For most of these field measurements the field 
itself remained very constant, as shown by the 
monitoring station. Variations were observed 
between daytime and nighttime records, and 
with weather; but over considerable periods (3-4 
hours) if the weather and daylight conditions 
were constant the field remained quite steady. 

In general it is considered that the results 
obtained thus far show sufficient promise to 
warrant a continuation of the research. The 
equipment has proven sufficiently rugged, sensi- 
tive and simple to operate over long periods of 
field measurements. The cost of operation in- 
volves only the depreciation of equipment, and 
station wagon expenses. The method of field 
intensity recording has been used to locate a 
geologic anomaly as determined by other inde- 
pendent means, and has generally located other 
geologic features. Accordingly, further work on 
the subject is being done. 


MODIFICATIONS 


In light of past experience, several modifica- 
tions are recommended for the equipment, and 
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will be incorporated into it shortly. The portable 
transmitter will provide a much more satis- 
factory field from the point of view of modulation 
as well as that of situation -and pattern. For 
continuous recording, the milliammeter is being 
equipped with speedometer drive. This will 
permit operation by one person rather than the 
operator and driver previously necessary. The 
loop antenna is being supplemented by a vertical 
antenna, the two being interchangeable by means 
of a switch. This will eliminate the necessity of 
rotating the loop when azimuths are not being 
taken. A switching system in the charging 
circuit will make that operation somewhat more 
convenient to perform. The monitoring unit is 
being made self-sufficient by powering it with 
another rotary converter operating from 12 volts 
d.c. The SX-28 Hallicrafters receiver is being 
replaced by a smaller, more portable one. 

Considerable attention is being devoted to 
frequencies below the broadcast band, as field 
experience here supports the general theory that 
low frequencies are more suitable. Experiments 
will also be performed on bands above the 
broadcast, up to about 12,000 ke. 
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A Note on a Paper by Faust and Beck 


WILLIAM M. STONE 
Iowa State College, Ames, Iowa 
(Received October 31, 1946) 


An infinite sum transformation is defined and applied to a system of linear difference 
equations discussed by Faust and Beck in their paper on single tuned amplifiers. Some trans- 
forms of the more common functions are given and points of superiority of the transform method 
over the classical methods of solution of difference equations are emphasized. 


N their paper, “Oscillation Conditions in 

Single Tuned Amplifiers,” Faust and Beck 
solved a simultaneous set of two first-order 
difference equations by classical methods. The 
purpose of this paper is to introduce the General- 
ized Laplace Transformation of Samuelson? and 
to illustrate its application to the same problem 
in circuit analysis. The advantages are that (1) 
a difference system is reduced to an algebraic 
system, (2) there are no extraneous constants 
introduced, and (3) no assumptions need be 
made as to the nature of the solution. 

The generalized Laplace transformation is 
defined as 


L{F(k)}e=f(s)=> F(n)s-"", (1) 


where F(k) must be of exponential order, 
F(k)< Ma", k>K, and defined for all positive 
integral values of k. The series is then uniformly 
convergent for all values of s>a. The subscript 
refers to the shifting operator, E™F(k) = F(k+m), 
and serves to distinguish the generalized Laplace 
transformation from the ordinary case. Some 
transforms may be found by elementary sum- 
mation but a more powerful method is readily 
available. By definition of the E operator 


LiE-F(k) \2=>, F(n+1)s—"—, 
0 


=s{ F(0)s7+ F(1)s?+---]—F(0), 
=sf(s)—F(0), (2) 
which may be generalized at once by mathe- 


! Faust and Beck, J. App. Phys. 17, 749 (1946). 
2 Samuelson, Bull. Am. Math. Soc., Abstract 52-3-86 
(1946). 
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matical induction to the form 


m—1 
L{E™- F(k)\p=s™f(s) — © s™--iF(j). (3) 


j7=0 


Following is a list of some of the more useful 
transform pairs: 


F(k) 
1 1 


f(s) 
(s—1) 


S5—da COSW 





a* cosw 
s*—2as cosw+a* 


a sinw 





a* sinwk — 
s*—2as cosw+a’ 


s—as coshw 





a* coshwk 
s?—2as coshw+a? 


a sinhw 





a* sinhwk 
s?—2as coshw+a? 


a*k™ a"n!/(s—a)"*}. 


Here k\ means the Newton factorial 
nomial, k(k—1)(k—2)---(k—n-+1). 

Using capital letters for the object function, 
small letters for the transforms, Faust and Beck's 
equations for voltage and current may be written 
as 


poly- 


E(k+1)—E(k)—ZiI(k)=0, (4) 
E(k+1)+8mZ2E(k) —Zol(k+1) 
+ZoI(k)=0. (5) 


By means of (2) these may be transformed at 
once to an algebraic system: 


(s—1)e(s) —Z1(s) = E(0), (6) 
(s+ 2mZ2)e(s) —Z2(s—1)i(s) =E(0)—1(0)Z>2. (7) 
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The solution of this system is found to be 


Zi Zi 
E00) - (: +=) |+10)2.-B0)— 
































(5) 2Z> 2Z2 (8) 
eas) =—— ’ 
Zi 
st—25(14+-—-) 41-02 
2Z2 
10) (: EO) ) 2 
+ load ean Zm 
s) I(0)Z2 (9) 
1(s)= : 9 
Z; 
st—25( 1+ )+1-enZs 
Substitute 
Zi Zi Zi . 4 
i+ |-+(—) +2nZs| 
2Z2 Z2 2Z2 
a=(1—gmZ,)', coshw= , sinhw= : 
a a 
=A a, = B/a, 
so that 
Z\ 
E(0)(s—A) ah aaa ia 
e(s) =— messin, (10) 
s*—2As+a? 
E(0) 
I(0)Z2 
i(s) = (11) 
s*—2As+a? 


Making the denominators conform with the list of transforms the object functions may be written as 














Zi 
E(k) =a*| E(0) coshwk-+——— * sinhwk |, (12) 
Z, E(0) 

I@)-—+ (1+gmZ2) 

2Z2 2 
I(k) =a*!? hoe coshwk + : sinhwk |. (13) 

The boundary conditions are 

I(n) =0, (14) 


where J(m) is the current in the last mesh, E(0) is the voltage in the input stage, Ei» is the voltage 
of the input generator whose internal impedance is Z,. Substitution of both conditions in (13) 
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and algebraic reduction leads to 


+1 
B coshwn +—- sinhwn 


+2 


E(0) =E;,— " (16) 
Zi. 1 
B coshwn +|= +2,(—+2.)| sinhw 
2Z> / 


<2 


1 
(— +2) sinhwn 
Z2 








I(0) = —E ix ' (17) 
Zi 1 
B coshaen-+| = +2,(—+2)| sinhwn , 


In practice the gain off all m stages is of most importance. If Z, is set equal to Z2 as is customary 
the gain of the m stages may be reduced to 


E(n) B(1—gmZ,)"? 
—= ' (18) 
E;, Bcoshwn+(A+gmZ2) sinhwn 





The well-known expansion formulas 





n 
coshwn = >-; ( ) (coshw)"~** (sinhw)?*, (19) 
. n , . 
sinhwon = Eo ) (coshes)=-# \(sinhw)?**!, (20) 
2k+1 
where i) represents the binomial coefficients as k takes on the values 0, 1, 2, ..., lead to the 
form of (18) 
E(n) (1—gmZ1)" 
= , (21) 
Ein n+1 n 
= )a n—2k Bk 9 Z, = jan -1B2k 
2k+1 2k+1 


which is equivalent to the general formula for the gain given by Faust and Beck. 
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Internal Friction of Zinc Single Crystals* 


Irvin H. Swirtt** AND JouN E. RicHARDsON*** 
Department of Physics, The State University of Iowa, Iowa City, Iowa 


(Received December 24, 1946) 


Measurements were made of the internal friction of 
single crystal rods when vibrating longitudinally at 45 
kilocycles per second. The decrement, or ratio of energy 
dissipated per cycle to twice the total vibrational energy, 
was used as a measure of the internal friction. At strain 
amplitudes below about 10-* or 10~* the decrement was 
independent of strain amplitude. The value of the constant 
decrement at low amplitudes, and the changes which 
occurred as the amplitude was increased, were both greatly 
dependent on the history of the specimen. The constant 
decrement for different crystals ranged from 710-5 to 
values higher than could be measured (about 2001075). 


For some of the initially high-decrement crystals, values 
decreased with time (room temperature anneal). In general, 
for annealed crystals the decrement, measured as a function 
of increasing strain amplitude, remained almost constant 
up to and beyond the static elastic limit and then rose 
rapidly. It returned to the initial value when the strain 
amplitude was again reduced. On a second run, however, 
the rise in decrement occurred at much lower amplitude 
than on the initial one. It is concluded that the decrement 
is the result of two separate processes. The changing 
decrement at higher strain amplitudes is qualitatively ex- 
plained with the use of dislocation theory. 





INTRODUCTION 


HE internal friction of single crystals of 

exceedingly pure zinc has been measured 
by Read,' and of less pure zinc by Read and 
Tyndall.? The decrement, which is the ratio of the 
energy dissipated per cycle to twice the total 
vibrational energy, was taken as a measure of the 
internal friction. Read found the decrement to be 
dependent on the past history of the specimen, 
strain or stress amplitude at which it was 
measured, time, etc. For the purer zinc decre- 
ments of the order of 10-5 were found at the 
lowest strain amplitudes but these, for most of 
the crystals, increased many fold with increasing 
amplitude. The present work consists of further 
observations of the behavior of the decrement of 
longitudinally oscillating zinc crystals under 
various circumstances and with partially con- 
trolled history. 


EXPERIMENTAL METHOD 


The zinc? used in this investigation was slightly 
less pure than that used by Read.' The crystals 


* The main contents of two dissertations submitted in 
partial fulfillment of the requirements for degrees of Doctor 
of Philosophy in the Department of Physics in the Gradu- 
ate College of the State University of Iowa. 

** April, 1943. 

*** Tune, 1946. 

t Now at the Naval Ordnance Test Station, Inyokern, 
California. 

1T. A, Read, Phys. Rev. 58, 371 (1940). 

2 T. A. Read and E. P. T. Tyndall, J. App. Phys. 17, 
713 (1946). ’ 

3 The zinc used was “Bunker Hill,” a commercial zinc 
with not more than 0.01 percent impurity. Read and 
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were grown by the Bridgman method in which a 
small casting is lowered through a furnace. The 
majority of the crystals were grown with the 
casting imbedded in dental investment com- 
pound and the rest were grown in the glass tube 
in which the original casting was made. The 
dental investment compound was easily softened 
by water after the crystals were grown and could 
then be removed with little danger of straining 
the crystals. The crystals had a cross-sectional 
area of 9 sq. mm and had fairly smooth surfaces. 

The decrements were measured by essentially 
the same method as that used by Read‘ in which a 
composite oscillator is constructed by cementing 
toa quartz crystal a zinc crystal of the same cross- 
sectional area and cut with a fine jeweler’s saw to 
a length such that the composite oscillator 
resonates at a frequency within less than 4 
percent of the resonant frequency of the quartz 
crystal. The resonant frequency of the quartz 
crystal used in the present work was 45 kilocycles 
when vibrating longitudinally in the fundamental 
mode. The composite oscillator was suspended by 
silk fibers placed at displacement nodes of vibra- 
tion between, but not touching, parallel plane 
electrodes. The assembly was then placed in a jar 
and the jar evacuated to a pressure of less than 


Tyndall (reference 2) report data from three crystals 
grown from this same zinc and two with additional im- 
purity. 

4 The writers are indebted to Dr. Read for information 
about his bridge circuit, etc., in addition to that in 
reference 1. 
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0.1 mm of mercury. The quartz crystal was 
etched with hydroflouric acid before using and 
had a decrement of 2.7 10~-°. 

The electrical apparatus consisted of an alter- 
nating-current bridge containing the composite 
oscillator as one arm and driven by a stable 
35-100 ke oscillator. The detector used with the 
bridge consisted of an r-f amplifier, pentagrid con- 
verter, and associated oscillator which changed 
the high frequency signal from the bridge to an 
audible signal, and several stages of audio ampli- 
fication. Its sensitivity with headphones was 
greater than that needed in any of the measure- 
ments. Either a vacuum thermocouple or a 
vacuum tube voltmeter were used to measure the 
voltage ¢ applied to the composite oscillator. The 
maximum strain amplitude, Uo, in the zinc 
specimen was then computed from the following 
formula :5 


2e 
Uy =——-(10K)!X 10%, (1) 
RV, 


where R is the electrical resistance of the com- 
posite oscillator at resonance, Vo is the velocity 
of sound in the zinc specimen, and K is the 
constant relating the inductance L of the com- 
posite oscillator and the masses M, and M; of the 
cylinders in the equation : 


L=K(M,+M)). (2) 


The subscripts 1 and 2 in Eqs. (1) and (2) refer to 
the quartz and specimen, respectively. 

The decrements of the zinc crystals were ob- 
tained from a knowledge of the inductance in the 
equivalent circuit, the resistance measured by 
the bridge, the resonant frequency of the com- 
posite oscillator, and the masses of the quartz and 
specimen. The constant K was evaluated with 
Eq. (11) in a paper by Cooke® and Eq. (2), from 
the observed variation of the reactance of the 
quartz crystal alone with frequency in the im- 
mediate neighborhood of the resonance fre- 
quency. It was determined to an accuracy of 1.5 
percent. The accurate measurement of changes in 





5 This formula is an obvious modification of formula 
(1), reference 2, for stress amplitude. It may be obtained 
by equating the expression for total energy in electrical 
terms to that in mechanical terms, and solving for Up. 

*W. T. Cooke, Phys. Rev. 50, 1158 (1936). Cooke’s K 
is eight times the K used in Eq. (2). 
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frequency required in the plotting of the reso- 
nance curve was obtained as follows. The crystal 
controlled multivibrators of a piezoelectric fre- 
quency standard were adjusted in frequency 
until a harmonic of one of them was within about 
50 cycles of the resonance frequency of the quartz 
crystal. The difference frequency obtained from 
this reference frequency and the frequency of the 
oscillator driving the bridge was amplified and 
used to drive a Cenco counter. The number of 
counts per second was then equal to the deviation 
of the unknown frequency from the crystal con- 
trolled reference frequency. 


EXPERIMENTAL RESULTS 


Before describing the quantitative results, 
some qualitative observations may be of interest. 
It was noted that some of the crystals could be 
made to give a high pitched “ringing” sound 
when set into flexural oscillations by suspending 
them from the mid-point and tapping lightly. 
Later work showed that crystals for which this 
ringing was audible for 2 or 3 seconds had 
decrements of the order of 10~* or less, whereas 
higher decrement crystals instead of ringing 
yielded a metallic or even wooden sound. A rough 
idea of the behavior of the decrement of a crystal 
could be obtained by subjecting it to this simple 
test. For example, eight crystals were tested just 
after they had been sawed to the desired length 
and then again after an anneal at 250°C for three 
hours. No noticeable difference was obtained at 
this time, but after storage at room temperature 
for several days a noticeable drop in the decre- 
ment, as judged by the increased time of ringing, 
was observed for some of the crystals while the 
rest showed little change. The crystals of the 
former class seemed to be persistent in their 
behavior. For example, crystal No. 1 (orientation 
15°) would always return to a low decrement 
after two or three days’ rest at room’ temperature 
without handling, regardless of its previous 
history. At one time when it was known to have a 
low decrement, it was accidently dropped onto 
the floor and bent considerably. It-was then 
straightened by hand and left idle for several 
days after which it appeared to ring as well as 
ever. The data in Fig. 1 (curve A) and Fig. 8 
(curve A) were taken on this crystal after it had 
rested nine days before it was mounted in the 
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apparatus and 73 hours afterward. Other crystals 
behaved similarily although none of these were 
subjected to such drastic treatment. Crystals 
which did not develop a ring in a few days after 
growth did not do so after months of lying idle. 
Several crystals of this type were also annealed at 
250°C and still did not develop a ring. These 
qualitative tests moreover showed very clearly 
how little handling was needed to increase the 
decrement. If a crystal which had a good ring was 
placed on a table and one end was elevated about 
a fourth of an inch and then dropped, the crystal 
thereafter would give only a metallic ‘clink’ 
when tapped. In fact, while the crystals were 
being subjected to the ringing test, the gentle 
taps needed to set up the oscillations a few times 
would noticeably reduce the time of ringing. One 
or two days of rest at room temperature after the 
above treatment was however sufficient to restore 
the ringing properties of the crystals. 

Nothing was found in the method of growth, 
amount of handling, etc., which gave any clue as 
to why some crystals had a low decrement while 
others had decrements which were too high to 
measure with the apparatus and which would not 
decrease with either a high or room temperature 
anneal. There was no apparent correlation of high 
and low decrement crystals with the degree of 
smoothness of their surfaces. 

The writers observed, as did Read,! that the 
decrement measured at low amplitudes immedi- 
ately after mounting in the apparatus was always 
much higher than after the crystal had been 
allowed to rest in the apparatus for some time 
undisturbed except for infrequent low amplitude 
measurements.’ Two examples of this behavior 
are given in Fig. 1. Curves A and B are for two 
different crystals. The decrease in decrement with 
time after mounting was caused by, at least 
primarily, the self-annealing of the crystals at 
room temperature. Part of the decrease, par- 
ticularly the initial rapid decrease, may have been 
caused by the drying of the cement between the 
specimen and the quartz. Calculations based on 
estimates of the thickness of the cement layer, 


7 The crystals of reference 2 made of the same zinc as 
the writers’ did not show this decrease of decrement with 
time. These crystals however were of larger cross section 
than the writers’ and the natural conclusion seems to be 
that they could be mounted without producing any 
detrimental strains. 
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TIME AFTER MOUNTING (HOURS) 


Fic. 1. Decrease of decrement with time after mounting, 
showing room temperature annealing. Curve A was taken 
on crystal No. 1. 


etc., indicate, however, that the decrement of the 
cement would have to be as large as 100 to 1000 
to cause the high decrements obtained. This 
seems improbable in view of the fact that the 
cement had dried for twenty or thirty minutes 
before the crystals were mounted. Also, Fig. 2 in 
Read’s paper! shows a similar decrease of the 
decrement with time after mounting and his 
crystals were fastened to the quartz with beeswax 
and rosin which would show no drying effect. 
Observations were made of the decrement as a 
function of maximum strain amplitude and time 
for thirteen crystals. It is not possible to present 
all of these data in full detail, but certain rather 
typical data are summarized and presented 
graphically. In Fig. 2 are shown curves of the 
decrement against strain amplitude for three 
crystals which had been allowed to rest in the 
apparatus for the following lengths of time: 
No. 2, two days; No. 3, three days; No. 4, thirty 
hours. In all cases the decrement had decreased 
to a nearly constant value before the data were 
taken. Crystal No. 4 (curve C) is particularly 
interesting as its decrement was practically con- 
stant to a value of maximum strain amplitude 
over five times that at the static elastic limit? 
(about 80X10~’) for this crystal. Subsequent to 
the data in Fig. 2, crystal No. 4 was stored at 
room temperature for over six months. It was 
then remounted in the apparatus and the data 
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Fic. 2. Typical decrement versus maximum strain ampli- 
tude curves. Curve A was for crystal No. 2 (@=71°); 
curve B was for crystal No. 3 (@=71°); and curve C was 
for crystal No. 4 (@=70°). 


shown in Fig. 3 were taken. The decrement had 
been allowed to reach a fairly constant value 
after mounting before the data were taken. The 
time given with each curve is the elapsed time 
after the taking of the first curve. The decrement 
was measured for increasing and then decreasing 
strain amplitude, as shown by the arrows. The 
decrement at the lowest strain amplitude de- 
creased from 95 X10-5 to 7010-5 in the course 
of 143 hours, but did not come down to the 
previously measured value of 7.0 10~® (Fig. 2). 
It may be seen from the lower branch of each of 
the curves shown in Fig. 3 that the strain 
amplitude required to increase the decrement by 
some amount, say 50 percent, increases con- 
siderably with elapsed time. In curve J the maxi- 
mum strain amplitude was at the static elastic 
limit (8010-7) before the decrement started to 
increase very rapidly. Where two curves were 
taken with little intervening time, such as B and 
C, G and H, and I and J, the ascending branch 
of the second curve of a pair was nearly identical 
to the descending branch of the previous curve. 
After these data were taken, the crystal was 
removed from the apparatus and allowed to rest 
for 20 days. It was then remounted in the appa- 
ratus and left idle for seven days after which the 
data shown in curve B, Fig. 4, were obtained. 
The decrement in this curve shows only a gradual 
rise to a strain amplitude nearly three times that 
at the elastic limit. The decrement at low strain 
amplitude is still however about nine times as 
large as the first previously measured value 
(curve C, Fig. 2). 
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Fic. 3. Decrement versus maximum strain amplitude Up 
for crystal No. 4 (@=70°). The time given by each curve 
is the elapsed time after the first curve was taken. These 
curves show the effect of self-annealing at room tem- 
perature. 


The curves shown in Fig. 5 demonstrate very 
well the phenomenon mentioned in connection 
with Fig. 3. These curves were taken on crystal 
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No. 6 (@=56°) which was known to be strained 
prior to mounting. Data were first taken to a 
strain amplitude of 5010-7 (curve 1) and the 
strain amplitude was then reduced to its lowest 
value. Curve 2 was then run toa strain amplitude 
of 8010-7, with curves 3 and 4 following in a 
similar fashion. From these curves it appears that 
previous runs affect subsequent runs only at 
strain amplitudes below the highest maximum 
strain amplitude obtained on the previous run. 
Thus, had the crystal been taken continuously 
from the minimum to the maximum strain 
amplitude (9410-7), it would probably have 
followed the lowest curve, which is a smooth 
continuation of curve 1. Also, while descending 
curves are not shown, they coincide quite well 
with the next ascending curve, as previously 
pointed out in connection with Fig. 3. 

The curves shown in Fig. 6 taken on crystal 
No. 5 (orientation 80°) illustrate some effects of 
room temperature annealing. The decrement of 
this crystal remained constant until a critical 
value of strain amplitude was obtained and then 
started to rise. As successive runs were taken, 
this critical value of strain amplitude increased. 
The data for curve A were taken soon after the 
crystal was mounted. Data for curve B were 
taken 53 hours after A was taken; and curve C, 
98 hours after B. Other sets of data were taken on 
the crystal during these two periods and there- 
fore they should not be considered strictly as rest 
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Fic. 4. Decrement versus maximum strain amplitude U» 
for crystals No. 9, 4, 10, and 5 corresponding to curves A, 
B, C, and D, respectively. All of these crystals had rested 
for a long period of time in the apparatus after mounting 
and before measurements were taken. The arrows on each 
curve designate the strain corresponding to the minimum 
static stress which will produce glide in the basal plane. 
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Fic. 5. Decrement versus maximum strain amplitude Uo 
for crystal No. 6 (@= 56°), which was known to be strained 
before mounting in the apparatus. The curves are shown 
for increasing strain amplitude only and were taken in 
the order designated. They show that the decrement was 
influenced by previous runs only at strain amplitudes 
lower than the highest strain amplitude obtained on 
previous runs. 
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Fic. 6. Decrement versus maximum strain amplitude U» 
for crystal No. 5 (@=80°). Curve A was taken soon after 
mounting, curve B 53 hours after A, and curve C 98 
hours after B. Succeeding curves were very little different 
from curve C. As the crystal annealed at room tempera- 
ture, the point of rise of the decrement occurred at higher 
strain amplitudes. 


periods. Finally the crystal refused to anneal any 
further at room temperature and data could be 
very nearly repeated on successive runs taken 
several days apart. These later curves were 
substantially identical to curve C, Fig. 6. At this 
point the crystal was annealed for several days 
at 260°C, mounted in the apparatus, and left idle 
for a week. The decrement then remained nearly 
constant for a range of stress amplitude 2} times 
the value at the static elastic limit, as shown by 
curve D, Fig. 4. It will be noticed that after the 
high temperature anneal and remounting, the 
decrement at low strain amplitudes was higher 
than before the anneal (Fig. 6). 

The two curves given in Fig. 7 indicate the 
diverse results that can be obtained when the 
conditions are very nearly the same. Curve A is 
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Fic. 7. Decrement versus maximum strain amplitude U>. 
Curve A was for crystal No. 7a and curve B was for crystal 
No. 7b. Each of these two crystals were sections cut from 
one long 60° crystal. 


for crystal No. 7a and curve B is for No. 7b. 
These two crystals were sections cut from the 
same original 60° crystal and annealed for two 
hours at 250°C. Crystal No. 7a was mounted four 
days after the anneal and curve A was taken 
after it had rested in the apparatus for one day. 
Crystal No. 7b was mounted eleven days after 
the anneal and curve B was taken after it had 
rested in the apparatus for two days. Curve A 
rises sharply at a low value of strain amplitude, 
while curve B shows a constant decrement as the 
strain amplitude is increased to a value of six or 
seven times that at which curve A starts to rise. 
The only known difference in the histories of the 
two crystals was in the time they were allowed to 
rest before and after mounting in the apparatus. 
Crystal No. 7b apparently had sufficient time to 
anneal at room temperature while No. 7a was 
still in a strained state. One other curve was 
taken on crystal No. 7b one day previous to curve 
B (one day after mounting), and it was nearly 
identical to curve B indicating that this crystal 
was in an annealed state (see discussion of 
Fig. 6). 

Curve A of Fig. 8 was taken on the crystal 
(crystal No. 1, orientation 15°) which was 
dropped, bent, and subsequently straightened, as 
mentioned earlier. After the above treatment it 
was left idle for five days and mounted in the 
apparatus. At this time its resonant frequency 
was found to be too low and it was necessary to 
cut it to a shorter length. Three days after the 
sawing operation it was remounted in the appa- 
ratus and left idle for about eight hours except 
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for occasional measurements of the decrement at 
very low strain amplitudes (curve A, Fig. 1). The 
data for curve A, Fig. 8 were then taken. 
Curve B represents data on crystal No. 8 
(orientation 73°). It had rested for several days 
before mounting, and 34 hours after, before the 
data were taken. These two crystals had about 
the same resolved shear stress in the basal plane 
for a given tensile stress. The two curves are 
similar in form, indicating that crystal No. 1 had 
recovered fairly well from the severe strain 
given it. ‘ 

The curves A, B, C, and D in Fig. 4 are for 
crystals which had been-allowed to rest for long 
periods of time in the apparatus after they were 
mounted. The histories of crystals No. 4 and 
No. 5 have been given earlier in this paper. 
Crystal No. 9 (orientation 50°) and crystal 
No. 10 (orientation 20°) were left idle in the 
apparatus for about the same periods of time as 
crystals No. 4and No. 5. The decrement of three 
of the crystals rose slowly with strain amplitude 
for values of strain amplitude considerably 
greater than the value of the static elastic limit, 
while the decrement of the other crystal began to 
rise sharply at about the static elastic limit (curve 
A). This crystal (orientation 50°) probably was 
still in a strained state as crystals with orienta- 
tions near 45° are much more easily strained than 
those with orientation angles near 0° or 90°. 

When the voltage across the composite oscil- 
lator was suddenly changed from one constant 
value to another, the decrement took an ap- 
preciable time to change from the previous value 





/{ 


+. 


_— 
Ne 


ar 





™— 
4 
ON 


w 
‘ 
—" 





a x10 

















: a 
5 
. ie) 


i} 
u, x 107 





Fic. 8. Decrement versus maximum strain amplitude Uo. 
Curve A was for crystal No. 1 and curve B was for crystal 
No. 8. 
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to the final value at the new stress amplitude. 
The curves in Fig. 9 illustrate this effect. For 
curves 1 and 2 the voltage change was an increase 
while for the others the voltage was decreased. 
Usually when the voltage across the crystal was 
changed, a time of less than 20 seconds was re- 
quired for the decrement to attain, say, 90 
percent of its total change. In two of the curves 
of Fig. 9, however, the decrement required about 
a minute to attain 90 percent of its change. No 
attempt was made to fit these data with a 
mathematical function because the data were too 
inaccurate, particularly in the region where the 
decrement was changing rapidly with time. In 
addition there were other complicating factors. 
The strain amplitude is proportional to the cur- 
rent through the crystal and hence as the 
decrement (proportional to resistance, which 
determines the current) changed with time after 
the voltage across the crystal was changed to the 
new value, the strain amplitude was not constant 
during the time indicated in Fig. 9. Also, because 
of the low decrements involved, the time re- 
quired for the strain amplitude to reach a steady 
state upon a change of the applied voltage from 
one value to another was large enough to make 
the effective zero of time uncertain. 


DISCUSSION OF DATA 


A number of points should be mentioned in 
regard to the measurements described above. 
These are: 

(1) The decrement measured at low amplitude 
decreased with time after mounting in the appa- 
ratus, approaching a constant value (Fig. 1). The 
final value reached was not, however, unique, in 
that a crystal after a second mounting would 
show a decline in decrement to a constant value 
but not necessarily to the same value found in the 
first experiment. 

(2) The decrement was independent of strain 
amplitude for all maximum strain amplitudes 
below a critical value of the order of 10-5 to 10~® 
(the bridge detector was sufficiently sensitive to 
balance the bridge at strain amplitudes much 
lower than 10~’, probably 10-*). The value of the 
low amplitude decrement was, for a given 
mounting in the apparatus, affected only by time 
(as shown in Fig. 1). 

(3) The value of the low amplitude decrement 
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TIME IN SECONDS 


Fic. 9. Curves showing for crystal No. 7a the change in 
decrement with time following a sudden change of driving 
voltage across the composite piezoelectric oscillator. The 
curves were taken in the order designated by the numbers 
on each curve. 


was higher for all crystals measured than for 
Read’s! spectroscopically pure zinc crystals. The 
lowest decrements found were about the same as 
for those of crystals made of the same zinc, as 
reported by Read and Tyndall.? It does not 
follow, however, that the decrement is a function 
only of the amount of impurities present. Other 
factors are much more important as is shown by 
the data on crystal No. 4. In Fig. 2 it had a 
decrement of 710-5 whereas after it had lain 
idle for six months and was remounted in the 
apparatus, the decrement at low amplitudes was 
70X10-5 (Fig. 3). 

(4) The value of strain amplitude at which the 
decrement started to increase became greater 
with elapsed time following mounting in the 
apparatus, i.e., increased as the value of the low 
amplitude decrement decreased. This effect is 
illustrated in Figs. 3 and 6. 

(5) If a crystal had been oscillated at a suffi- 
cient amplitude momentarily to increase the 
decrement, then on measurements made soon 
after, the decrement was larger at all strain 
amplitudes between the value at which the 
decrement started to rise and the highest value 
obtained on the previous run. The amount that a 
run influenced the decrement obtained on a suc- 
ceeding run decreased with elapsed time between 
the runs. The effect of time between runs is 
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illustrated particularly well by Fig. 5 but is also 
evident from the series of curves in Fig. 3. 

(6) The decrement lagged in time behind the 
oscillation amplitude when the amplitude was 
changed suddenly. It is apparent from these 
curves that a finite time was required for the 
to reach a stable value after the 
oscillation amplitude was changed. 


decrement 


The shape of the curves of decrement versus 
strain amplitude suggests that the observed 
decrement was composed of a constant decre- 
ment plus an amount which varied with strain 
amplitude. The portion of the decrement (to be 
designated by decrement A) which was constant 
with strain amplitude was measured alone at 
small strain amplitudes. Decrement B is that 
portion which was zero or negligible at zero 
strain amplitude and increased with increase of 
strain amplitude. It is possible on this basis to 
give a qualitative explanation of the behavior of 
the decrement in the experiments described 
above. The dislocation theory is used together 
with several assumptions to be listed later. 

The dislocation theory’ is based upon a certain 
type of local imperfection in the crystal lattice 
called a line dislocation. The local imperfection is 
concentrated about a line in a slip plane and is 
composed of such an arrangement of atoms that 
movement of the line imperfection through the 
crystal along a slip plane gives the same effect as 
if the material above the slip plane had been 
shifted a distance of one atomic spacing with 
respect to that below. The dislocations can be 
either positive or negative, one being comple- 
mentary to the other in the sense that it is the 
mirror image of the other. A dislocation experi- 
ences a force in a shear stress field which reverses 
in sign with change in sign of the dislocation. 
Hence since a dislocation sets up a shear stress 
field of its own, dislocations of unlike sign are 
attracted to each other and those of like sign are 
repelled. Also a dislocation is attracted to the 
crystal surface by the force which would be given 
’ by an image dislocation of opposite sign outside 
of the surface. Dislocations can be formed in 


’ A summary of the dislocation theory is given in the 
article by Frederick Seitz and T. A. Read, J. App. Phys. 
12, 100 (1941). The forces between certain types of dis- 
locations and other related matters are discussed in J. S. 
Koehler, Phys. Rev. 60, 397 (1941). 
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pairs in the interior of a crystal or singly at the 
surface as the result of shear stress. 

In addition to the properties of dislocations 
summarized briefly above, the following as- 
sumptions will be used in the interpretation to be 
given for the results: 

(a) The dislocations present as a result of 
previous low frequency, or “‘static”’ strains tend 
to form a stable lattice-like array within the 
crystal specimen.$ 

(b) The low amplitude decrement (A) is a 
function only of the total number of dislocations 
present at a given time in the specimen. 

(c) The number of dislocations, if any, that are 
formed in the specimen as a result of the oscil- 
lating stress is a negligible fraction of the total 
number present. Presumably the short time 
during which stress was applied did not allow 
much opportunity for dislocations to form either 
pairwise in the interior, or at the surface. 

(d) The value of decrement B is a function of 
the number of dislocations present which are 
“free’”’ to dissipate energy. Whenever a disloca- 
tion is held within a dislocation lattice it is 
assumed to mean that the dislocation is pre- 
vented, or at least restricted, from participation 
in the energy dissipation process. 

A crystal which had been recently mounted in 
the measuring apparatus would be expected to be 
in a strained state as a result of handling during 
mounting. Hence there would be a large number 
of dislocations of both signs distributed through- 
out the crystal. With the aid of thermal fluctua- 
tions, the forces on and between dislocations 
would cause them to migrate. Some dislocations, 
particularly those near the surface, would diffuse 
out of the crystal, others would disappear by 
cancellation of positive and negative dislocations 
(if the pair was originally in the same slip plane), 
and the remaining ones would migrate to posi- 
tions such that the crystal had a minimum of 
potential energy, i.e. they would arrange them- 
selves in some sort of orderly pattern. After a 
long time the number of dislocations in the 
crystal would be essentially constant, all being 
held in the lattice. Under assumption (b), that 
decrement A is a function only of the total 
number of dislocations present, the curves of 
Fig. 1 would follow. Furthermore, the final value 
reached by the decrement several days after 
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mounting would not be expected to be the same 
for a different mounting in the apparatus. Each 
mounting would entail a different kind and 
amount of the unavoidable strain, resulting in a 
different number and distribution of dislocations 
immediately thereafter. Hence the stable lattice 
reached would not be unique. This explains such 
data as that in Fig. 2 (crystal No. 4) where the 
low amplitude decrement was 7 X10~° and later 
after the crystal was remounted the decrement 
was 7010-5, as shown in Fig. 3. The process by 
which energy is dissipated in an amount pro- 
portional to the total number of dislocations and 
independent of amplitude is difficult to determine 
without more detailed information. One possi- 
bility is that irreversible heat flow occurs as the 
result of oscillating stress in an elastically 
inhomogeneous medium. In this event, the dislo- 
cations would presumably provide the elastic 
inhomogeneity. Quantative calculations should 
show whether or not this possibility is reasonable. 

Decrement B, a function of strain amplitude, 
could be interpreted according to the following 
picture. The dissipation represented by decre- 
ment B would be assumed to be the result of 
movement of the dislocations. This movement 
would result in a plastic change in length of the 
specimen, with consequent loss of oscillation 
energy. The data could then be interpreted by 
assuming that the dislocations did not move 
appreciably when they were attached to the 
dislocation lattice and that they could be broken 
away from the lattice by the application of a 
sufficiently large stress. The amount of stress 
required would vary according to how tightly 
each dislocation was bound in the lattice. Hence 
as the oscillation amplitude was increased 
steadily from zero, more and more dislocations 
would be freed from the restraining influence of 
their neighbors and be made available to dissi- 
pate energy, with a consequent rise of the 
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decrement. The decrement, during the subsequent 
decrease of amplitude would then be larger at a 
given amplitude than on the ascending curve 
because more dislocations would be free than 
before. Examples of this hysteresis-like effect are 
given in Figs. 2, 3, 5, 7, and 8. This picture, 
furthermore, would explain the effect shown in 
Figs. 3 and 5 in which an ascending branch of a 
curve was nearly identical to the descending 
branch of a curve taken soon before it. It would 
be due to the fact that the same number of 
dislocations were free to dissipate energy on both 
curves. The dislocations that had been freed from 
the dislocation lattice by the oscillating stress 
would, upon removal of the stress, gradually take 
either their old positions, or new positions in the 
lattice. This process would, however, take a 
finite time as thermal energies would be expected 
to be an important factor. The data given in 
Fig. 3 show how elapsed time after a run erases 
the effect of a previous run. Also, the increase of 
the point of rise of the decrement shown in Fig. 6 
would be caused by the fact that as the crystal 
had more time to anneal, the dislocations which 
were loosely bound in the lattice would disappear 
from the crystal, leaving only those which were 
tightly bound and required a large stress to make 
them available to dissipate energy. 

The number of dislocations that are freed from 
the lattice would be expected to depend some- 
what on the time that the crystal was oscillated 
at a given amplitude, i.e., some time would be 
required to establish an equilibrium state. The 
curves of Fig. 9 show that several seconds were 
required for a stable value of the decrement to be 
reached when the oscillation amplitude was 
changed from one constant value to another. 

In conclusion, the writers wish to express their 
appreciation to Professor E. P. T. Tyndall for his 
guidance and help during the progress of this 
work. 


425 








Letters to the Editor 








Measurement of Surface Tension 


ROBERT WEIL 


South-West Essex Technical College and School of Art, 
Walthamstow, England 


January 27, 1947 


UGDEN'’S! method was employed for measuring the 
surface tension of distilled water at a temperature of 
25°C in order to determine how the value would be affected 
by imperfections in the circularity of the capillaries. The 
procedure was as follows: six capillaries were drawn out 
from pieces of glass tubing so that the diameter at the tip 
was between .04 and .06 cm. While still soft the tips were 
squeezed to different degrees; thus the cross sections as- 
sumed elliptical shapes. It was argued that the ratio R of 
the minimum to the maximum diameter was a measure of 
the deviation from the true circularity. This ratio, as 
determined by means of a cathetometer, varied from .537 to 
.914. The diameter of the large tube was .2814 cm. 


The maximum bubble pressures were determined for’ 


each specimen. A high rate of bubbling was used at first 
with smaller rates following it. The mean of three readings 
for each rate was taken: the divergences were negligible. 
No automatic counting device being available, it was not 
possible to exceed a rate of about 140 to 150 bubbles per 
minute. Graphs were drawn plotting pressure against rate 
and extrapolated to infinitely long periods of bubble 
formation. This extrapolation seemed justified since, after 
some practice, it became possible to reduce the rates to as 
small values as two bubbles per minute. A mean value of 
the measured radii was used in order to calculate the 
surface tension. 

Two interesting results have been found: for .9<R<1 
values for the surface tension agree very well with data as 
obtained by other workers with capillaries whose cross 
section approximates to a circle to a greater degree than 
those used in this investigation. 

It would appear that Martin’s* conclusions obtained from 
work on similar deviations in Jager’s original method do not 
apply here. 

Secondly, when plotting the pressure p against the rate 
of bubbling, it was found in all cases that the values of p 
increased with the rate of bubbling up to about 60 bubbles 
per minute or more, and then definitely decreased for 
higher rates. It was found, as might be expected, that, on 
raising the rate to a greater value than could be measured 
visually, the pressure rose considerably. The cause for this 
kink in the curves is not quite obvious. It is a matter of 
_ some conjecture whether it can be accounted for by the fact 
that at small rates actual maxima and minima are recorded 
in the positions of the liquid in the gauge, while at great 
frequencies it is the mean positions which are observed. In 
some unpublished work L. B. Wood points out the necessity 
of extrapolating to infinitely long periods of bubble 
formation to obtain reliable and consistent results. It would 
seem that there are occasions when it is not sufficient to 
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maintain the rate of bubbling at less than two per second, 
as is suggested by Adam. 

It was not possible to extend the measurements to 
capillaries with different diameters and different wall 
thicknesses: it is certainly desirable that this should be 
done. 

' Sugden, J. Chem. Soc. p. 858 (1922); p. 27 (1924). 


2 Martin, Sitz. Akad. Wiss. Wien 123, 2491 (1914). 
% Adam, The Physics and Chemistry of Surfaces, p. 375. 


Letter to the Editor 


FRANCIS J. MURRAY 


The Institute for Advanced Study, School of Mathematics, 
Princeton, New Jersey 


January 8, 1947 


AY I respectfully call your attention to the fact that 
the suggestion of Mr. Eaglesfield on linear equation 
solvers, on page 1125 of the current issue of the Journal of 
Applied Physics (December, 1946) duplicates a paper of 
mine presented to the American Mathematical Society, 
September 1945, and abstracted in the November issue of 
the Bulletin of the American Mathematical Society, Vol. 51 
(1945), page 883. 

The content of the paper appears in full in my book 7he 
Theory of Mathematical Machines, a current publication of 
the Columbia University Press. 

In this paper, I prove the convergence in the case in 
which = d,? is used. In the case of a machine which uses 
maxd; as an indication, one can give examples where the 
machine will reach a minimum not zero. 

I have constructed a number of models upon the = d,? 
principle and a large one is now being constructed by the 
Watson Scientific Computing Laboratories at Columbia 
University. 





Here and There 








Personnel 


The National Bureau of Standards recently announced 
the retirement of two longtime employees. Harvey L. 
Curtis, physicist at the Buread since 1907, had deferred his 
retirement in 1945 to continue direction of interior ballistics 
research. Frederick J. Bates, sugar physicist who had also 
been chief of the Optics Division since 1941, has given more 
than 43 years of continuous service to the Bureau. 


On December 30 Carroll L. Wilson resigned as vice 
president and director of National Research Corporation in 
order to assume his duties as general manager of the 
Atomic Energy Commission. 


Ross Gunn has resigned as superintendent of the Physics 
Division, Naval Research Laboratory, to act as director of 


physical research for the U. S. Weather Bureau. 


In January three appointments to the new Department 
of Engineering Sciences and Applied Physics at Harvard 


JOURNAL OF APPLIED PHYSICS 





d 


ie 


of 


it 





University were announced. Leon N. Brillouin was ap- 
pointed Gordon McKay Professor of Applied Mathematics, 
and James B. Fisk and Ronold W. P. King were named 
Gordon McKay Professors of Applied Physics. 


Chester Peterson has been appointed physicist in the 
Resistance Measurements Section of the National Bureau 
of Standards. 


The following men have new appointments to the De- 
partment of Physics at the University of Minnesota: 
Charles A, Critchfield, Frank Oppenheimer, Edward J. 
Lofgren, Joseph Weinberg, and Otto H. Schmitt. 


Division of High-Polymer Physics 


The following are officers for 1947 of the Division of 
High-Polymer Physics of the American Physical Society: 
Chairman, Lawrence A. Wood; Vice Chairman, Hubert M. 
James; Secretary, W. James Lyons; Treasurer, Robert S. 
Spencer. Professor James and Dr. Spencer were also elected 
to the executive committee. 


General Electric Fellowships at Union College 


For the third year General Electric Company is offering 
50 six-week, all-expense fellowships to high school science 
teachers in seventeen northeastern states, for summer 
study at Union College in Schenectady, New York. 


Leonard Loeb Honored 


School and Society reports that Leonard B. Loeb, pro- 
fessor of physics, University of California (Berkeley), has 
been chosen to deliver the 38th annual Kelvin lecture at the 
Institution of Electrical Engineers, London, England, on 
April 24. Professor Loeb will speak on “Electrical Dis- 
charge through Gases.” 


Junior Chamber Puts Physicist in Top List 


The U. S. Junior Chamber of Commerce published in 
January its selection of the nation’s ten outstanding young 
men of the year. Included among them was Philip 
Morrison of Cornell University, atomic physicist, “for his 
outstanding work in the development of the atomic bomb 
and his struggle to educate the people to the danger and 
promise of atomic energy.” 


Electronics Fellowships at MIT 


A number of graduate and advanced research fellowships 
are being offered by the Massachusetts Institute of 
Technology, Cambridge, Massachusetts, for study and re- 
search in the field of electronics. They wili be known as 
Industrial Fellowships in Electronics and are sponsored 
jointly by a group of industrial organizations concerned 
with the advancement of electronics and its applications. 
Applicants should communicate with the Director, Re- 
search Laboratory of Electronics. Under normal circum- 
stances application must be made at least four months 
prior to the intended date of entrance. 


VOLUME 18, APRIL, 1947 


Work to Start Soon on Schenectady Atomic Laboratory 


Construction of the Knolls Atomic Power Laboratory to 
be built near Schenectady, New York, by the Atomic 
Energy Commission will start this spring. The laboratory 
is being erected under the supervision of the General 
Electric Company, which will operate it for the government 
when completed. It is predicted the new building will be 
occupied by the middle of 1948, 


Acta Crystallographica 


A special international organization, formed largely for 
this purpose, is planning the publication of a journal to 
perform the services formerly rendered by the Zeitschrift 
fiir Kristallographie. Acta Crystallographica, as it will be 
called, will be published either in England or America, and 
while French, German, and Russian will be permitted 
languages, the bulk of the journal will be in English. It is 
expected that there will be six issues per year totaling ap- 
proximately 1000 pages. The price to individual subscribers 
is tentatively set at $10 a year. 

The American representatives in the international 
organization are extremely anxious to obtain a good idea 
of the probable circulation of the journal in North America. 
Any reader of this notice who would, or who might, 
subscribe to Acta Crystallographica is urged to send a note 
or postcard to that effect addressed to Dr. Henry A. 
Barton, American Institute of Physics, 57 East 55 Street, 
New York 22, New York. 


Reorganization of American Institute of Physics 


Reorganization of the American Institute of Physics and 
plans for a new semi-popular journal devoted to physics 
and its relation to society have been announced by the 
Institute. Effective immediately, the 7000 members of five 
societies (The American Physical Society, The Optical 
Society of America, The Acoustical Society of America, 
The American Association of Physics Teachers, and The 
Society of Rheology) will become members of the Institute. 
Heretofore the Institute has been an organization of these 
five Member Societies without any individual memberships. 
Provision is also made for associate membership in the 
Institute, open to others interested in physics. No change 
is made in the organization, functions, and activities of the 
five societies. 

The new journal, for which initial financing is now being 
sought, will contain news about physicists, their meetings, 
their work, and other activities, and about similar features 
in related fields; information about the action of govern- 
mental and other non-physics agencies affecting physics; 
articles and letters presenting the views of physicists on 
problems of general importance; references to current 
literature in physics; book reviews, etc. Technical journals 
in the field of physics already being published by the 
Institute and its Member Societies will remain unchanged. 


Division of Electron and Ion Optics 


The Division of Electron and Ion Optics will hold a 
meeting in connection with the Washington meeting of the 
American Physical Society on May 1-3, 1947. There will be 
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a program of invited papers in the general field of “The 
Physics of Surface Phenomena.” Contributed papers will be 
welcomed. For further information address J. R. Pierce 
(secretary-treasurer), Bell Telephone Laboratories, Inc., 
463 West Street, New York 14, New York. 


The Bulletin of Mathematical Biophysics 


The following is the table of contents for the March 1947 
issue of The Bulletin of Mathematical Biophysics: 
Suggestions for a Mathematical Biophysics of Some Psychoses 

N. RASHEVSKY 
A Problem in Mathematical Biophysics of Interaction of Two or More 

Individuals Which May Be of Interest in Mathematical Sociology 

N. RASHEVSKY 
Mathematical Theory of Motivation Interactions of Two Individuals: 

I.—ANATOL RAPOPORT , 

The Mechanism of the Middle Ear: I. The Two Piston Problem 

MArtTINuS H. M, Esser 


The University of Chicago Press, Chicago, Illinois 
Volume 9 Number 1. 


Symposium at Ohio State 


A Symposium on Molecular Structure and Spectroscopy 
will be held at the Mendenhall Laboratory of Physics at 
The Ohio State University from Monday, June 9, through 
Saturday, June 14. There will be discussions of the in- 
terpretation of molecular spectroscopic data as well as 
methods of obtaining such data. In addition, there will be 
sessions devoted to microwave and Raman spectroscopy. A 
dormitory will be available for those who wish to reside on 
the campus during the meeting. For further information, or 
for a copy of the program when it becomes available, write 
to Professor Harald H. Nielsen, Mendenhall Laboratory of 
Physics, The Ohio State University, Columbus 10, Ohio. 


International Nickel Cooperative Educational Program 


A broadening of International Nickel Company’s co- 
operation with universities and colleges in the United 
States and Canada in the field of engineering education 
through the distribution of technical literature has been 
announced by T. H. Wickenden, manager of the develop- 
ment and research division of The International Nickel 
Company, Inc. 

The program will make available useful material for 
classroom instruction in training students in scientific 
fields. It has been offered to, and accepted by, a number of 
institutions in the United States and Canada which give 


‘accredited courses in mining, metallurgy, chemical en- 


gineering, and one or more other engineering courses. As 
rapidly as possible the program will be offered to all 
engineering schools. 

The development and research division of the company 
will furnish to each institution an exhibit containing ap- 
proximately 50 specimens of nickel-containing materials; 
a portable metals identification kit containing approxi- 
mately 35 specimens of important metals and alloys for 
qualitative identification of metals and alloys; literature 
and other information concerning nickel and its alloys; and 
other data. Motion pictures of the company’s mining, 
smelting, and refining operations will be made available. 
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New Books 








Circuit Analysis by Laboratory Methods 


By Cari E. SKRODER AND M. STANLEY HELM. Pp. 
288+xvi, 6X9 in. Prentice-Hall, Inc., New York, 
1946. Price $5.35. 


Some twenty years ago Professor Harry E. Clifford 
planned a system of laboratory reports which was used 
under his cognizance at M.I.T. and Harvard, and the form 
of which is still used there and in many other schools. The 
purpose of this system was to encourage the student to plan 
and think about the experiment rather than treat it as a 
mere manual exercise. The modus operandi was to require 
a preliminary report which involved the following sections: 
a. A statement of the problem, b. The theory involved, 
c. The procedure, d. The circuit diagrams, e. The results 
expected together with a sample calculation, and f. A 
proposed analysis of the results. 

After having this preliminary report, together with a list 
of needed instruments, accepted by the instructor, the 
student was then free to perform the experiment and was 
expected to perform it with little additional guidance. The 
experiment was finished only when a final report, including 
the completion of Sections e and f, was accepted. This 
method was and still is in sharp contrast with the cook book 
method prevalent in some schools. This latter method is 
sometimes reduced to a point where the student is asked 
only to read instruments and to put the readings in 
appropriate squares on a data sheet provided by the 
instructor. There is no question as to which method does 
the better job of teaching. 

In the book Circuit Analysis by Laboratory Methods, the 
authors attempt to follow the preliminary report system 
and at the same time reduce some of the drudgery involved 
in merely copying sections of theory from references. This 
they do admirably. However, by the title one is rather led 
to expect that the book will contain an entirely new and 
different set of experiments. This is not so, since many of 
the experiments are good old standbys, but the emphasis 
has very definitely been shifted from machines to circuits. 
Only three out of the thirty-five experiments are on 
rotating machine characteristics. Eight of the remainder 
are on series and parallel circuits including resonance, which 
acknowledges their importance both in engineering and as 
a teaching mechanism. Three are on polyphase circuits, and 
three on reactive power. The remainder are distributed 
among instruments, Kirchhoff’s laws, power, resistance, 
etc. 

Although the text explaining the philosophy of the 
experiments is good, it is expected that the course would be 
taken concurrently with a theory course in circuits. Since 
there should be some correlation between this and the 
laboratory, one would hope to find a classroom text which 
covers the same material in the same order. This might not 
be easy, although it might be accomplished by a certain 
amount of reorganization. 

The chapter describing the report to be expected of the 
student is good, but it could have been rendered more vivid 
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by including an example of what the authors consider to be 
a good report on one of the experiments included in the 
later chapters. 

It seems that the authors have lost some crispness in 
trying to make the book so general that it can be used in all 
laboratories with all kinds of equipment. Somehow the 
sentence, “The inductance of the coil and the capacitance 
of the condenser should be such that resonance will be 
obtained at a frequency about in the middle of the fre- 
quency range designated by the instructor,”’ does not seem 
nearly as authoritative as, “Choose the inductance and 
capacitance so that resonance will be obtained at 2000 cps.” 

Eric A. WALKER 
The Pennsylvania State College 


The Diffraction of X-Rays and Electrons by 
Free Molecules 


By M. H. Pirenne. Pp. 160+-xii, Figs. 82, 548} in. 
‘Cambridge Series of Physical Chemistry. Cambridge 
University Press, Cambridge, England, 1946. 


The thirteen chapters cover the following topics: 
Coherent scattering, incoherent scattering, fluorescence 
radiation, experimental measurements of the scattering of 
x-rays by atoms, diffraction of fast electrons by atoms, 
diffraction of x-rays by free molecules, intermolecular 
interferences in gases and liquids, thermal movements in 
molecules, the chemical bond, use of x-ray and electron 
diffraction for determination of structure, Fourier analysis 
of interference measurements of free molecules, measure- 
ment of x-ray scattering functions of gases, and a list of 
molecules studied by x-ray diffraction. An excellent 
bibliography contains 187 references. 

The book is essentially a report of the work of the Debye 
school in the field of x-ray scattering by free molecules. The 
material is well presented and the book makes easy and 
interesting reading. Theoretical considerations and deriva- 
tions are summarized, the results and conclusions are 
clearly stated, and the step-by-step details omitted. The 
first twelve chapters could well serve as a complete coverage 
for twelve semi-popular lectures. Most of the work pre- 
sented in the book is more than ten years old, and much of 
it is familiar to workers in the field. However, even the 
experienced x-ray diffraction worker will find that the 
careful authoritative presentation clears up for him a few 
misunderstood points. The book should serve a useful 
purpose as a general report and reference for this phase of 
the x-ray diffraction field. 

B. E. WARREN 
Massachusetts Institute of Technology 





New Booklets 





—_—- 





Tracerlab, Inc., 55 Oliver Street, Boston 10, Massa- 
chusetts, has published the first issue of Tracerlog, the 
purpose of which is “to publish technical information on 
laboratory, medical, and industrial uses of radioactivity, 
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not to mention an occasional ‘commercial’ plug for Tracer- 
lab’s own products and services.”’ It will be published 
monthly and mailed free of charge on request to physicists, 
chemists, physicians, and engineers interested in radio- 
activity and its applications. 


Tube Department of Radio Corporation of America has 
announced the publication of two 16-page booklets. Form 
PG-101 is entitled ‘‘Power and Gas Tubes for Radio and 
for Industry” and Form 1275-C is named “‘Receiving Tubes 
for Television, FM, and Standard Broadcast.” Either book- 
let may be obtained for 10 cents a copy from RCA tube 
distributors, or by sending 10 cents direct to Commercial 
Engineering, Tube Department, Radio Corporation of 
America, Harrison, New Jersey. 


Miniature Precision Bearings, Inc., Keene, New Hamp- 
shire, offers a four-page bulletin on miniature ball bearings 
for precision instruments and mechanisms. It describes in 
detail five series: radial, super-light radial, pivot, angular 
contact, and thrust. 


Acme Scientific Company, 1448 West Randolph Street, 
Chicago 7, Illinois, has issued a four-page folder describing 
in detail Acme’s new mold polishing service. The company 
feels that the type of surface which it can provide on 
irregular contours may be useful to the readers of this 
journal, 


Eberbach and Son Company, Ann Arbor, Michigan, 
publishes a 16-page house organ called Announcer of 
Scientific Equipment, of interest to research, control, and 
development laboratories. 


Rex Rheostat Company, P. O. Box 232, Baldwin, Long 
Island, New York, announces its new Catalogue No. 3, 
8 pages. The company manufactures slide-contact rheostats 
and resistors. 


Leeds and Northrup Company, 4934 Stenton Avenue, 
Philadelphia 44, Pennsylvania, has published a new 
catalog, ND46(1), which provides information about a 
pyrometer which supplements the Micromax line of 
instruments for those applications which require split- 
second response to temperature changes, or unusually rapid 
concentration on a single chart of data from many points. 
16 pages. 


Walker-Jimieson, Inc.,. 311 South Western Avenue, 
Chicago 12, Illinois, makers of radio and electronic supplies, 
have issued their 1947 Reference Book and Buyer's Guide. 
106 pages. 


Radio City Products Company, Inc., 127 West 26 Street, 
New York 1, New York, is announcing its new Catalogue 
No. 129, which describes a comprehensive line of radio, 
electrical, and electronic test instruments. 24 pages. 


American Documentation Institute, 1719 N Street, 
N.W., Washington 6, D. C., has published its first Catalogue 
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of Auxiliary Publications in Microfilms and Photoprints. To 
be distributed free so long as this first edition is available. 
52 pages. 


Dow Corning Corporation, Midland, Michigan, has 
announced the third edition of its DC Silicone Catalogue. 
12 pages. 


RCA Victor Division, Camden, New Jersey, has pub- 
lished two new booklets for the wider dissemination of 
technical information about electron tubes. Each contains 
16 pages, and the titles are Receiving Tubes for Television, 
FM, and Standard Broadcast, and Power and Gas Tubes for 
Radio and for Industry. 


Castle Films Division of United World Films, Inc., 30 
Rockefeller Plaza, New York 20, New York, is announcing 
a new catalog listing many new films released during the 
past year for school and industrial training purposes. 


The Ruberoid Company, 500 Fifth Avenue, New York 
18, New York, has issued a 40-page illustrated booklet by 
Dr. Oliver Bowles entitled A sbestos—The Silk of the Mineral 
Kingdom. Dr. Bowles is chief of the Nonmetal Economics 
Division of the U. S. Bureau of Mines and is a prominent 
authority on the subject of asbestos. 


Philips Research Laboratory has published a 20-page 
index listing all articles which have appeared in Philips 
Technical Review during the period January 1935 (first 
issue) to June 1942 (last issue prior to suspension of 
publication during the German occupation of Holland). It 
is available on request from Elsevjer Book Company, Inc., 
215 Fourth Avenue, New York 3, New York. The index, 
distributed free as a service to readers, alphabetically 
catalogs by subject all articles which appeared during that 
period. 
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